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Abstract 

In financial markets, the information that traders have about an asset is reflected in its 
price. The arrival of new information then leads to price changes. The 'information- 
based framework' of Brody, Hughston and Macrina (BHM) isolates the emergence of 
information, and examines its role as a driver of price dynamics. This approach has 
led to the development of new models that capture a broad range of price behaviour. 
This thesis extends the work of BHM by introducing a wider class of processes for the 
generation of the market filtration. In the BHM framework, each asset is associated 
with a collection of random cash flows. The asset price is the sum of the discounted 
expectations of the cash flows. Expectations are taken with respect (i) an appropriate 
measure, and (ii) the filtration generated by a set of so-called information processes that 
carry noisy or imperfect market information about the cash flows. To model the flow 
of information, we introduce a class of processes termed Levy random bridges (LRBs), 
generalising the Brownian and gamma information processes of BHM. Conditioned on 
its terminal value, an LRB is identical in law to a Levy bridge. We consider in detail 
the case where the asset generates a single cash flow X T at a fixed date T. The flow 
of information about X? is modelled by an LRB with random terminal value X?- 
An explicit expression for the price process is found by working out the discounted 
conditional expectation of X T with respect to the natural filtration of the LRB. New 
models are constructed using information processes related to the Poisson process, the 
Cauchy process, the stable-1/2 subordinator, the variance-gamma process, and the 
normal inverse-Gaussian process. These are applied to the valuation of credit-risky 
bonds, vanilla and exotic options, and non-life insurance liabilities. 



3 



Acknowledgements 

I am very grateful to Lane P. Hughston, my supervisor, for his help and support. 
With his breadth of knowledge, Lane's teachings have stretched beyond mathematical 
finance to such eclectic subjects as quantum mechanics, Italian opera, and Indian 
literature. Many thanks go my co-author Andrea Macrina whose enthusiasm for this 
work often pushed me forward. I would like to express my gratitude to members 
of the mathematical finance groups at Imperial College London and King's College 
London, and members of the London Graduate School in Mathematical Finance, with 
particular thanks going to my teachers M. Davis, H. Geman, W. Shaw, and M. Zervos; 
and also to P. Barrieu, D. Brody, M. Pistorius, and R. Norberg. I am grateful to 
F. Delbaen, and M. Schweizer, and other members of the mathematical finance group 
at ETH Zurich for their hospitality and support. I also wish to thank D. Taylor for 
his hospitality at the University of the Witwatersrand. This work was supported by 
an EPSRC Doctoral Training Grant and a European Science Foundation research visit 
grant under the Advanced Mathematical Methods in Finance programme (AMaMeF). 
Part of this work has undertaken when I was student at King's College London. Finally, 
for their unfailing support throughout my studies, I must thank my family — Mum, Dad, 
Charlotte and Thomas — and Christine. 



4 



The work presented in this thesis is my own. 



A.E.V. Hoyle 



Contents 



1 Introduction and summary 

2 Levy processes and Levy bridges 

2.1 Levy processes 

2.2 Stable processes 

2.3 Levy bridges 

2.4 Stable bridges 

2.5 Brownian motion and Brownian bridge 

2.6 Gamma process and gamma bridge 

2.7 Variance-gamma process and variance-gamma bridge 

2.8 Stable-1/2 subordinator and stable-1/2 bridge 

2.9 Cauchy process and Cauchy bridge 

2.10 Normal inverse-Gaussian process and normal inverse-Gaussian bridge 

2.11 Poisson process and Poisson bridge 

3 Levy random bridges 

3.1 Denning LRBs 

3.2 Finite-dimensional distributions 

3.3 LRBs as conditioned Levy processes 

3.4 The Markov property 

3.5 Conditional terminal distributions 

3.6 Measure changes 

3.7 Dynamic consistency 

3.8 Increments of LRBs 



9 



5 



4 Information-based asset pricing 

4.1 BHM framework 

4.2 Levy bridge information 

4.3 European option pricing 

4.4 Binary bond 

5 Variance-gamma information 

5.1 Variance- gamma random bridge .... 

5.2 Simulation 

5.3 VG equity model 

5.4 VG binary bond 

6 Stable-1/2 information 

6.1 Stable-1/2 random bridge 

6.2 Insurance model 

6.3 Estimating the ultimate loss 

6.4 The paid-claims process 

6.5 Reinsurance 

6.6 Tail behaviour 

6.7 Generalized inverse-Gaussian prior . . 

6.8 Exposure adjustment 

6.9 Simulation 

6.10 Multiple lines of business 

7 Cauchy information 

7.1 Cauchy random bridges 

7.2 Simulation 

7.3 Cauchy binary bond 

8 Normal inverse-Gaussian information 

8.1 Normal inverse-Gaussian random bridge 

8.2 Simulation 

8.3 NIG equity model 

8.4 NIG binary bond 



Contents 



7 



9 Poisson information Il22 

9.1 Poisson random bridge 123 

9.2 Mixed Poisson processes 132 

9.3 Simulation 134 

9.4 Compound Poisson random bridge 135 

9.5 rath-to-default baskets Il37 

A Some integrals Il39 

A.l Proof of Proposition [2X2 139 

A.2 Proof of Proposition l2iT5l 140 

A.3 Proof of Proposition EO 142 

A.4 Proof of Proposition \U4 



Bibliography 



146 



List of Figures 



2.1 Stable- 1/2 bridge simulations 

5.1 VG binary bond simulations: defaulting bonds 

5.2 VG binary bond simulations: non-defaulting bonds 

5.3 VG binary bond simulations: varying the rate parameter 

6.1 Craighead curves: truncated Weibull time-changes 

6.2 Simulations from the stable- 1/2 random bridge reserving model . . . . 

7.1 The Cauchy binary bond price plotted as a function of the information 
process 

7.2 Cauchy binary bond simulations: defaulting bonds 

7.3 Cauchy binary bond simulations: non-defaulting bonds 

8.1 NIG binary bond simulations: defaulting bonds 

8.2 NIG binary bond simulations: non-defaulting bonds 

8.3 NIG binary bond simulations: varying the rate parameter 

9.1 Poisson random bridge simulations 



8 



Chapter 1 

Introduction and summary 



The formation of prices in financial markets has long been a concern for economists. 
Macroeconomic factors, market microstructure, and investor preferences constitute an 
inexhaustive list of constituents that play a part. It is a daunting task to formulate 
a parsimonious pricing model that incorporates even this short list of effects, and is 
capable of delivering useful and timely results. It is not surprising then that analysis 
of price formation tends to focus on a single or small number of relevant elements at 
any one time. 

The information that traders and investors have about an asset is reflected in its 
price. Information about an asset might include information about any of the effects 
that influence the formation of its price. The arrival of new information leads to 
changes in the price of the asset. Qualitatively speaking, if information about an asset 
arrives infrequently and in large lots, then its price process will exhibit large jumps. 
Conversely, if information arrives smoothly and steadily, then the impact of information 
arrival over short time-scales will be modest. Thus, the emergence of information as 
driver of price dynamics presents itself as interesting avenue of investigation. 

The objective of this thesis is to provide a framework for the derivation of price 
dynamics of assets (or, indeed, the valuation dynamics of liabilities) through the mod- 
elling of information flow. It has become commonplace in the mathematical finance 
literature to develop pricing models under the risk-neutral measure. Two situations 
when this may be appropriate are (a) when the market is incomplete and there exists 
a multitude of equivalent martingale measures, and the selection of any one for pric- 
ing is made subjectively; and (b) when pricing models are calibrated to market prices 
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1 Introduction and summary 



(of options), since these prices are (theoretically) expectations under the risk- neutral 
measure. We adopt a similar approach. Throughout this work we build models under 
a 'suitable' measure Q. At time t, we define the price of a contingent claim X? due at 
time T as P^ Eq [X^ | J- t ], where P t r is a discount factor, and {Ft} is the market fil- 
tration. For this price to agree with the finance-theoretic arbitrage-free price, then the 
measure Q needs to be interpreted as the T-forward measure in the case of stochastic 
interest rates, or the risk-neutral measure in the case of deterministic interest rates. 
If Xt is the size of a future liability that needs to be accounted for ('booked'), then 
it may not be appropriate to use the market (i.e. arbitrage-free) price for valuation. 
This is particularly true if the market for such a claim is illiquid. In order to remain as 
general as possible we refrain from being precise in the interpretation of Q, and leave 
that to the judgement of the implementer. 

In a stochastic model, it is the filtration {F t } that encodes the emergence of infor- 
mation. Choosing an asset price process to be geometric Brownian motion, for example, 
implies that the process is adapted to a Brownian filtration. Although this approach 
of implicitly choosing a filtration by specifying the law of a price process is common in 
mathematical finance, we wish to avoid it and to specify {Ft} directly. In particular, 
we postulate the existence of a market information process {£,tr} which generates {J~t}- 
Then prices are derived by taking conditional expectations of cash flows with respect 
to this filtration. 

Earlier, we were careful to include the valuation of liabilities within our remit. We 
will consider in detail how the methods we develop can be applied to the calculation of 
reserves that an insurance company should set aside to meet future claims. A non-life 
insurance company may underwrite various risks for a particular year in return for 
premiums. The company 'incurs' claims over the one-year period, which means that 
losses are triggered that the company is liable to cover. However, there may be a delay 
between the date a loss is triggered and the day it is reported to the company, the 
total size of a claim may not be known when the claim is reported, and a claim may 
not be paid by a single cash-flow on a single date. The insurer may be paying for these 
losses for many years to come. The problem is: How much money should the insurer 
reserve at a given time to cover all future claim payments? This has implications for 
the company's accounting, tax liability, solvency, capital adequacy, and investment 
strategy. 
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Chapter [2] begins with a brief introduction to Levy processes, strictly stable pro- 
cesses, and Levy bridges. A Levy bridge is a stochastic process defined over a finite 
time horizon, and is a Levy process whose terminal value is known from the outset. 
We provide a proof of the Markov property for Levy bridges. For the remainder of the 
chapter we examine particular examples of these processes. First are the well-known 
Brownian motion and the Brownian bridge; these are Gaussian processes with con- 
tinuous sample paths. Brownian motion without drift is a stable process. Then are 
the gamma process and gamma bridge, which are increasing processes. The variance- 
gamma (VG) process is closely related to Brownian motion and the gamma process. 
In particular, a VG process is constructed by subordinating a Brownian motion by an 
independent gamma process, or by taking the difference of two independent gamma 
processes. We list some properties of the VG process and then examine VG bridges. 
We are able to utilise the relationship between the VG process and Brownian motion 
and the gamma process to derive two constructions of the VG bridge. These construc- 
tions write the VG bridge in terms of Brownian bridges, gamma bridges, and random 
volatility terms. The second stable process we examine is the stable- 1/2 subordinator. 
This process can be constructed as the hitting times of a Brownian motion. We choose 
to examine this subordinator over other stable subordinators because its density is 
known in closed form. It will become apparent that this is a desirable feature in this 
work. The stable-1/2 subordinator has heavy tails; the expected value of the process 
at any future date is infinite. However, all positive moments of the stable-1/2 bridge 
exist. By subordinating a Brownian motion with an independent stable-1/2 subordina- 
tor we can construct a Cauchy process. The Cauchy process is the third stable process 
we consider. It too is heavy tailed — its first moment does not exist. We find that the 
first two moments of the Cauchy bridge do exist and are finite. One of the effects of 
pinning the end point of the Cauchy process is to temper its wild behaviour. The last 
process we consider that has a continuous state space is the normal inverse-Gaussian 
(NIG) process. This process is constructed by subordinating a Brownian motion by an 
inverse-Gaussian (IG) process. It is a very similar process to the VG process. It is not 
surprising then that NIG bridges are similar to VG bridges. Finally, we come to the 
Poisson process and the Poisson bridge. These processes have a discrete state spaces. 
The state space of the Poisson process is N , and Possion bridges are restricted to a 
subset of Nq. We show how the Poisson bridge can be written as a Poisson process 
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under a time change that is not independent of the process. 

In Chapter [3] we define Levy random bridges (LRBs). An LRB is a process defined 
over a finite time horizon whose bridge laws are the bridge laws of a Levy process. It 
can be interpreted as being a Levy process conditioned to have a fixed marginal law 
at a fixed future date. The motivation for LRBs is the desire for Markov processes, 
defined over a fixed time interval, over which we have a distributional degree of freedom 
for the terminal value. We shall see that such processes are useful for asset pricing in 
the information-based framework. We derive various properties of LRBs including: 
that LRBs are Markov processes; that the law of an LRB is equivalent to the law of 
a Levy process (at least for all times up to the LRB's termination time); that LRBs 
have stationary increments; that the joint distribution of the increments of an LRB 
have a generalized multivariate Liouville distribution; and, if the path of an LRB is 
split into non-overlapping portions, then each portion is itself an LRB. The marginal 
characteristic function often proves convenient in the analysis of a Levy process. This 
is not the case for an LRB. However, we are able to provide an expression for the 
transition law of a general LRB. 

The information-based framework of Brody, Hughston & Macrina is described in 
Chapter HI In this framework, cash flows are functions of independent A-factors. 
'Information processes' generate the market filtration, and reveal the value of the X- 
factors — thus, they reveal the value of the cash flows. Cash flows are priced by discount- 
ing their expected value given the market information. A general multi-factor set-up is 
described that allows for a rich dependency structure between cash flows. Much of the 
analysis is presented for a single A-factor market where the only information process 
is an LRB. The final value of the LRB is set to be the value of the A-factor. Using a 
single factor may sound restrictive, but this includes some well-known one- dimensional 
exponential Levy models as special cases (the Black-Scholes model was recovered from 
an information-based model by Brody et al. [19]). Using Bayesian methods, we are able 
to derive the dynamics of the price of a cash flow. From this, we can price European 
call options on the cash flow price. An expression for the call price is given in a general 
LRB- information model. The material in Chapters |3] and H] appears in 52]. 

In Chapter [5] we examine a one-parameter VG random bridge and apply it to 
information-based pricing. We derive two terminal value decompositions of the VG 
random bridge using the decomposition of the VG bridge. We show that a three 
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parameter VG process can be recovered by scaling a standard VG random bridge which 
las an asymmetric VG terminal law. This allows the VG equity model of Madan et al. 



671 ] to be derived as a special case of a single X-factor information-based model, when 



information is provided by a standard VG random bridge. We price a binary bond in 
a VG information model and include a rate parameter a which acts like a volatility 
parameter. We provide two algorithms for the simulation of sample paths of the VG 
random bridge, and provide plots of simulated binary bond prices. 

We develop a non-life reserving model in Chapter [H] using a stable- 1/2 random 
bridge to simulate the accumulation of paid claims, allowing for an arbitrary choice of 
a priori distribution for the ultimate loss. Taking a Bayesian approach to the reserving 
problem, we derive the process of the conditional distribution of the ultimate loss. 
The 'best-estimate ultimate loss process' is given by the conditional expectation of the 
ultimate loss. We derive explicit expressions for the best-estimate ultimate loss process, 
and for expected recoveries arising from aggregate excess-of-loss reinsurance treaties. 
Use of a deterministic time change allows for the matching of any initial (increasing) 
development pattern for the paid claims. We show that these methods are well-suited 
to the modelling of claims where there is a non-trivial probability of catastrophic loss. 
The generalized inverse-Gaussian (GIG) distribution is shown to be a natural choice 
for the a priori ultimate loss distribution. For particular GIG parameter choices, the 
best-estimate ultimate loss process can be written as a rational function of the paid- 
claims process. We extend the model to include a second paid-claims process, and 
allow the two processes to be dependent. The results obtained can be applied to the 
modelling of multiple lines of business or multiple origin years. The multidimensional 
model has the attractive property that the dimensionality of calculations remains low, 
regardless of the number of paid-claims processes. An algorithm is provided for the 
simulation of the paid-claims processes. The material in this chapter appears in 53]. 

In Chapter [7] we derive properties of Cauchy random bridges, and describe a method 
for simulating sample paths. Modelling the information process as a Cauchy random 
bridge, we examine the pricing of a binary bond in an information-based model. We 
derive an explicit expression for the price of a call option on the bond price. 

Chapter [8] closely follows Chapter [5j Only, in this case, we examine the NIG random 
bridge, and use it to model an information process, as opposed to using the VG random 
bridge. Since the NIG and VG processes are similar, the information-based models their 
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respective random bridges produce are similar. 

In Chapter [9] we provide an example of an LRB with a discrete state-space, the 
Poisson random bridge. Poisson random bridges are counting processes, and we derive 
expressions for the waiting times between jumps in terms of the probability generating 
function of the terminal distribution. We show that a Poisson random bridge can be 
written as a Poisson process with a state-dependent intensity, and we derive an explicit 
expression for the intensity process. When the terminal distribution of a Poisson 
random bridge is a negative binomial distribution, we show that all of the increment 
distributions of the process are negative binomial. We then generalise this result to 
show that a Poisson random bridge with a mixed Poisson terminal distribution is a 
mixed Poisson process. That is, the distribution of any increment of the process is 
a Poisson distribution with a mixed mean. By making the jump sizes of the PRB 
random we construct the compound Possion random bridge. We derive an expression 
for the characteristic function of compound Poisson random bridge. Finally, we price 
an nth-to-default credit swap in a model where defaults occur at the jump times of a 
Poisson random bridge. In this credit swap the buyer pays a premium in return for a 
lump-sum payment on the event of the nth default from a basket of credit risks. 



Chapter 2 

Levy processes and Levy bridges 



We fix a probability space (0, Q, J 7 ), and assume that all processes and nitrations under 
consideration are cadlag. Unless otherwise stated, when discussing a stochastic process 
we assume that the process takes values in R, begins at time 0, and the filtration is 
that generated by the process itself. We work with a finite time horizon [0,T]. 

2.1 Levy processes 

This section and the next summarise a few well-known results about one-dimensional 
^evy processes and stable processes, further details of which can be found in Bertoin 
13j|, Kyprianou [61], and Sato A Levy process is a stochastically-continuous pro- 
cess that starts from the value 0, and has stationary, independent increments. An 
increasing Levy process is called a subordinator. For {L t } a Levy process, its charac- 
teristic exponent ^ : R — > C is defined by 

E[e iAit ] = exp(-ttf(A)), A G R. (2.1) 

The characteristic exponent of a Levy process characterises its law, and its form is 
prescribed by the Levy-Khintchine formula: 

1 f°° 

\ff(A) = iaA + -a 2 A 2 + / (1 - e lxX + ix\t { \ x \ <1} )U(dx), (2.2) 
^ J -00 

where a G R, a > 0, and IT is a measure (the Levy measure) on R\{0} such that 

(1 A |x| 2 ) il(dx) < 00. (2.3) 
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2.1 Levy processes 



There are particular subclasses of Levy processes that we shall consider, defined as 
follows: 

Definition 2.1.1. Let {L t }o<t<T and {M t }o<t<T be Levy processes. Then we write 

1. {L t } G C[0,T] if the density of L t exists for every t G (0, T], 

2. {M t } G T> if the marginal law of M t is discrete for some t > 0. 

Remark 2.1.2. If the marginal law of M t is discrete for some t > 0, then the marginal 
law of M t is discrete for all t > 0. The density of L t exists if and only if its law is 
absolutely continuous with respect to the Lebesgue measure. In general, the absolute 
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continuity of L t depends on t; thus C[0,T 2 ] C C[0, Ti] for T\ < T 2 . See Sato 
chap. 5] for further details on the time dependence of distributional properties of Levy 
processes. 

We reserve the notation ft{x) to represent the density of L t for some {L t } G C[0, T\. 
Hence f t : R — > K + and Q[L t G dx] = f t (x)dx. We reserve Qt(a) to represent the 
probability mass function of M t for some {M t } G T>. We denote the state-space of 
{M t } by {cii} C R. Hence Q t : {aj [0, 1] and Q[M t = ai} = Q t {(Xi). We assume that 
the sequence {a{} is strictly increasing. 

The transition probabilities of Levy processes satisfy the convolution identities 

/oo 
ft-,(x - v)Uv) dy for {L t } G C[0,T], (2.4) 

-oo 

and 

oo 

Qt(a n ) = ^2 Qts( a n ~ a m)Qs(a m ) for {M t } G V, (2.5) 

m=— oo 

for < s < t < T. These are the Chapman-Kolmogorov equations for the processes 
{L t } and {M t }. 

The law of any cadlag stochastic process is characterised by its finite-dimensional 
distributions. The finite-dimensional densities of {L t }o<t<T exist and, with the under- 
standing that xq = t = 0, they are given by 

n 

h G dxi, . . . , L tn G dx n ] = J"[ [fu-ti-i ( x i - X i-l) & X i] > ( 2 -6) 

i=l 



2.2 Stable processes 
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for every n G N+, every < t± < • • • < t n < T, and every (x\, . . . , x n ) G W 1 . With the 
understanding that ak = to — 0, the finite-dimensional probabilities of {M t } are 

n 

Q[M tl =a kl ,...,M tn = a kn ] = JI^-iK ~ ^-J, (2-7) 

i=i 

for every n G N+, every < t± < ■ • • < t n , and every (ki, . . . , k n ) G Z n . 



2.2 Stable processes 

The (strictly) stable processes form a subclass of the Levy processes. We say that a 
Levy process {S°} is a stable process with index a (or stable-a process) if its charac- 
teristic exponent satisfies 

tf(JfeA) = M(A), (2.8) 

for every > and every A 6 1; a is restricted to values in (0,2]. The process {5"f} 
satisfies the scaling property 

{k-^S^o = {S?} t > for k > 0. (2.9) 

Equation (12. 8 j) and the Levy-Khintchine formula restrict the characteristic exponent 
and the Levy measure of to take an explicit form which depends on a. When 

ae (0,1)U(1,2), 

*(A) = K\\\ a (l - i/3sign(A) tan(7ra/2)), (2.10) 
where k > 0, and /3 G [—1, 1]. In this case the Levy measure can be written 

dx for x > 0, 

(2.11) 

/t |x| a dx for x < 0, 
where k + and k~ are non-negative numbers satisfying 

If /3 = 1 then the process exhibits only positive jumps, if /3 = the process is symmetric, 
and if j3 = — 1 the process exhibits only negative jumps. If a = 2 then \I/(A) = |cr 2 A 2 , 
and n(dx) = 0. In this case {S'f } is a Brownian motion (without drift). If a = 1 then 
^(A) = iaA + c|A|, for c > 0, and II(dx) = cx _2 dx. In this case {S'f} is a Cauchy 
process with drift. 
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2.3 Levy bridges 



Excluding the case when a = 2, stable processes are heavy-tailed processes. The 
fractional moments of the stable random variable (a < 2) satisfy 

E[\S?\ p ]<oc iip<a, (2.13) 

E[\S?\ p } = oo iip>a. (2.14) 

Hence the second moment of is infinite for a < 2, and the expected value of 
does not exist or is infinite for a < 1. 

The density of exists and is continuous for any a, and so {S^} G C[0,T] for 
any T > 0. However, this density can be expressed in terms of elementary functions 
only when {5*"} is a Brownian motion, a Cauchy process, or a stable subordinator 
with index a = 1/2. We will examine each of these special cases later in this chapter. 
(See Feller 0, XVII. 6] for examples of series representations for the density of a stable 
random variable with arbitrary index a.) 

If {(Sf } is a stable subordinator then the Laplace transform of exists and is given 

by 

E[exp (-ASJ*)] = exp (-Kt\ a ) for A > 0, (2.15) 
where k > 0, and a must be further restricted to < a < 1. 



2.3 Levy bridges 

A bridge is a stochastic process that is pinned to some fixed point at a fixed future time 



Bridges of Markov processes were constructed and analysed by Fitzsimmons et al. 41 1 
in a general setting. In this section we focus on the bridges of Levy processes in the 
classes C[0, T] and D. In particular we have the following: 

Proposition 2.3.1. The bridges of processes in C[Q,T] andT> are Markov processes. 

Proof. We need to the show that the process {L t } 6 C[0,T] is a Markov process when 
we know that Lt = x, for some constant x such that < fr{x) < oo. (It will be 
explained later why the condition that < fr(x) < oo is required to ensure that the 
law of the bridge process is well defined.) In other words, we need to show that 

Q [L t < y | L tl = xx, . . . , L tm = x m , L T = x] = <Q> [L t < y | L tm = x m , L T = x] , (2.16) 



2.3 Levy bridges 
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for all m G N+, all (xi, . . . , x m , y) G R m+1 , and all < ti < • • • < t m < t < T. It is 
crucial to the proof that {L t } has independent increments. Let us write 



\ — L u — L ti _ 1 , 

&i Xi Xi—i, 



(2.17) 
(2.18) 



for 1 < i < m, where to = and xo = 0. Then we have: 

Q [L t < y | L tl = xi, . . . , L tm = x m , L T = x] 

— Q[L t - L tm < y - x m | Ai = Si, . . . , A m = S m , L T - L tm = x - x m ] 
= Q [L t - L tm < y - x m | L T - L tm = x - x m \ 

J t ~ L tm < y — x rn | Lt — L tm = x — x m , L tm = x rn ] 
H<y\L T = x, L tm = x m ] . (2.19) 

The proof for processes in class V is similar. □ 

Let {Lt} G C[0,T], and let {-^|r}o<t<T be an {L t }-bridge to the value z G K at 
time T. For the transition probabilities of the bridge process to be well defined, we 
require that < fr(z) < oo. By the Bayes theorem we have 



1$ G dy 



L, sT — X 



Q [L t Edy\L s = x,L T = z] 

= Q [L t G dy,L T G dz\L s = x] 
Q [L T Edz\L s = x] 
ft-s(y-x)f T -t(z-y) 

= 7 — t v d v, 

f T -s{z - x) 

for < s < t < T. We define the marginal bridge density ftriy, z) by 

ft{y)fr-t(z - y) 



(2.20) 



ftrfaz) = 



In this way 



Lg G dy 



t(z) _ 
Lj sT — X 



h{z) 



ft- 9 ,T-,(y - x; z - x) dy. 



(2.21) 



(2.22) 



The condition < fr(z) < oo is enough to ensure that 



V ^ ft-s,T- s (y -L ( ^;z- L { ^> 



(2.23) 
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2.3 Levy bridges 



is a well-defined density for almost every value of L^. To see this, note that 



oo />oo 



t-s,T- 



,(y -x;z-x) 



— oo J — oo 



G dx 



dy 



fr(z 

From f T2T24j) it follows that 



/oo /*oo 
/ fts,T~s(y -x;z- x)f SjT {x; z) dx dy 
oo J —oo 

fr-t(z - y) 

f T ~t(z -y)f t (y) dy= I. 



oo 

1 



ft- s {y - x)f s {x) dx dy 



/oo 
ft-a,T-*(y -L ( ^;z- 1$) dy = 1 
-oo 



1. 



(2.24) 



(2.25) 



M$ = a,- 



M$ = a, 



Let {M 4 } G D, and let {M^} < 4 < T be an {M t }-bridge to the value at time T, 
so Q[Myy = 0^] = 1. For the transition probabilities of the bridge to be well defined, 
we require that Q[Mr = Ofc] = Qr^k) > 0. Then the Bayes theorem gives 

Q[Af t = o i |Af a = a i ,M r = oj fc ] 

_ Q [M f = Qj, M T = dk\ M s = dj] 
Q[M T = a fc | M s = 
Qt-sjo-j ~ «i)Qr-tK ~ Qj) ,„ 9fi ^ 

where s, t satisfy < s < t < T. Note that if Qria-k) — 0, then the ratio (12.261) is not 
well defined when s = 0. 

We provide sufficient conditions for the integrability of Levy bridges: 

Proposition 2.3.2. If there exists a constant C < oo such that \x\ l+ ^ft(x) is bounded 
for \x\ > C and all t G (0, T\, then 



\x 



l+2a 



ffr{x\ z) dx < oo, 



for every a G (0, 0). In other words, 

, s l+2a 



E 



L 



< oo (0 < t < T). 



(2.27) 



(2.28) 



Similarly, if there exists a constant C such that \ai\ 1+/3 Q t (ai) is bounded for \i\ > C 
and all t G (0, T], then 

l+2a 



E 



M 



W 

tT 



< oo (0 < t < T). 



^2.29) 
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Proof. We prove the proposition in the continuous case. The discrete case is similar. 
The cases t — and t = T are trivial, so we will assume that t G (0, T). 
Fix a G (0, 0) and assume that 

|x| 1+/3 / t (:c) < K < oo for \x\ > C and all t G (0, T]. (2.30) 

First we prove that / \x\ a f t (x) dx < oo for t G (0, T). We have 

kr/t(^)da;= / \x\ a f t (x) dx + / |a;| a / t (:r) drr 

-oo J-C J(-oo,-C*lU[C*,oo) 



<2C° + 2Kl 



oo ^.a 



/3 -a 

For ?/ 6 R, we can generalise this to 



2K 

2C a + C a ~P < oo. (2.31) 



/OO f*OQ 
\x + y\ a f t (x)dx< / (\x\ + \y\) a f t (x)dx 
-oo J —oo 

= f (M + \y\Tft(x) dx + f (\x\ + \y\Yf t {x) dx 

J-y J (~oo,-y]U[y,oo) 

/oo 
\x\ a f t (x) dx < oo. (2.32) 
-oo 



Finally, we have 

r-oo r C 



/oo ply n 

\x\ 1+2a f tT (x; z)dx= / \x\ 1+2a f tT (x; z) dx + / \x\ 1+2a f tT (x; z) dx 

-oo J-C ^(-oo,-C]U[C,oo) 

< 2C 1+2a + — — / \x\ a f T -t{z - x) dx 
K f°° 

= 2C 1+2a + — - / \y - z\ a f T - t (y) dy < oo. (2.33) 

□ 



2.4 Stable bridges 

Bridges of stable processes inherit a scaling property: 

Proposition 2.4.1. Let {S t } be a stable process with index a, and let k > be a 
constant. If {S^} is a bridge of {S t } to the value z at time T, and {S^ T } is a bridge 
of {S t } to the value ( = k x l a z at time kT, then 

{S;?} = {&. (2.34) 
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2.5 Brownian motion and Brownian bridge 



Proof. Denote the density of St by ft(x). From the scaling property of stable processes, 
we have 



It follows that 



From (I2.20p we have 



G dy 



S S T ~ X 



St ^ k~ 1/a S kt . 



f t (x) = k l ' a f kt (k 1/a x) 



ft-s{y - x)f T -t(z - y) 



(2.35) 
(2.36) 



dy 



f T -s(z ~ X) 

,/ Jkt-ks (k 1/a y - k^x) f kT _ kt (C - k 1/a y) 
fkT-ks (C - k l / a x) 

K D kt,kT fc U i/ ft D fcs,fcT — x 



(2.37) 



□ 



2.5 Brownian motion and Brownian bridge 
2.5.1 Brownian motion 

Brownian motion is a Levy process, and is a Gaussian process (i.e. all of its finite- 
dimensional distributions are multivariate normal). Gaussian processes are charac- 
terised by their mean and covariance functions. In the case of a one-dimensional 
Brownian motion {B t }, these are 

E[B t ] = 6t, Cov[B s , B t ] = a 2 min(s, t), (2.38) 

where 9 G R is the drift of {B t }, and o > is the diffusion coefficient. The characteristic 
exponent of {B t } is 

(A) = -i8X + ^a 2 X 2 . (2.39) 

The density of B t is 

The sample paths of Brownian motion are continuous but nowhere differentiable. When 
8 = 0, {B t } is a stable process with index a = 2, and satisfies the scaling identity 

{k- l ' 2 B kt } ^ {B t } for k > 0. (2.41) 



2.6 Gamma process and gamma bridge 
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When 8 = and o = 1 we say that {B t } is a standard Brownian motion (or Weiner 
process) . 

2.5.2 Brownian bridge 

A Brownian bridge is also a Gaussian process. Let {Ptr}o<t<T be a standard Brownian 
bridge to the point z G R. The mean and covariance functions of {13$} are 

E[/3«] = Iz, Cov[£M $>] = min(S) t ) - |. (2.42) 
It follows that 

{k- 1/2 &T Z) } = {/$} forA:>0. (2.43) 
Let {W t } be a standard Brownian motion, and define the process {W^} by 

W$ = W t + ^(z-W T ) (0<t<T). (2.44) 

Calculating the mean and covariance functions for this process verifies that it is a 
standard Brownian bridge to the value z at time T. It is also notable and easily 
verified that {W$} is independent of Wt- 



2.6 Gamma process and gamma bridge 
2.6.1 Gamma process 

A gamma process is a subordinator with gamma-distributed increments. The law of a 

gamma process is uniquely determined by its mean and variance at time 1. Both of 

these quantities are positive. Let {'jt} be a gamma process with mean 1 and variance 
m _1 > at time 1; then 

E[ 7t ] = t, Var[ 7t ] = t/m. (2.45) 

The density of j t is 

^(*) = 1 { , >0} ^^-V™ , (2.46) 
where T[z] is the gamma function, defined as usual for x > by 

POO 

T[x] = / u'Vdu. (2.47) 
Jo 
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2.6 Gamma process and gamma bridge 



Hence {jt} G C[0, T] for any T > 0. Due to the scaling property of the gamma 
distribution, if k > then the process {^7*} is a gamma process with mean k, and 
variance m~ l K 2 at £ = 1. The characteristic exponent of {7 t } is 

#(A) = m log[l - iX/m], (2.48) 

and so the characteristic function of 7 t is 

E[e iA7t ] = (1 - i\/m)~ mt . (2.49) 

In the limit m — > oo this characteristic function is e lA *, which is the characteristic 
function of the Dirac measure centred at t. It follows that {74} {£} asm-) 00. 

It should be noted that we find it convenient here to use a somewhat different 
parametrisation scheme for the family of gamma processes from that presented in 



Brody et al. 20]. In their scheme the 'basic' gamma process {74} is characterised by a 
single parameter m, with units of inverse time, such that Epyt] = nit and Varpyt] = mt. 
The 'general' gamma process is then obtained by considering a 'scaled' process {ftjt} 
where k > is a constant. Clearly E[«7t] = Kmt and Var[/«7 t ] = n 2 mt. Thus, if 
the mean rate fi and variance rate a 2 of a gamma process is specified, then we have 
m = /i 2 /er 2 and k = a 2 /fi. The scheme introduced in (12.451) above is equivalent to the 
choice k = vrC x in the BHM scheme. The advantage of the choice (12.451) for present 
purposes as the basic process is that it gives {7^ } the dimensionality of time, and hence 
makes it suitable as a basis for a time change. 

2.6.2 Gamma bridge 

Gamma bridges exhibit a number of remarkable similarities to Brownian bridges, some 
of which have been presented by Emery & Yor [33|] . Let {^tT}o<t<T be a gamma bridge 
with final value 1 associated with the gamma process {7*}. The transition law of {7*t} 
is given by 

Q [jrr e dy | 7 sT = x] = Q [j t e dy \ 7 S = x, 7 T = 1 ] 

= 9t-l(y - x)gj^ t (i -y) 

= t{ x<y< i} , l ~ X r , 77 . -rdy, (2.50) 

(1 — x)Ys[m{t — s), m[T — t)\ 



2.6 Gamma process and gamma bridge 
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for < s < t < T and x > 0. Here B[a, f3] is the beta function, defined for a, f3 > by 

B[a,fl= / 1 x- 1 (l-x)^dx = ^M. (2.51) 

If the gamma bridge {jtr} has reached the value x at time s, then it must yet travel a 
distance 1 — x over the time period (s,T\. Equation (I2.50P shows that the proportion 
of this distance that the gamma bridge will cover over (s, t] is a random variable with 
a beta distribution (with parameters a = m(t — s) and = m(T — t)). The conditional 
characteristic function of ■jtr is 

E [e iA7tT | 7sT = x]= M[m(t - s),m(T - s),i(l - x)A], (2.52) 

where M[a, (3, z] is Kummer's confluent hypergeometric function of the first kind, which 
can be expanded as the power series 0, 13.1.2] 

„ , r n n a a;(a + 1) z a(a + l)(a + 2) z 

M la , M = l + r+ ^^ + A_>^ _ + .... (2 .53) 

For brevity, we will later refer to ( I2.53P as 'Kummer's function M[a, (3, z}\ Taking the 
limit as m — > oo in (12.521) . we have 

fc=o ^ ' 

= exp^i^^(l-x)A^), (2.54) 

which is the characteristic function of the Dirac measure centred at (l—x)(t—s) / (T—s). 
It then follows from the Markov property of gamma bridges that { ltT } ^ {t/T} 
as m — > oo. It is a property of gamma processes that the renormalised process 
{lt/lT}o<t<T is independent of jt (indeed, this independence property characterises 
the gamma process among Levy processes). This leads to the remarkable identity 

-) = Ht}. (2.55) 
It J 

The identity (12.551) can be proved by showing that the process on the left-hand side is 
Markov, and then verifying that its transition law is the same as (I2.50p . which can be 
done using the results in Brody et al. [20|. Two further properties of gamma bridges 
follow immediately from (I2.55p . The first is that the bridge of the scaled gamma process 
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2.7 Variance-gamma process and variance-gamma bridge 



{fir/t} is, for any k > 0, identical in law to the bridge of the unsealed process. Hence, 
for fixed m and fixed terminal value, the bridge of the gamma process defined by f!2.45p 
is identical to the bridge of the basic BHM gamma process. The second property is 
that a {7t}-bridge to the value z > at time T is identical in law to the process {z^tr}- 



2.7 Variance- gamma process and variance- gamma 
bridge 

2.7.1 Variance-gamma process 

A variance-gamma (VG) process is a Brownian motion with drift, subordinated by an 
independent gamma process. Letting {W t } denote a standard Brownian motion, we 
define the process {Vt} by 



V t = aW~ /t + 07 t a > and 9 G 



(2.56) 



Then {V t } is a VG process in its most general form (on the real line). The mean of 
the gamma process {7t} at t = 1 was fixed as unity, but (in terms of the law of the 
VG process) varying this is equivalent to an appropriate change of the parameters a 
and 9. When 9 = we say that {V t } is a symmetric VG process; if, in addition, o = 1 
then we say that {V t } is a standard VG process. The characteristic exponent of {V^} 
is given by 

\P(A) = mlog 



m a 2 \ 2 

1 — + 

m 2m 



where m~ l is the variance of 71. This can be decomposed as 
\I/(A) = m log 

where 





+ m log 


[l + ^1 


m 




m 



(2.57) 



(2.58) 



fi + = X - (V9 2 + 2m<r 2 + 9 s ) 



(2.59) 



and 



2 + 2ma 2 - 9 



(2.60) 
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The right-hand side of ( I2.58P is the sum of two characteristic exponents. The first 
corresponds to a gamma process, and the second corresponds to a decreasing Levy 
process whose absolute value is a gamma process. It follows that, for {7$ } and 
{7$ } independent copies of {7*}, we have 

„,(+) „,(-) 



{Vt} = (rf'-Hr)- (2-61) 
The characteristic function of Vt is 

In the limit as m — > 00, this characteristic function tends to 

exp (ieXt - \(T 2 X 2 t) . (2.63) 

It follows that 

{V t } ^ {aW t + et] as m -» 00. (2.64) 



The distribution of V is norma 



with a gamma-mixed mean and variance. The 



density of V t can be shown to be 67] 



— , mi 1 

2 fyi rnt ^ x / cr f x 



- 2 4 



(2.65) 

where is the modified Bessel function of the third kind. We see that {Vt} G C[0, T] 
for all T > 0. Each of the following facts about K v [z] can be found in Abramowitz & 
Stegun (J 9.6]: 

1. < K v (z) < 00 for 1/ G E and 2 > 0, (2.66) 

2. = (2.67) 

3. lim^!i = l, (2.68) 
z->o - log z 

K \z\ 

4. lim "[ J =1 for v > 0, (2.69) 
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2.7 Variance-gamma process and variance-gamma bridge 



n the case when v is half an odd number we have following expression for K v [z] from 
1 10.2.15]: 

^n+|W = \/|; e "*i> + for n G N, (2.71) 



where (to, n) is Hankel's symbol, 



to, n 



T[m + I + n) 
n\ T[m + \ — n] 



(2.72) 



Writing 



(m,0,l). 



mi 1 

2 m m * / x 2 \ 2 4 
7r T [mt] \ 2m J 



mt- 



V2 



mar 



(2.73) 



we have 



Am,8,<r) ( \ _ 1 6>z/<t 2 fr. \ 1— 2mt ,(m) / , ,\ 

it [X) — —e {K(m,0,a)) Jt [K(m,6,a)% / & ) 



(2.74) 



where 



k 



2ma 2 



(2.75) 



Here ft m \x) is the density of the standard VG random variable W^t). Since K u [z] is 
finite for z > 0, ft m ' e ' a \x) is finite away from zero. In other words, < ft m ' 9 ' a \x) < oo 
for all t > and all x G M\{0}. When x = 0, we have 



1 Am" / # 2 \ - 

+ oo 



r[mt - 1/2] 
T[mt] 



for t > (2m 



,-i 



^2.76) 



for < t < {2m)- 1 . 



2.7.2 Variance gamma bridge 

Let {V t } be a VG process with parameter set {m, 9, a}; and let {V^ } be a {Vj}-bridge 
to the value z G 1R at time T. Since we must have that the density of Vt is positive 
and finite at z, it follows from (I2.76P that either z ^ 0, or T > (2m) -1 . The transition 
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law of {vjj)} is given by 



V t P e dy 



sT 



k fl-l 



f^\z-x) 



dy 



x) 



))f$$tG(z-y)) 



a 



dy 



a 



k U^edy 



a 



rr(kz/a) 

k UST 



X 



(2.77) 



for < s < t < T, k = k( m fi,a), and {Uf£ x } is a {iy(7 t )}-bridge to kz/a at time 
T. This shows that a bridge corresponding to any VG process is identical in law to a 
scaled bridge of a standard VG process. We focus now on bridges of the standard VG 
process. It follows from ( I2.44p that 

law 



{W t } 



0<t<l 



{A + Wl}o<Kl, 



(2.78) 



where {f3 t } is a Brownian bridge, independent of W%, ending at the value at time 
t — 1. Also, from the time-scaling property of Brownian motion, 



{W(Xt)} l = {VXW(t)} 



(2.79) 



for X a positive random variable independent of {W^t}. We can use these identities, 
with some properties of gamma bridges, to derive a terminal-value decomposition of a 
standard VG process: For t G [0,T], we have 



{W{lt)} 



W 



"fr- 



it 



law 



law 



law 



law 



7t 

{y/jrW^tT)} 



{itTy/jT W(l) + v / 7r/3(7tT)} 
{j tT W(j T ) + a/7t 0(jtT)} ■ 



(2.80) 
(2.81) 
(2.82) 
(2.83) 



In the above, (I2.80P holds since {jt/lr}, and 7r are independent; (I2.8ip follows 

from (12775]) : ( 12782]) follows from fB77B|> : ( 12783]) also follows from (12773]) . Note that in 
( I2.83P the Brownian bridge {/3(i)}o<t<i, the gamma bridge {7t:r}o<t<T, and the random 
vector (j T , W(jt)) are independent. The joint density of (72-, W^t)) is 



exp 



If! 
2 j/ 



(2.84) 
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2.7 Variance-gamma process and variance-gamma bridge 



gp\y). (2.85) 



Hence, given W(jt) = z, the conditional density of 7t is 

1 / l2 2 ' 

V l— —r^ exp 

ft\z)V2^y V 2y 

Some simplification shows that this is the generalized inverse- Gaussian (GIG) density 
faiciy] mT - 1/2, \z\,V%m), where 

I'gig{ x 'i A, 5, 7) = 1 {X>0} (D ^^^ ^exp (-K^ 1 + 7 V)) ■ (2.86) 

It follows from ( I2.83P that the standard VG bridge satisfies 

{l/W} ^ w |z 7tT + v^/3(7 tT )} , (2.87) 

where S 2 has a GIG distribution with parameter set {mT — 1/2, \z\, \j2m\. 
From f)2.6ip we have that 

= V(7t-7 t )}, (2.88) 

where {74} and {7<} are independent, identical gamma processes with parameter m, 
and // = (m/2)- 1/2 . We can use this to derive an alternative representation of the 
standard VG bridge. Define the process {Gt} by setting Gt = jt — It- Then we have 

{G t } = {it - it} 

law r _ _ 

= XlTltT - iTltTS 

= {G T ltT + It (jtr - Jtr)} , (2.89) 

where {7^} and {jtr} are identical gamma bridges, independent of each other, and 
independent of the random vector ( 7r , Gt)- The joint density of (7^, Gt) is given by 

(y,z)^g^\ y )g^\ y + z). (2.90) 

Given that Gt = z, the conditional density of 7^ is 

where x + denotes the positive part of x, so (— z) + = max(0, — z). Then we have 

{U$} = {zjtr +Y Z » (7tT - Itr)} , (2.92) 
where Y z > is a random variable with density ( I2.9ip . 



2.8 Stable-1/2 subordinator and stable-1/2 bridge 
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2.8 Stable-1/2 subordinator and stable-1/2 bridge 

2.8.1 Stable-1/2 subordinator 

Let {St} be a stable-1/2 subordinator . The characteristic exponent of {St} is 

M/(A) = -||A| 1 / 2 (l-isign(A)), (2.93) 

where c > is a parameter related to k in f[2TT0]) by c = k^2. Then {S t } satisfies 
the scaling property {k~ 2 Skt} = {St}, for k > 0. The random variable St has a 'Levy 
distribution' with density 

M ' ) = 1 ™v£g*^{-rf)- (2J4) 

We call c the 'activity parameter' of {S t }. The density (I2.94p is bounded for alH > 
and is strictly positive for all x > 0, hence {St} G C[0, T] for all T > 0. Integrating 
(I2.94p yields the distribution function 

/ f t (y)dy = 2<5>[-ctx~ 1 / 2 }, (2.95) 
Jo 

where $>[x] is the standard normal distribution function. The random variable St has 
infinite mean, indeed E[Sf] < oo if and only if p < 1/2. The density of 1/5^ is 

x ^ 1{x>0} v2Wm x ~ 1/2 exp ( - ^ cVx ) ' (2,96) 

Thus the increments of {S t } are distributed as reciprocals of gamma random variables. 
For {Wt} a standard Brownian motion, define the exceedence times {r t }t>o by 

r t = mi{s : W s > ct}. (2.97) 

From Feller 0, X.7], we then have 

{St} = {r t }. (2.98) 

2.8.2 Stable-1/2 bridge 

Fix z > and let {S^}o<t<r be a bridge of the process {S t } to the value z and time 
T. We call {5^- } a stable-1/2 bridge. The density function of the random variable 



32 



2.8 Stable-1/2 subordinator and stable-1/2 bridge 



My; z) 



ft{y)h-t{z-y) 



f T (z) 



1 



{0<y<z}- 



1 c 2 (Ty-tz) 2 

1 ct(T-t) e *P{—2 U-v) 



V27T T (y- yt/zf' 2 

This density is bounded, and has bounded support, so 



(2.99) 



E 



< oo for p > 0. 



(2.100) 



Remark 2.8.1. From Proposition ^.^. 4 stable-1/2 bridges satisfy the following scaling 
property: For k > a constant, and {S^} a stable-1/2 bridge, we have 



( q(z)\ law j ,-2 q(k 2 z)\ 

r tT h<t<T ~ r ^W <*<t 



(2.101) 



Integrating the density (12.991) yields the following (details of the proof can be found 
in Appendix lA.lj) : 

Proposition 2.8.2. For y e [0,2], the distribution function of the random variable 
Sfl) is given by 



F tT {y-z) = $ 



c(Ty - tz) 

Vv z ( z - y) 



\- — \ e 2c 2 KT-t)/z $ 
T 



c((2t-T)y-tz) 
^yz(z - y) 



. (2.102) 



Remark 2.8.3. When t = T/2, the second term in the distribution function \2.102\) 
vanishes. The distribution function is then analytically invertible, and we obtain the 
identity 



g(z) law 1 . j 

T/2 ' T 2 ^ ' y/c?T 2 /z + Z 2 
where Z is a standard normal random variable. 



(2.103) 



Corollary 2.8.4. 



{Sfl}} {Ifz} as c — >■ oo. 



(2.104) 



Proof. Fix z > 0. It is sufficient to show that 



lim F tT (y; z) = l{T y >tz} for Lebesgue-a.e. y e (0, z), 



(2.105) 
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since this is equivalent to 



lim < 

c—too 



StT T Z 



< 6 



Define a by 



a 



1 for all t e [0, T] and any e > 0. 
(2t - T)y - tz 



y/yz(z - 2/L 

and note that a > for y £ (0,2;). The inequality [l, 7.1.13] states 



f oo 

- 2 I -t 



e~ l dt < 



x 



+ sjx 2 + 4/vr 



(z > 0), 



(2.106) 
(2.107) 

(2.108) 



from which we deduce 



e ^t(T-t)/z^_ ac] < e 2A(T- ( )/ 2 /2 



-a 2 c 2 /2 



7i"ac+ v^c 2 + 2/tt 



2 (Ty-tz) 2 

2 ex P I - c fefefe 



7Tac+ y / a 2 c 2_ +2/V' 
Since the left-hand side of (12.1091) is positive, we see that 

lim e^-^i-ac] = 0. 

Then we have 



(2.109) 



(2.110) 



lim Ftx{y] z) = lim $ 



c(Ty - tz) 



X 1 / c— »oo 



s/yzjz^y) 

= ^{Ty-tz>0} ~ \l{ Ty =tz}, 

which completes the proof. 

The proof of the following proposition can be found in Appendix IA.21 

(z) 

Proposition 2.8.5. Define the incomplete first moment of S" tT by 

M tT (y;z) = f uf tT (u;z)du (0<y<z). 
Jo 

Then we have 



■ac\ 



(2.111) 



□ 



M tT (y; z) = ^z{ <I> 



c(Ty - tz) 
y/yzjz - y) 

?(*) 



_ e 2c 2 t(T-t)/z $ 



c{(2t-T)y-tz) 
^Jyz(z-y) 



(2.112) 



(2.113) 



and the second moment of S^ T is given by 

2 



E 



S, 



^\l-c(T-t)e^ x l-^[- C Tz-^]\. 



(2.114) 



34 



2.8 Stable-1/2 subordinator and stable-1/2 bridge 



Corollary 2.8.6. 



E[S«] 
Var[5W] 



*(r - 1) 

J>2 



e 22 



27TC 2 T 2 



C 2 T 2 



(2.115) 
(2.116) 



Remark 2.8.7. In general, we have 



E 



<7 



5 



(^) 

sT 



;r 



T-t t-s 

7^ X + — Z, 

T-s T-s 



(2.117) 



and 



E 



5 



(2) 



X 



t-s 
T-s 



[z-x) 2 < l-c(T-t)e~ 



z — x) 



2tt 



z - X 



T-s ' 
\Jz — X 



(2.11* 



/or < s < t < T. 




Figure 2.1: Simulations of the stable-1/2 bridge demonstrating the influence of the 
activity parameter c. Qualitatively speaking, increasing the value of c decreases the 
frequency of large jumps, and increases the frequency of small jumps. 



2.9 Cauchy process and Cauchy bridge 
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2.9 Cauchy process and Cauchy bridge 
2.9.1 Cauchy process 

The Cauchy process is a stable process with index a = 1. Let {Z t } be a driftless 
Cauchy process. The characteristic exponent of {Z t } is 

#(A) = c|A| for XeR, (2.119) 

where c > is a scale parameter. The scaling property of {Z t } can be written 



{k^Zkt} '= {Z t } for k > 0. (2.120) 



The density function of Z t is 



ct 

n(y 2 ~c 2 t 2 )' 



Mv) = —* i ,2,2V ( 2 - 121 ) 



which is symmetric about y = 0, and the distribution function is 

f y 11 \V^ 

/ / t (x) dx = - + - arctan — . (2.122) 

We have E[|Z t |] = oo, but E[|Z t | p ] < oo for < p < 1. Let {5^} be a stable subordinator 
with index a = 1/2 (as described in the previous section). From Feller 37], X.9], if 
is a standard Brownian motion then we have 

{Z t } 1 ^ {W(S t )}. (2.123) 

2.9.2 Cauchy bridge 

Fix z G R and let {Z^} be a bridge of the Cauchy process {Z t } terminating at the 
value z at time T. The density of the random variable Z%p is 

f ( \ ft(y)fT- t (z-y) 
My;?) = — 



fr{z) 

ct(T - t) z 2 + c 2 T 2 

7tT (y 2 + c 2 t 2 )({z - y) 2 + c 2 (T - t) 2 ) ' 



(2.124) 



Integrating (12.1241) yields the following (details of the proof can be found in Appendix 
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2.9 Cauchy process and Cauchy bridge 



Proposition 2.9.1. The distribution function of Z^ T is 

. , 1 (T - t)(c 2 T(T - 2t) + z 2 ) 
FtT{y; Z)= 2 + nT(c 2 (T-2t) 2 + zi) ^ 



t(c 2 T(T - 2t) - z 2 ) 
H — r-r^ — arctan 



+ 



ttT(c 2 (T - 2t) 2 + z 2 ) 

ct(T - t)z 
ttT(c 2 (T - 2t) 2 + z 2 ) 



log 



z — y 
[cjT^t)\ 

y 2 + cH 2 



(z-y) 2 + c 2 (T-ty 



(2.125) 



where < t < T , c> 0. 



Remark 2.9.2. It follows from Proposition 2.4-1 that, for fixed k > 0, the Cauchy 
bridge {Z^} exhibits the scaling property 



k -l Z (kz) 



law f y(z) 
kt < kT( 0<t<T ~ I J 0<t<T 



(2.126) 



The first moment of a Cauchy process is not well defined, and the second is infinite. 
However, the first two moments of the Cauchy bridge exist and are finite. Hence 
conditioning a Cauchy process on its final value can temper its behaviour. The proof 
of the following proposition can be found in Appendix IA.41 

(z) 

Proposition 2.9.3. The first two moments of Z\f exist, and are given by 



E 



E 



Z 



w 

tT 
2 



T 
t 



z 2 + c 2 T(T - 1)). 



(2.127) 
(2.128) 



Remark 2.9.4. In general we have 



E 



J tT 



E 



/(*) 

J tT 



Zj sT — X 



7(2) _ ™ 
Zj sT — X 



T-t t-s 



-X 



T-s T-s ' 

l^x 2 + t -^{z 2 + c 2 {T-s){T-t)\ 
1 — s 1 — s 



(2.129) 
(2.130) 



for < s < t < T. 

Remark 2.9.5. The third moment of Z\ T does not exist; however 



E 



7(2) 



PI 



< 00 for p G (0, 3). 



(2.131) 



2.10 Normal inverse-Gaussian process and normal inverse-Gaussian bridge 
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2.10 Normal inverse- Gaussian process and normal 
inverse-Gaussian bridge 

2.10.1 Inverse- Gaussian process 

The inverse- Gaussian (IG) process is a subordinates: with Levy exponent 



tf(A) = c( v / 7 2 -2iA- 7 ), (2.132) 

for c > and 7 > constants. Let {X t } be an inverse- Gaussian process. The density 
function of X t is 

cte lct 1 

q t (x) = l {x>0 } r— 3/2 exp {-\{cH 2 x~ + j 2 x)) 



ct 1 I I 2 f c ^ 



t{x>o} v^^ exp \^v~i t] '■ (2 - 133) 



For fcGl, the moment = E[X t fc ] can be written 



,(*) 



ct\ k K k _ 1/2 [yct] 



7/ Ky 2 [icb] 



= ][l^ ct (~) +2 ( 2 - 134 ) 

Using the series representation of K n+ u 2 [z] given in (I2.7ip . the first four integer mo- 
ments simplify to 

mP = —, (2.135) 

7 

= — (1 + 7C*), (2-136) 

/■yO 

m f > = 4(3 + 37ct + 7W), (2.137) 

m (4) = + + g 7 2 c 2 t 2 + 7 3 c 3 t 3N (2.138) 

7 7 

The variance of X t is ct/'-f 3 . Hence, similar to the gamma process, the law of an IG 
process is characterised by its mean and variance at t = 1. 

For {Wt} a Brownian motion, define the exceedence times {ft} by 

T t = inf{t > : c~ x W t + c"V > t}. (2.139) 
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2.10 Normal inverse-Gaussian process and normal inverse-Gaussian bridge 



Then we have 



{X,} ^ w {rj. (2.140) 



Proposition 2.10.1. A bridge of the process {X t } to a fixed point z > at time T is 
identical in law to a bridge process of the stable- 1/2 subordinator {S t }. 

Proof. The (one-dimensional) marginal density of the IG bridge is 

q t (y)q T -t(z - y) 



Qt{z) 

1 c 2 (Ty-tz) 2 



ftT(y,z) 

= W47^= { T ( K 2/ \ 3/2 ' ■ (2-141) 

V2vr I (y- y 2 /z) 

This is identical to the marginal density of the stable-1/2 bridge given in (12.991) . It 
follows from f!2.22j) that the transition law of the IG bridge is identical to the transition 
law of the stable-1/2 bridge. □ 



2.10.2 Normal inverse- Gaussian process 

We construct the normal inverse- Gaussian (NIG) process by subordinating a Brownian 
motion by an independent IG process. Again let {X t } be an IG process, but this time 
set c = 7. This ensures that E[X f ] = t, so {X t } has units of time. Then we define the 
NIG process {Y t } by setting 

Y t = aW(X t ) + 9X t , (2.142) 

for a > and ffeK constants. In terms of the law of the NIG process, allowing the 
mean rate of {X t } to differ from unity is equivalent to an appropriate change in the 
parameters a and 9. We say that {Y t } is a standard NIG process if 9 = and a = 1. 
The characteristic exponent of {Y t } is 



(A) = cVc 2 + a 2 X 2 - 2i9X - c 2 . (2.143) 

The density of Y t is 



fM"\x) = *- ( fW , J f a * + d \ K i o--W(0 2 + c 2 a 2 )(cW + x 2 ) . (2.144) 
an V c z o~ z t z + x z 



2.11 Poisson process and Poisson bridge 
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Note that tf c ' e,a) (x) is bounded and continuous everywhere. Hence {Y t } G C[0,T] for 
all T > 0. The density of the standard NIG variable W(X t ) is 



ft\x) = fr' i \x) 

c 2 te cH 



c(c,0,l)/ 



n^cH 2 + X 2 
After some rearrangement, we find 

k 



cV 'cH 2 + x 2 



a 



where k > 1 and a > are given by 

k 2 = c~We 2 + c 2 , 



a 



2 = cVe 2 + c 2 . 



(2.145) 



(2.146) 



(2.147) 



2.10.3 Normal inverse- Gaussian bridge 

Let {Y t } be an NIG process with parameter set {c, 6, cr}, and let {Y^ '} be a {F t }-bridge 
to the value zGMat time T. The transition law of {Y^)} is 

f^\y-x)f^r\z-y) 



YP G dy 



y( z ) - T 
1 sT — x 



/t- s '(z-x) 

kftl (*(*-*))/&(*(*-*)) 



(7 



ft!. Hi*-*)) 



dy 



\UP M e dy 



k 



sT 



X 



(2.148) 



where (i) < s < t < T, (ii) k and a are given by (12.147p . and (iii) {U^^^} is a 
standard NIG bridge to kz/a at time T with parameter a. Thus, the three parameter 
NIG bridge is identical in law to a scaled standard NIG bridge. 



2.11 Poisson process and Poisson bridge 
2.11.1 Poisson process 

Let {N t } be a Poisson process with intensity A > 0. Then {N t } is a Levy process and, 
for each n G No, Q[N t = n] = Qt{n) where 

e~ xt (Xt) n 

Qtin) = 1 { „>o} (2.149) 
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2.11 Poisson process and Poisson bridge 



The characteristic exponent of {N t } is 

*(u) = A(l -e iu ). 



(2.150) 



Since {N t } has state space No and is increasing in increments of 1, it is a so-called 
counting process. The moment generating function of N t exists, and is given by 



E [e uNt ] = exp {Xt{e u - 1)) for uGl. 



(2.151) 



We have 

E[JV t ] = Var[iV t ] = Xt. (2.152) 

All higher moments of N t are finite, and can be found by differentiating (I2.15ip . Define 
the sequence {Tj}°^ of jump times of {N t } by 



Ti = inf {t > : N t > i} 



(2.153) 



and the let {Wi} be the sequence of waiting times W{ = Ti — Tj_x- The Wj's are 
independent, identically distributed random variables with 

(2.154) 



Q[Wi >t]= Q[N t = 0] = e~ xt , 
i.e. the waiting times of Poisson processes are exponentially distributed. 



2.11.2 Poisson bridge 

Let {iV t ^ }o<t<T be a Poisson bridge to the value k G N + at time T (we are excluding 



the case where = 0). The distribution of N^' is 



(fc) 



'N {k) - 7 



j,N T = k] 



\it = k] 

QtU)Q T -t(k-j: 



i 



QAk) 
f k 



{0<j<k} 



T 



T 



(2.155) 



for < £ < T. Hence has a Binomial ( fc, t/T) distribution (note that this holds 
for any choice of A). The probability generating function of is 

k 



E 



z tT 



t t 

T + f Z 



(z > 0). 



(2.156) 



2.11 Poisson process and Poisson bridge 
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All moments of N^) exist, and the first two are 



E 



(k) 



tr 



E 



tT 



1 



t 

T 



(2.157) 
(2.158) 



We denote the jump times of {N^} by {T^}^ =1} and the waiting times by {W^}^. 
It is well known that the jump times are distributed as ordered, independent uniform 
random variables (see, for example, Sato 83] or Mikosch j74|). In particular, for Um < 
■ • ■ < Uq.) an ordered sequence of independent standard uniform random variables, we 
have 



T, 



k+l 



law 



(2.159) 
(2.160) 



T(U(i), . . . , U(k)), 

where {T«} are the jump times of the Poisson process {N t }. It is also well known that 
the jump time T k+ i is independent of the vector 



T 



T, 



fc+i 



(Tu-..,T k ). 



The distributional identity ( I2.159P is equivalent to the following: 



w, 



(k)\ law 



T 



(W 1 ,...,W k ), 



(2.161) 



(2.162) 



where {Wi} are the waiting times of {N t } (i.e. they are independent, identically- 
distributed exponential random variables). The ordered uniform random variable Uu\ 
has a beta distribution with parameters a = i and /3 = k — i + 1. Hence we have 

Qpf ] < t] = I t/T [i, k — i + l]. (2.163) 

Here I z [a,/3] is the regularized incomplete beta function which is defined by 



I z [a,P] 



yy-'dy" 



for a, (3 > 0. 



Define a process {N$} for < t < T by 



iV, 



(k) 



IT 



N(±T k+1 ) , N, 



(k) 

TT 



k. 



(2.164) 



(2.165) 



Then {N^} < t <T is a counting process that jumps at the times ^-j(Ti, . . . ,T k ), and 
is independent of T k+1 . It follows that 



{N t T }o<t<T — {N t T }o<t<T- 



(2.166) 



Chapter 3 

Levy random bridges 



The idea of information-based asset pricing is to model the flow of information in 
financial markets and hence to construct the market filtration explicitly. Let Xt be a 
random variable (a market factor), with a given a priori distribution. The value of Xt 
will be revealed to the market at time T. We wish to construct an information process 
{£tr} such that £tt = X T . We can then use the filtration generated by {^t} to model 
the information that market participants have about Xt- One problem to overcome is 
how to ensure that the marginal law of £tt is the a priori law of Xt- 

Two explicit forms for the information process have been considered in the litera- 
ture. The first is 

&t = ~*t + At (0 < t < T), (3.1) 
where \0tT}n< t<r is a Brownian bridge starting and ending at the value (see By 



55 



64 



18 



8 ll|). The second is 

^ T = X T -ftT (0<t<T), (3.2) 

where Xt > and \itT\o<t<T is a gamma bridge starting at the value and ending at 
the value 1 (see Q). These forms share the property that each is identical in law to 
a Levy process conditioned to have the a priori law of Xt at time T. The Brownian 
bridge information process is identical in law to a conditioned Brownian motion, and 
the gamma bridge information process is identical in law to a conditioned gamma 
process. 

With this as motivation, in this chapter we define a class of processes that we 
call Levy random bridges (LRBs). An LRB is identical in law to a Levy process 
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3.1 Defining LRBs 
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conditioned to have a prespecified marginal law at T. Later we shall use LRBs as 
information processes in information-based models. 

3.1 Defining LRBs 

An LRB can be described as a process whose bridge laws are Levy bridge laws. In 
the definitions below we define LRBs by reference to their finite-dimensional distri- 
butions rather than as conditioned Levy processes. This proves convenient in future 
calculations. 

Definition 3.1.1. We say that the process {L t T}o<t<T has the law LRBc([0, T], {ft}, v) 
if the following are satisfied: 

1. Ltt has marginal law v. 

2. There exists a Levy process {L t } £ C[0,T] such that L t has density ft(x) for all 

te(o,T]. 

3. v concentrates mass where fx{z) is positive and finite, i.e. < fr{z) < oo for 
v-a.e. z. 

4- For every n £ N +; every < ii < - • • < t n < T, every (x\, . . . , x n ) £ W 1 , and 
v-a.e. z, we have 

Q [L tl ,T <xi,..., L tnjT < x n | L TT = z] = Q [L tl <x 1 ,...,L tn < x n \L T = z}. 

Remark 3.1.2. Conditions 2 and 3 in Definition \ 3.1.1\ are sufficient conditions for 
the right-hand side of the equation in condition 4 to make sense. 

Definition 3.1.3. We say that the process {M tT }o<t<T has the lawLRBz>([0, T], {Qt}, P) 
if the following are satisfied: 

1. Mtt has probability mass function P. 

2. There exists a Levy process {M t } £ T> such that M t has marginal probability mass 
function Qt{a) for all t £ (0, T\. 

3. The law of Mtt is absolutely continuous with respect to the law of Mt, i.e. 

if P{a) > then Q T (a) > 0. 
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3.2 Finite-dimensional distributions 



4- For every n G N +; every < t± < • • • < t n < T, every {k\, . . . , k n ) G 7U 1 , and 
every b such that P{b) > 0, we have 

Q [M tltT = a kl ,..., M tn)T = a kn \ M TT — b] — Q [M tl = a kl , . . . , M tn = a kn \ M T = b] . 

Definition 3.1.4. For a fixed time s <T, if the law of the process {t] s +t}o<t<T-s is of 
the type LRBc([0, T — s], ■ , ■), resp. LRB-p([Q, T — s], • , • ), then we say that {r] t } s <t<T 
has law LRB c ([s,T], ■ , ■), resp. LRBx>([s,T], •, •). 

If the law of a process is one of the Li?i?-types defined above, then we say that it 
is a Levy random bridge (LRB). 



3.2 Finite-dimensional distributions 

For the rest of this chapter we assume that {L t r} and {M t r} are LRBs with laws 
LRBc([0, T], {ft}, v) and LRB?>([0,T], {Q t }, P), respectively. We also assume that 
{L t } is a Levy process such that L t has density ft(x) for t < T, and {M t } is a Levy 
process such that M t has probability mass function Qt(a^) for t < T. 
The finite-dimensional distributions of {L t r} are given by 

■n 

Q [L tuT G dxi, . . . , Lf n>T G dx n , L TT G dz] = J^J [f u _ u ^(xi - x^i) dxi] ipt n (dz; x n ), 

(3.3) 

where the (un-normalised) measure ip t {dz;^) is given by 

Mte;0 = v(dz), (3-4) 

^(d Z; o = /r "; ( "; e) ^), (3.5) 

for < t < T. Given f)3.3p . Kolmogorov's extension theorem ensures the existence of 
LRBs (up to a modification). It follows from the definition of LRBc([0, T], {ft}, v) and 
equation ( I2.25P that 

f tT {x; z) = —— (3.6) 

Jt{z) 

is a well-defined density (as a function of x) for t < T and u-a.e. z. Then from (13. 3p 
the marginal law of L t r is given by 

POO 

it G dx) = f t (x)^tO&; x)dx = / ftr(x; z) v{dz) dx. (3.7) 



3.2 Finite-dimensional distributions 



45 



Hence the density of L t r exists for t <T, and 

< i/j t (M.;x) < oo for Lebesgue-a.e. x G Support (/t). (3.8) 

In particular, we have 

< Vt(R; Ut) < oo and < f T -t(x - L tT ) < oo (3.9) 

for a.e. value of LtT- If v({z}) = 1 for some point z G R, i.e. Q[Ltt = z] = 1, then 
is a Levy bridge. If i/(cU) = fr(z) dz, then = {L t } for t G [0, T\. 

In the discrete case, the finite-dimensional probabilities of {M^} are 

n 

Q[M tl>T = a kl , . . . , M tniT = a kn ,M TT = z] = JJ [Qu-t^Mh -flfci-J] 0t„(^«fc„), 

(3.10) 

where the function </>t(z; ^) is given by 



»o(^0 = ^), (3.H) 
^;fl = gr -' ( '7 P(*), (3.12) 



for < t < T. If P is identical to Q T , then {M tT } l = {M t } for t G [0, T]. 

The existing literature on information-based asset pricing exploits special properties 



Brownian and gamma bridges. See Emery & Yor 33j for an insight into how remarkable 
these bridges are. The methods we use do not require special properties of particular 
Levy bridges. However, we will often use the Brownian and gamma cases as examples, 
and the results we obtain agree with previous work. 

Many of the results that follow are proved for the LRB {Ltr}, which has a contin- 
uous state-space. Analogous results are provided for the discrete state-space process 
{M tT }; details of proofs are omitted since they are similar to the continuous case. 
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3.3 LRBs as conditioned Levy processes 



3.3 LRBs as conditioned Levy processes 

It is useful to interpret an LRB as a Levy process conditioned to have a specified 
marginal law v at time T. Suppose that the random variable Z has law v, then: 

Q [L tl e dxi, ...,L tn G dx n , L T edz\L T = Z] 

= Q [L h G dxi, ...,L tn G dx n | L T = z] u(dz) 
fr-tn^i (z — Xn-l 



j [ [fu-u-ifa - Zi-i) dxi] v{dz). (3.13) 
Hence the conditioned Levy process has law LRB c ([0, T], {ft}, v). 

3.4 The Markov property 

In this section we show that LRBs are Markov processes. The Markov property is 
a key tool in the application of LRBs to information-based asset pricing. As will be 
seen below, the Markov property of an LRB follows from the Markov property of the 
associated Levy bridge processes. Note that if we stop the LRB {L t x} at time s > 0, 
and then restart it, the terminal law will no longer be the a priori law u, but instead 
an updated law conditional on L s t- 

3.4.1 Continuous state-space 

Proposition 3.4.1. The process {L tT } < t < T is a Markov process with transition law 



it G dy L sT = x\ = — — — rft- s {y - x) dy, 
ip s (R;x) 

orr c a i t i ^sidy-x) 

Q[L TT G dy | L sT = x\ 



(3.14) 



where < s < t < T. 

Proof. To show that {L t r} is Markov, it is sufficient to show that 

Q [L t T < y | L tl , T = xi, . . . , L tmtT = x m ] = Q [L tT < y \ L tm , T = x m ] , (3.15) 

for all m G N+, all {x u . . . , x m , y) G M m+1 , and all < t x < ■ ■ ■ < t m < t < T. When 
t = T we apply the Bayes theorem to (13. 3p and obtain 

Q [L TT G dy | L h , T = Xl , . . . , L tm , T = x m ) = ^ J ^ Xm l (3.16) 



3.4 The Markov property 
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We need now only consider the case t < T. Proposition 12.3.11 shows that Levy bridges 
are Markov processes; therefore, 

Q[L t <y\ L tl =x 1 ,...,L tm = x m ,L T = x] = <Q> [L t < y | L tm = x m , L T = x] . (3.17) 



It is straightforward by Definition 13.1.11 part HI to show that LRBs are Markov processes. 
Indeed we have: 

Q [L tT < V | L tuT = xi, . . . , L tm)T = x m ] 

/oo 
Q [L t T < V I L tl ,T = x u . . . , L tmtT = x m , L TjT = x] u(dx) 
-oo 

/oo 
Q[L t <y\ L tl = Xx, . . . , L tm = x m , L T = x] u(dx) 
-oo 

[L t <y\ L tm = x m , L T = x] u(dx) 



oo 
oo 



— oo 
oo 



Q [L tT < y I L tm:T = x m , L T:T = x] u(dx) 

— Q [LtT < y I L tmi T — x m ] . (3.18) 
The form of the transition law of {L^} appearing in (I3.14p follows from (13.31) . □ 



Example. In the Brownian case we set 

ft{z) - 



r7 hxt 



exp 



z 
2t 



(3.19) 



for t > 0. Thus ft{x) is the marginal density of a standard Brownian motion at time 
t. Then we have 



IT 



G dy | L sT = x] 



uydz) e 2 *- s 



T — t poo 

e 

J —CO 



1 \ (z-x)2 z 2 

2 T-s T 



idz) VMt-s) 



dy, (3.20) 



and 



Q[L T t e dy | L sT = x] 



v{dy) 



I: 



1 (z-x)- 1 z - 



2 T-s T 



] z/(d^) X 



^P-"' J z/(dy) 
- e r^[ xz -^ z2 } v(dz] 



. (3.21) 



Example. In the gamma case we consider the one-parameter family of processes 
indexed by m > described in Section 12.6.11 We set 

,mt 



m 



mt—l n —mx. 



(3.22) 
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so ft(z) = g[ m \x) for g^{x) given by f )2.46p . These densities are the increment 
densities of the gamma process with mean 1 and variance mT x at t — 1. Then 



Q[L tT E dy | L sT = x] 

1 -W ^-"^-"-'^^fe-^C-,-^, (3 ,3) 



B[m(T - t), m(t - s)] - a;)m(T--)-l z l-mT j,( dz ) 

and 

i li„>x>(y-a;) m(T_sKl y 1_mT ^(dw) 
^ £ ^ ' = = 'f(' - ^-.mAt • ^ 

where B[a,/5] is the Beta function. 
3.4.2 Discrete state-space 

The analogous result to Proposition 13.4.11 for the discrete case is provided below — the 
proof is similar. 

Proposition 3.4.2. The process {M t T}o<t<T has the Markov property, with transition 
probabilities given by 

X]/~L-oo ( t ) t{. a k] a j) 

Q [M t T = aj | M sT = en] = ^ °° zQt- s {aj - a*), 

Q [M TT = a, | M sT = ai } = =J^j^ r , 
where < s < t < T. 

□ 



3.5 Conditional terminal distributions 

Let {Ft} and {J 7 ^ 1 } be the nitrations generated by {L t f} and {M t j-}, respectively. 



s 

P s be the J 7 ^ 1 -conditional probability mass function of the terminal value Mtt- 



Definition 3.5.1. Let v s be the -conditional law of the terminal value Ltt, and let 

-M 



We have uq(B) = v(B), and Po(a) = P(a). Furthermore, when s > 0, it follows 
from the results of the previous section that 

^ ) = Q[ilTeB |^] = |||M, (3 , 6) 
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and 



P s (a k ) = Q [M TT = a k | Tf } = ^J^f^ y (3-27) 

When the a priori qth moment of Ltt is finite, the conditional gth moment is finite 
and given by 

/•oo 

\z\ q v s {dz). (3.28) 

Similarly, when the a priori gth moment of Mtt is finite, the J r j Vf -conditional gth 
moment is finite and given by 



W\ q P s {a k ). (3.29) 



k=— oo 



When they are finite, the quantities in (13.281) and (I3.29p are martingales with respect to 
{J^i} and {J 7 ^ 1 }, respectively. If q G Z then J \z\ q u(dz) < oo ensures that J z q i>(dz) 
is a martingale, and ^ |afc| 9 -P(a/c) < oo ensures that a q k P(a k ) is a martingale. 

When the terminal law v admits a density, we denote it by p(z), i.e. u(dz) = p(z) dz. 
In this case the LtT-conditional density of Ltt exists, and we denote it by 

v t (dz) f T -t(z - L tT )p(z) 

P,W = — = MXWMz) ' (3 ' 30) 



3.6 Measure changes 

This section leads to Proposition 13.6.11 which states that there exists a measure L 
equivalent to Q under which the LRB {Lit} is a Levy process. To begin the analysis 
it proves convenient to assume the existence of such a measure L, which we do, and 
we further assume that under L the density of L tT is ft(x). 

Writing ip t = ip t (M; L tT ) , we can show that {^}o<i<r is an L-martingale (with 
respect to the filtration generated by {L tT }). In particular, for times < s < t we 
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have 



fr-t( z — L tT ) 



u(dz) 



7 L 



Mz) 

fr-t{z — L s t — {L t T — L s t)) 



OO 

OO POO 



f T -t(z - L sT - y) 



u(dz) 



L 



sT 



y=—oo J z— — oo 
oo 1 r-oo 



Mz) 



v(dz) ft- s {y) dy 



POO j /*OC 

/ -r-TT / fr-t{z- L sT -y)f t „ s (y)dyv(dz) 

J Z = — OO JT\Z) Jy= — oo 



f T - s (z - L sT ) 

Mz) 



v{dz) 



(3.31) 



Since ipo — 1, we can define a probability measure L rb by the Radon- Nikodym derivative 

dL rb 



dL 



ip t for < t < T. 



(3.32) 



It was noted in Section !3Ti2l that < ip t < oo, so L rb is equivalent to L for t < T. For 
s, t satisfying < s < t < T, the transition law of {L t r} under L rb is 

L rb [L tT e dy | ^] = E Lrb [l {LtTedy} | ^] 

= ^ E L [lp t t{L tT €dy} | L sT ] 

fr-t(z - y) 



Mz) 



v(dz) ft- 8 (y - L sT ) dy 



L S T. 



ft-s{y — L s t) dy. 



(3.33) 



We see that {L tT } < t<T is a Markov process under the measure L rb . Furthermore, by 
virtue of Proposition 13.4.11 {L tT } is an LRB with law LRB c ([0, T], {ft}, v). 
We can restate this result with reference to the measure Q as the following: 



Proposition 3.6.1. Let L be defined by 

dL 



(3.34) 



* i 



fort G [0, T). Thenh is a probability measure. Under h, {L t T}o<t<T is a Levy process, 
and L t T has density ft(x). 



3.7 Dynamic consistency 
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□ 

In the case of a discrete state space a similar result is obtained: 
Proposition 3.6.2. Let L be defined by 

-l 

(3.35) 

for t G [0, T) . Then L is a probability measure. Under L, {Mtr}o<t<T is a Levy process, 
and M t T has mass function Qt(a). 

□ 

3.7 Dynamic consistency 

In this section we show that LRBs possess the so-called dynamic consistency property. 
For {L tT }, this property means the process {r) t } defined by setting 

rn = Ut -L sT (s<t< T) (3.36) 

is an LRB for fixed s and L s t given. Defining the filtration {F^} by 

F? = a (L sT , {r]u}s<u<t) , (3.37) 

we see that 

s<u<T) \F t L ], (3.38) 

for < s < t < T and F an arbitrary measurable functional. Suppose two market 
participants, trader A and trader B, watch the evolution of {L tT }; trader A watching 
from t — and trader B watching from t — s. The filtration of trader A, {F^}, is 
larger than the filtration of trader B, {F?}, but they have a common view of the future 
evolution of {L t r}- This is the Markov property. The dynamic consistency property is 
stronger. It states that the filtration of trader B can be regarded as being generated by 
an LRB, in this case {r) t }, plus some information about the current state of the world, 
in this case L s t- The a priori terminal law of {i] t } will not be v, but instead it will be 
an updated measure that depends on the current state L s t- 

Later we shall model the market filtration as being generated by a set of LRBs. 
Through the dynamic consistency property, we can consider each market participant's 



dL 



ttM 



<Pt{a k -M tT ) 
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filtration to be generated by a set of LRBs, regardless of the time in which they enter 
the market, and without their views being inconsistent with other participants. 

The dynamic consistency property was introduced in Brody et al. 17| with regard 



to Brownian random 



random bridges in 20] 



Dridges, and was shown by the same authors to hold for gamma 



Fix a time s < T. Given L s t, we define a process {r] t } by (13.361) . We shall show 
that {rj t } is an LRB. At time s, the law of rjr is 

u*(A) = u s (A + L sT ) for all A G B(R), (3.39) 

where A + y denotes the shifted set 

A + y = {x : x -y G A} . (3.40) 

Given the terminal value rjx, the finite-dimensional distributions of {r] t } are given by 

Q [r] s+tl G dx 1 , . . . , i] s+tn G dx n | L sT , r] T = z] 

= Q [L s +ti,T — L sT G d^i, . . . , L s+tniT — L sT G dx n \ L sT , L TT — L sT = z] 
= Q [L s+h - L s G dxi, . . . , L s+tn - L s G dx n \ L s , L T - L s = z] 
= Q [L tl G dxi, . . . , L tn G dx n | L r _ s = z] 

= fT ~ S : tn { ) = Xn) ]I fu-u-, {x t - x^) , (3.41) 

for every n G N+, every = t < t± < ■ ■ ■ < t n < T — s, and every (xi, . . . , x n ) G M n , 
where Xq = 0. Then we have 

Q G dxi , . . . , r/ s+in G dx n , r] T edz\ L sT ] 

= fT ~ S - tn {Z ~ Xn) Ufu-t^ fa-x^) v*{dz). (3.42) 

Comparison of this expression to (13. 3p shows that the process {i]s+t}o<t<T-s has the 
law LRB c ({0, T - s], {f t }, v*), and so the law of {^} s < t < T is LRB c ([s, T], {f t }, v*). 
In the discrete case, we define {q t } by setting 

Vt = M tT -M sT (s<t<T). (3.43) 

Then, given M sT , {rj t } has the law LRBv([s,T], {Q t },P*), where P* is defined by 

P*(a) = P s (a + M sT ). (3.44) 



3.8 Increments of LRBs 
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3.8 Increments of LRBs 

The form of the transition law in Prop osit ion 13 . 4 . ll shows that in general the increments 
of an LRB are not independent. The special cases of LRBs with independent increments 
are discussed later. A result that holds for all LRBs is that they have stationary 
increments: 

Proposition 3.8.1. For s, t, u satisfying < s < u < T and < t < T — u, we have 
Q [L u+t>T -L uT <z\ L sT ] = Q[L s+t:T -L sT <z\ L sT ), (3.45) 

and 

Q [M u+t , T — M uT < z | M sT ] = Q[M s+t , T -M sT <z\M sT }. (3.46) 

Proof. We provide the proof for {L t x}- The proof for {M t r} is similar. Throughout 
the proof we assume that t < T — u. The case t = T — u follows from the stochastic 
continuity of {L t rp}. First we consider the case s = 0. It follows from (13.141) that 



: u+t , T e dy, L uT e dx] = ipu+t(^; y)ft(y - x)f u (x) dxdy. (3.47) 
Then we have 

j u+ t, T - L uT e dz, L uT e dx] = tpu+tiM; z + x)f t {z)f u (x) dx dz 

dw f t (z) f u (x) dxdz. 



f T -(u+t){w -Z-X) 



fr(w) 

(3.48) 



Integrating over x, and changing the order of integration yields 

f°° f°° dw 
j u+t,T ~ L uT Edz}= / / T _ (u+t) (w - z - x)f u (x) dx f t (z) dz. 

J w=— oo J x=— oo /TV / 

f T -t(w - z) 



w=—oo 



dwf t (z)dz 

Jt{w) 
i> t (R,z)f t (z) dz 

tredz]. (3.49) 



For the case s > 0, we use the dynamic consistency property. For s fixed and 
L sT given, the process {7] u t}s<u<t = {L u t — L sT } s < u < T is an LRB with the law 
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LRB c ([s,T}, {f t }, v*), where v*{A) = v s {A + L sT ). We have 

Q [L u +t,T - L uT G dz | L sT ] = Q [r] u+tt T - VuT G dz | L sT ] 

= Q[7/tr G dz\L sT ] 
f T -t(w - z) 



h-s{w 
f T -t(w - z + L sT 



f T - s (w - L sT 



u*(dw) f t - s (z) dz 

Vs(dw) fts(z) dz 



1 



f T -t{w -z + L sT ) 



t/j t (R; z + L sT 
tfj s {R; L sT ) 

Q[L tT - L sT e dz | L sT 



fr{w) 
fts(z) dz 



u(dw) fts(z) dz 



(3.50) 
□ 



When {Ltr} is integrable, the stationary increments property offers enough struc- 
ture to allow the calculation of the expected value of L t r'- 



Corollary 3.8.2. J/E[|L tT |] < oo for allt G (0,T] then 

E M L, T ) r _ s _„ , T _ s 
and ifE[\M tT \] < oo for all t G (0,T] then 



T — t t — s 

L s t + m E [Ltt I L s t\ 



E [Mfr | M sT ] = - — -M sT + - — -E [M TT | M sT ] 



[s<t), 



[s<t). 



(3.51) 



(3.52) 



T-s T-s 

Proof. We only provide the proof for {L t x} since the proof for {M t j-} is similar. The 
case t = T is immediate, so we assume that t < T. First we consider the case s = 0. 
Suppose that t = —T, where m, n G N + and m < n. We wish to show that 



E[L tT ] = -E[L TT \. 



n 



(3.53) 



Writing L(t,T) = L t ? for clarity, define the random variables {Aj} by 

A l = L (iT, T)-L (^T, t) . (3.54) 

It follows from Proposition 13.8.11 that the Aj's are identically distributed, and by as- 
sumption they are integrable. Hence we have 



E[Ai] = — E 
n 



i=l 



-E\L TT ] 
n 



(3.55) 
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Then, as required, we have 

E[L(fT,T)]=E 



1=1 



m 



n 



E[L 



TT\ 



(3.56) 



For general t, choose an increasing sequence of positive rational numbers {q{\ such 
that lim^oo qi = |j. From (I3.56P we have 



E[L(t, T)] = E[L( gi T, T)] + E[L(t, T) - L( gi T, T)) 
= qi E[L(T, T)] + E[L(t, T) - L( qi T,T)\. 



(3.57) 



By stochastic continuity, in the limit i — > oo the law of L(t,T) — L(g,T, T) tends to a 
Dirac measure centred at 0. Hence 



E[L(t, T)) = lim qi E[L(T, T)} + E[L(t, T) — L( qi T, T)} 

i—^oo 

= |e[L(T,T)]. 



(3.58) 



For the case where s > 0, we use the dynamic consistency property. For s fixed 
and L s t given, the process 

VtT = L tT -L sT {s<t< T) (3.59) 
is an LRB with law LRBc{[s, T], {ft}, v*), where u*(A) = v s {A + L s t)- Then we have 



E [L tT \L s t} = L sT + E[r] tT \ L sT ] 

t-s ro ° 



L sT + 
L sT + 



T-s 
t - s 
T-s 



zu*(dz) 

[z - L sT ) u s (dz) 



T — s t — s 

-L s t + m IE [Ltt I L s t\ 



T-s 



T-s 



(3.60) 
□ 



We have shown that the increments of LRBs are stationary, and so it is then natural 
to ask when the increments are independent, i.e. when is an LRB a Levy process? The 
answer lies in the functional form of ^(M; y). 

For < s < t < T, the likelihood that L t T = y given that L s t = x is 



q{t > y]S > x) = W^) ft - s{y ~ X) - 



(3.61) 
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If {L t r} has stationary, independent increments then 



q(t, y; s, x) = q(t-s,y- x\ 0, 0). 



(3.62) 



Therefore the ratio 



^ S (K; x) 

is a function of the differences t — s and y — x. Thus, if we have that 



(3.63) 



aexp(by + ct), 



(3.64) 



for constants a, b and c, then {Ltr} is a Levy process. There are constraints on a, b 
and c since (I3.6ip is a probability density. When b = c = we have v(dz) = fr( z ) 
which is the case where {L t x} == {L t }. Note that the subclass of LRBs that are Levy 
processes is small. 

Example. In the Brownian case we consider a process {W^t} with law 



where ft{x) is the normal density with zero mean and variance t given by (13.191) . In 
other words, {W 7 ^} is a standard Brownian motion conditioned so that Wtt is a normal 
random variable with mean 9T and variance T. In this case, we have 



Simplifying the expression for the transition densities of the process {W^t} allows one 
to verify that {Wtr} is a Brownian motion with drift 9. It is notable, by Girsanov's 
theorem, that the process {i/) t (M.; Wt)} is the Radon-Nikodym density process that 
transforms a standard Brownian motion into a Brownian motion with drift 9. Hence 
we can alternatively deduce that {W^t} is a Brownian motion with drift 9 from the 
analysis in Section I3T61 

Example. In the gamma case, we consider a process {T^} with law 



LRB c ([0,T],{f t }J T (z-9T) dz) 




(3.65) 



LRB c {[^Tl{f t },K- 1 f T {z/ K)dz), 
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where ft(x) = g ( [ a \x) is the gamma density with mean t and variance t/m defined by 
( I2.46p . and k > is constant. Then {T^} is a gamma process with mean unity and 
variance rrT 1 at t — 1, conditioned so that IYt has a gamma distribution with mean 
kT and variance K 2 T/m. We have: 

J-oc fAz) K 
= n- mt exp (m(l - K- l )y) . (3.66) 

The transition density of {r t T} is then 

(y _ x \m(t-s)-l e -m(y-x)/K 

Q\T tT e dy I T sT = x} = lt v>x \ , . u ^ r , rr- dy. (3.67) 

Hence {T^} is a gamma process with mean k and variance k 2 jm at i = 1. 

3.8.1 Increment distributions 

Partition the time interval [0, T] by = to < ti < t<i < ■ ■ ■ < t n = T. Then define the 
increments {Aj}™ =1 and {aj}™ =1 by 

A, = L tuT - L ti _ uT (3.68) 
Qij = ti — tj-i- (3.69) 



Assume that v has no continuous singular part 83[. Denoting the Dirac delta function 
centred at z by 5 z (x), x G R, we can then write 

oo 

u{dz) = J2 Vi8 Zi (z)dz+p(z)dz, (3.70) 

i=— oo 

for some {cij} C R, {z;} C R+, and p : R — > R + . Here is the density of the 

continuous part of u, and is a point mass of v located at Zi. From (13. 3p . the joint 
law of the random vector (A 1; . . . , A ra ) T is given by 



n \ n 



Q[Ai G dyi . . . , A n G dy n ] = / IJ^yi II /«*^*) ( 3 - 71 ) 

\i=l / i=l 

where 

/(*) = m • (3.72) 
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Equation (13.7ip shows that (Ai, . . . , A n ) T has a generalized multivariate Liouville 



distribution as defined by Gupta & Richards 51] . The classical multivariate Lion 



vil 
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e distribution is obtained when ft(x) is the density of a gamma distribution (see 
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_ .36J). A survey of Liouville distributions can be found in Gupta & Richards 
47]. Barndorff-Nielsen & J0rgensen [9] construct a generalized Liouville distribution 
by conditioning a vector of independent inverse- Gaussian random variables on their 
sum. 

In the discrete case, the joint distribution of increments also has a generalized 
Liouville distribution. Define the increments {D{\ by 

A = M UtT - M ti _ uT . (3.73) 

Then we can write 

/ n \ n 

Q[£>! ed Vl ...,D n e dy n ] = Q ( J> ) HdQM, (3.74) 

where 



j=i / i=i 



V°° P( ai )5 a .(z) 
Q{z) = ^ . (3.75) 

3.8.2 The reordering of increments 

We are able to extend the Markov property of LRBs. If we partition the path of an 
LRB into increments, then the Markov property means that future increments depend 
on the past only through the sum of past increments. We will show that for LRBs the 
ordering of the increments does not matter for this to hold — given the values of any set 
of increments of an LRB (past or future), the other increments depend on this subset 
only through the sum of its elements. 

Let Ti be a permutation of {1,2, ... ,n}. We define the partial sum by 

m 

S m = ^ A tt(») for m = 1,2, . . . ,n, (3.76) 

where the {Aj} are defined as in f)3.68p . and we define the partition = tg < t\ < 
■ ■ ■ < tl = T by 

3 

*J+i = ^<M<) for j = 1,2,..., Ti-l. (3.77) 

i=i 
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Proposition 3.8.3. We may extend the Markov property of {L t r} to the following: 



A^m+i) < y m +i, . . . , A„.( n ) < y n | A 7r(1) , . . . , A^^j ] = 

A^( m +i) < y m +i, ■ 



.,K(.n)<y n \Sl]. (3.78) 



If v has no singular continuous part, then 



A^m+i) G dy m+1 , . . . , A,r( n ) G dy n | 5^] = 



i=rrt+l 



Proof. Define the increments {A J 1 "} by 



A- — -Z^+7r t -1 — L+tz nr. 
t l n< 1 l n-l> 1 



(3.80) 



The law of the random vector (Aj , . . . , A^_ x , Af ) T is given by 



Ai G dyi , . . . , A-_ x G dj/ n _i, ]T Af G dz 



i=i 



i/(ds) 



n—1 \ n—1 



i=i / i=i 

This is also the law of (A 7r (i), . . . , A^(„_i), A 7r (j)) T , hence 



(A. (1)) ...,A. (n) )^(A-,...,A-). 



(3.82) 



The Markov property of LRBs gives 



A^ +1 < y m+1 ,...,Al<y n \A*,...,AZ l ] 



A^ +1 < y m +i, ...,A*<y n 



i=l 



(3.83) 



and so we have 



A^m+i) < y m +i, ■ ■ ■ , A^ < y n | A 7r(1 ), . . . , A,r( m ) ] = 

A,r( m+ i) < y-m+l, ■ 



.,A An) <y n \Sl]. (3.84) 



This proves the first part of the proposition. 
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For the second part of the proof we assume that v takes the form f)3.70p . Note that 
= Y^hLi ^Ht an d that the density of is 

x^f t ,Jx)^jR;x)= [°° k^)h-^-x) vm {3 g5) 

Jz=-oo Jt{Z) 

The elements of the vector (L^t, ^-m+ii • • • 5 A n) T are non-overlapping increments of 
{Ljt}, and the law of the vector is given by 

Q [L t ^T G dx, A n m+1 G dy m+1 , . . . , A£ e dy n ] = 

/ n \ n 

f\ x+ yi)ft7n( x ) dx [ I fa v(i) {yi) dyi. (3.86) 



i=m+l / i=m+l 



Thus we have 



[ A m+i e dy m+1 , . . . , A£ G dy n I L t ^ T = x] 

_ g [A* +1 G dy m+1 , . . . , A% G dy ra , L^t G dx] 
Q [L t ^ T G dx] 

/ (x + X]j=m+1 2/*) rii=m+l /«*•(<) (l/i) 



(3.87) 



□ 



We note that Gupta & Richards [5jJ prove that if (A 1; A 2 , . . . , A n ) T has a general- 



ized Liouville distribution then equation ( I3.78P holds. 

We can use Proposition 13.8.31 to extend the dynamic consistency property. In par- 
ticular we have the following: 

Corollary 3.8.4. a. Fix times si,Ti satisfying < T± < T — s±. The time- shifted, 
space-shifted partial process 

rfil = L Sl+ttT - L SuT , {0<t< Tx), (3.88) 

is an LRB with the law LRBc([0,Ti],{f t },u^), where is a probability law 
on R with density fa (x)ipTi (Mi x )- 

b. Construct partial processes {^y.}, i = 1, • • • , n, from non- overlapping portions of 
{L tT }, in a similar way to that above. The intervals [sj, Sj + Tj], i = 1, . . . , n, are 
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non- overlapping except possibly at the endpoints. Set rj^j,. = rjj, T . when t > TJ. 
Ifu>t, then 



(1) _ (1) (") _ (n) 

^u.Ti %,Ti — Xl ' • • • ' ^w,T n "t,T„ — x r. 



TV 



(1) _ (1) (n) _ in) - 

Vu,Ti ^i,Ti — X l' • • • ' Vu,T n Vt,T n — X r 



i=l 



where 



(3.89) 



(3.90) 



Remark 3.8.5. 27ie partial processes of Corollary\3.8.4\ are dependent, and we have 



7)$ G da; 



I]® G da; 



(3.91) 



for < s < t < T. 



We state but do not prove a discrete analogue of Proposition 13.8.3^ which is as 
follows: 

Proposition 3.8.6. We may extend the Markov property of {M tT } to the following: 

Q [Ar(m+1) < Vm+1, • • • , Ar(n) < 2/n | Ar(l)> • • • ; Ar(m) ] = 

Q [D <m+1) < y m+1 , D <n) < y n \ R^] , (3.92) 
where R 7 ^ = YlT=i Ar(?) ■ Furthermore, 

«r n n I 1 ^ ("^ m Si=m+1 Vi) TT ^ / \ 

Q [D n{m+1) = y m+1 , D <n) = y n \ D m \ = — ^ m ^ [[ Q a ^ ) {y l ). 



fc=_oo^(afc;-Rm 



i=m+l 



(3.93) 



Corollary 13.8.41 can be extended to include LRBs with discrete state-spaces. 



Chapter 4 



Information-based asset pricing 



We apply LRBs to the modelling of information flow within the information-based 
framework of Brody, Hughston & Macrina (BHM). The approach was applied to credit 
risk in Brody et al. 17| . and this was extended to include stochastic interest rates in 
'Jlutkowski & Yu 8l|. A general asset pricing framework was proposed in Brody et al. 



191 ] (see also Macrina |64h, and there have also been applications to inflation modelling 
Hughston & Macrina |55j), insider trading (Brody et al. [18]). insurance (Brody et al. 
20]), and interest rate theory (Hughston & Macrina jssj] ) . 

We model a financial market as a collection of cash flows occurring on fixed dates. 
We assume that each cash flow can be expressed as a function of independent market 
factors, which we call X-factors. Under the assumptions that interest rates are deter- 
ministic and that the pricing kernel is given, the no-arbitrage price of a cash flow is its 
discounted expected value under the risk-neutral measure. 

Each X-factor is revealed to the market through a factor information process. We 
model a factor information process as an LRB whose terminal value is the value of an 
X-factor. The market information process is then a collection of all factor information 
processes, and this process generates the market filtration. Hence, we are modelling the 
information flow in the market by our choice of X-factors and factor information pro- 
cesses. The dynamics of cash flow prices are derived by taking conditional expectations 
with respect to the market filtration. 

Stylistic path properties of price processes can be incorporated into the model by an 
appropriate choice of factor information processes. For example, if a cash flow depends 
on an X-factor whose information process is a pure-jump LRB, then the cash flow's 
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price process will be a pure-jump process. Dependence between two cash flows can be 
modelled by expressing them as functions of a common set of X-factors. 

To simplify matters, we examine the case of a cash flow that pays an amount equal 
to the value of a single X-factor. We derive the price process for the cash flow, and 
we derive an expression for a call option on this price. Although this model is simple 
in structure, because the class of information (LRB) processes is large, the results are 
quite general. In the special case that the X-factor can take only two values, we recover 
a generalisation of the binary bond model of Brody et al. Vj\ . All the results of this 
chapter are for a general LRB information process. We consider specific examples in 
the subsequent chapters. 



4.1 BHM framework 

We fix a finite time horizon [0,T] and a probability space (Q, J 7 , Q). We assume that 
the risk-free rate of interest {r t } is deterministic, and that r t > and f°°r u du = oo, 
for all t > 0. Then the time-s (no-arbitrage) price of a risk-free, zero-coupon bond 
maturing at time t (paying a nominal amount of unity) is 

Pst = exp (- r u du^j (s<t). (4.1) 

For t < T, the time-t price of an integrable contingent cash flow H?, due at time T, is 
given by an expression of the form 

H tT = P tT E[H T \T t ), (4.2) 

where {J-'t} is the market filtration. The sigma-algebra T t represents all the information 
available to market participants at time t. In order for equation (14. 2 j) to be consistent 
with the theory of no-arbitrage pricing, we must interpret Q to be the risk-neutral 
measure. 

In such a set-up, the dynamics of the price process {H t j<} are implicitly determined 
by the evolution of the market filtration {J 7 *}. We assume the existence of a (possibly 
multi-dimensional) information process {^tr}o<t<T such that 



J~t — °~ ({£st}o<s<*) 



(4.3) 
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So {£tr} is responsible for the delivery of all information to the market participants. 
The task of modelling the emergence of information in the market is reduced to that 
of specifying the law of the information process {£tr}- 

4.1.1 Single X-factor market 

We assume that the cash flow H T can be written in the form 

H T = h(X T ), (4.4) 

for some function h(x), and some market factor X T . We call X T an X-factor. We 
assume that {C,tx} is a one-dimensional process such that £tt = X T . Then we have 

H tT = P tT E[h(X T ) | T t ] = P tT E[/i(6rr) | F t \, (4.5) 

which ensures that H TT = H T . In the case where {£<t} is a Markov process, we have 

H tT = P tT E[h(i TT )\i tT \. (4.6) 

4.1.2 Multiple X-factor market 

In the more general framework, we model a financial asset that generates the N cash 
flows H Tl ,H T2 , . . . , H Tn) which are to be received on the dates 7\ < T 2 < • • • < T N , 
respectively. At time T k) we assume that the vector of X-factors X Tk € M" fe (n k G N+) 
is revealed to the market, and we write 

X Tk =(x%,X%,...,X^) T . (4.7) 

We assume the X-factors are mutually independent, and that 

Hr k = h k (X Tl , X Ta , . . . , X Tk ) , (4.8) 

for some h k : R™ 1 x R™ 2 x • • • x R nfc — > R which we call a cash-flow function. For each X- 
factor Xy.\ there is a factor information process J) } such that g J) = X^ for t > Tj, 
and the factor information processes are mutually independent. Setting T = T/v, we 
define the market information process {£*t} to be an R ni+n2+ '" +nAr -valued process with 
each of its elements being a factor information process. The market filtration {Ft} is 
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generated by {£tr}- By construction, FiT k is J-j-measurable for t > T k . The time-t 
price of the cash flow Hx k is 

(k) ( P t , Tk E[h k (X Tl ,X T2 ,...,X Th )\F t ] iort<T k , 
H t = < (4.9) 
I for t > T k . 

Here we adopt the the convention that cash flows have nil value at the time that they 
are due. In other words, prices are quoted on an ex-dividend basis. In this way the 
price process {H^} is right-continuous at t — T k . The asset price process is then 

n 

H tT ^J2 H t k) (0<t<T). (4.10) 

k=l 

4.2 Levy bridge information 

We consider a market with a single factor, which we denote Xt- This X-factor is the 
size of a contingent cash flow to be received at time T > 0, so we take h(x) = x. For 
example, Xt could be the redemption amount of a credit risky bond. Xt is assumed to 
be integrable and to have the a priori probability law v (we also exclude the case where 
Xt is constant). Information is supplied to the market by an information process {£,tr}- 
The law of {£tr} is LRB c ([0, T], {ft}, v), and we set £tt = X T . We assume throughout 
this chapter that the information process has a continuous state-space; the results can 
be extended to include LRB information processes with discrete state-spaces. Indeed, 
in Chapter [9] we give a detailed example of a discrete LRB, and apply it to the pricing 
of credit derivatives. 

Since the information process has the Markov property, the price of the cash flow 
Xt is given by 

X tT = P tT E [X T I &r] (0 < t < T). (4.11) 

We note that X T is J-^-measurable and X T t = X T) but X T is not J-'j-measurable for 
t < T since we have excluded the case where Xt is constant. For t G (0,T), the 
^-conditional law of Xt as given by equation (I3.26P is 

"< (d2) = *»fap (4 - 12) 

where 

Vt(d*;0= fT 'f Z ~° dz. (4.13) 

JT{Z) 
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Then we have 

/oo 
zu t (dz). (4.14) 
-oo 

When v admits a density p(z), the J-j- conditional density of X? exists and is given by 

fr-t(z - &t)p(z) 



Pt{z) 



(4.15) 



Example. In the Brownian case the price is 



X, 



tr — *tr~ 



The following SDE can be derived for {X tT } (see [ill, [jij 64 



(4.16) 



Ml) 



dX tT = r 4 X tT dt H ; dW t , 



T-t 



(4.17) 



where {Wt} is an {J-" t }-Brownian motion. 



Example. In the gamma case we have 

JZiz-^^-'z^HdzY 



XtT — PtT 



(4.18) 



4.3 European option pricing 

We consider the problem of pricing a European option on the price X t T at time t. For 
a strike price K and < s < t < T, the time-s price of a t-maturity call option on X t j- 
is 

C st = P st E [(Xtr - K) + | U] • (4.19) 
The expectation can be expanded in the form 

E Q [(X tT - K) + | = Eq [{PtrE[X T | 6t] - ^) + | U] 



Er 



/ {P tT z-K)v t {&z)\ 

J — oo / 



(4.20) 
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Recall that the Radon-Nikodym density process 

dL 



(4.21) 



defines a measure L under which {£,tT}o<t<T is a Levy process. By changing measure, 
we find that the expectation is 



1 



E- 



( f (PtTZ-K)i/; t (dz;CtT 

\J — oo 

(P tT z - K) 



■Eir 



fr-tjz — 6t) 

Mz) 



v(dz) 



. (4.22) 



Equation f)3.9p states that < fr-s( z — £st) < oo, and we have 

U s (dz) = n^ZsT)- l fT -j Z ~^ T) U(dz). 

JT{Z) 



(4.23) 



Thus we can write the expectation in terms of the ^-conditional terminal law v s in 
the form 



{P tT z - K)- -— u s {dz i 

j JT-s{Z - t, sT ) 



sT 



/oo / /»oo 
/ (P tT z-K)- 
-oo \J — oo 



f T -t{z - x) 



f T -s(z ~ £ S T 

We defined the (marginal) Levy bridge density ftr{x] z) by 

f T -t(z - x)f t (x) 



v s (dz)\ f t „ s {x-U)dx. (4.24) 



ftrix; z) 



Mz) 



From this we can define the £ S T-dependent law fi s t(dx; z) by 

jJ, st {dx] Z) = f t -a,T-s{x ~ £sT, Z - £ sT ) dx. 



(4.25) 



(4.26) 



Thus fi s t(dx; z) is the time-t marginal law of a Levy bridge starting at the value £ s t at 
time s, and terminating at the value z at time T. Defining the set B t by 



B, 



x G 



x e 



— oo 
oo 



(P tT z - K) l T ~f f v s {dz) > 

{P tT z - K) h -f Z ~ x) v{dz) > o) , (4.27) 
Jt{z) J 
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where the second equality follows from (I4.23p . the expectation reduces to 

/oo 

(P tT z-K)^ st {B t - Z )v s (dz). (4.28) 
-oo 

Then the option price is given by 

/oo 

(P tT z - K)fi st (B t ; z) u s (dz). (4.29) 
-oo 

We can write X t x = A(t,C, t T), for A a deterministic function. The set B t can then 
be written 

5 ( = {eel: A(U) >K}. (4.30) 

We see that if A is increasing in its second argument then B t = (£ f *, oo) for some critical 
value Q of the information process. A is monotonic if the {£*t} is a Levy process. 



Example. In the Brownian case we have 



A(t,x) = P tT 



J —OO 



u(dz) 



IJU2 



/OO — i — XZ 
QT-t I 2 T 
— OO 



u(dz) 



(4.31) 



It can be shown that the function A is increasing in its second argument (see 

,Q); 

hence B t = (£ t *,oo) for the unique £ 4 * satisfying A(i, £*) = K. A short calculation 
verifies that fi st (dx; z) is the normal law with mean M(z) and variance V given by 

T-t 



M(z) 



t-s 

-£sT + ~ Z, 



V 



t-s 



(T-t). 



(4.32) 



T-s T-s T-s 

This is the time-t marginal law of a Brownian bridge starting from the value £ s t at 
time s, and finishing at the value z at time T. We have 



Vst{B t ;z) = 1 - $ 



£ - M(z) 



M (z) - Q 



(4.33) 



where is the standard normal distribution function. The option price is then 



C 



sT 



^(dz)+P s ^ / $ 



M(z) - c 



u s (dz). (4.34) 



Example. In the gamma case we have 



A(t,x) = P tT 



J™(Z - x )m(T-t)-l z 2-mT 

J™(Z ~ X )m(T-t)-l z l-mT u ^ " 



(4.35) 



4.4 Binary bond 
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The monotonicity of A(t,x) in x is proved for m(T — t) > 1 by Brody et al. 20 ]. 
The authors also give a numerical example where A(t,x) was not monotonic in x for 
m(T - t) < 1. For all t e (0, T), we have 



// st (dx; z) = l{e, r<!B <*} 
where fc(z) is the normalising constant 



x m(t-s)-l / \ m(T-t)-l 

-1 I x ~ Sat ^ f Z ~V 



sT 



dx, (4.36) 



fc(z) = (z-£ sT )B[m(*-s),m(T-*)]. 



(4.37) 



So /i st (dx; z) is an (z — £ s <r)-scaled, £ S T-shifted beta-law with parameters a = m(t — s) 
and (3 = m(T — t). This is the time-t marginal law of a gamma bridge starting at the 
value £ s t at time s, and terminating at the value x at time T. When m(T — t) > 1, a 
critical £ t * exists such that A(t, = K. Then £? t = (£ t *, oo), and 



Atat(5 t ; z) = 1 - J 
= I 



sT 



m(£ — s),m(T — t) 



z-Ct 



z-C 



; m(T — t), m(t — s) 



sT 



(4.38) 



where I[z\ a, f}\ = I z [a, (3] is the regularized incomplete beta function. The option price 
is then given by 



Cst — PsT / Zl 



z-Q 



Z - isT 

+ PstK 



; m{T — t), m(t — s) 
z-Ct 



Z - £sT 



u s (dz) 

■; m(T — t),m(t — s) 



v s {dz). (4.39) 



4.4 Binary bond 

The simplest non-trivial contingent cash flow is Xt € {ko,ki}, for ko < k%. This is 
the pay-off from a zero-coupon, credit-risky bond that has principal k\, and a fixed 
recovery rate k /ki on default. Assume that, a priori, Q[Xx = k ] = p > and 
tx — kx] — 1 — p. Then 



Q[X T = k | 6t] 



/r(fco) /r-tQi -&r) 1 ~P 
/r(fci) /r-tOo - 6t) P 



-i 



(4.40) 
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and 



tr\ 



1 + 



/r(fci) /r-t(fco - ftr) P 



-i 



(4.41) 



/r(fco) /r-t(fci - &r) 1 -p, 
The bond price process {X t x} associated with the given terminal cash flow is given by 

X tT = P tT (k Q[X T = k I & T ] + h Q[X T = h | & T ]) (0 < t < T). (4.42) 



Example. In the Brownian case we have 
Q[X T = k | &r] = (l + ex P 

and 



I — p 

V 



-i 



(4.43) 



Q[X T = h | & T ] = 1 + exp 



1 fci - fc t 

2 T-t 1t 



(fc + fci) - 26t) 



P 



1 — J9 



(4.44) 



Writing pj = Q[X T = ki \ fa], n °te that 

Var[X T | fa] = (ki - A; ) 2 piPo 

= -(fc - k p - fcipi)(fci - A; Po - fciPi) 
= -(A; -X tr )(A; 1 -X tT ). 

Thus, recalling f !4.1T|> . we see that the SDE of {X t r} is 

AV v A+ PtAh - X tT )(h - X tT ) 
dX tT = r t X tT dt i dW *> 



(4.45) 



(4.46) 



with the initial condition X 0T = P 0T (k p + kx(l — p)). For K G (Ptrko, Ptrki), we are 
able to solve the equation A(t, x) = K for x. We have 

A{t, x) = P tT (h Q[X T = k | fa = x} + h Q[X T = k\ \ fa = x]) 

= P tT (h - (hi - fco) Q[X T = k \fa = x]), (4.47) 

so the solution to A(t, x) = K is 

p K - P tT k Q 



The price of a call option on X t T is 



C st = P st J2(PtTk-K)$ 



i=0 



1-pPtrh-K 

M ^-^®lx T = ki \t sT ) 



V 



(4.48) 



(4.49) 



Chapter 5 



Variance-gamma information 



We presented some of the properties of the VG process and the VG bridge in Section 
12.71 We now consider the (standard) variance gamma random bridge (VGRB). 

In a similar way to a VG bridge, we have two terminal value decompositions of 
a VGRB. The first decomposition writes a VGRB in terms of its terminal value, a 
Brownian bridge, a gamma bridge, and a random volatility factor. The second de- 
composition writes a VGRB in terms of its terminal value, two gamma bridges, and a 
random volatility factor. These lead to two efficient algorithms for the simulation of 
VGRBs. 

We demonstrate the flexibility of VG information models by way of two simple 
examples. The first example recovers the equity model of Madan et al. 67]. In this 
case we consider a single V-factor Xt which, after a linear transformation, represents 
the log-return on a non-dividend-paying stock over the time period [0, T\. We assume 
that, a priori^ Xt is an asymmetric VG random variable, and that the market filtration 
is generated by a standard VGRB whose terminal value is pinned to the value of Xt- 
Under these conditions, we show that the stock price process is an exponentiated 
asymmetric VG process. 

The second example is an application to credit risk. We price a binary bond when 
the market filtration is generated by a VGRB. Following Brody et al. 17| . we include 
a rate parameter a which controls the speed at which information is released. Through 
simulations, we explore the effect of o on the sample paths of bond price processes in 
the case of default. When a is high, information about the impending default enters 
the market quickly, and the price of the bond becomes small well before the maturity 
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date. When a is low, information about the default is slow to reach the market, 
and there is a sudden collapse in the bond price close to the maturity date. Further 
simulations are performed that demonstrate the effect of the VG parameter m on bond 
price trajectories, where m _1 was the variance of the gamma time change at t = 1. 
When m is small, bond price processes exhibit more large jumps. When m is large, 
bond price processes are smoother, and look more like the prices generated in [17] using 
Brownian information. 



5.1 Variance- gamma random bridge 

Let {Vtx} be a standard VGRB, so in the transition law (13.141) we take ft{x) = ft m \%) 
which we defined in (12.731) to be 

iiSfe)""'^*^]- <«> 

for m > 0. We assume that Vtt has marginal law v. To ensure that < fx(x) < oo 
v-&.s., we require that Q[V^t = 0] = z^({0}) = or T > (2m)~ x . From (I2.87p . we can 
decompose {Vtr} as 

V tT = VttIit + V^PM (° < * < T )> ( 5 - 2 ) 

where (i) {(3(t)} is a Brownian bridge to the value at time 1, (ii) {^tr} is a gamma 
bridge to the value 1 at time T (with parameter m > 0), and (iii) given Vtt, is a 
GIG random variable with parameter set {mT — 1/2, |Vtt|, \/2m}. Also, from (I2.92p . 
we can decompose {V t r} as 

V t T = VrTltT + Y T fx (7^ - 7^) (0 < t < T), (5.3) 

where (i) /i = (m/2)~ 1//2 , (ii) {7^} an d {7iT} are independent gamma bridges to the 
value 1 at time T, and (iii) given Vtt, Yt is a random variable with density 

2mt 

V » h,H-v TT ^—^—(yVrr 4- >.-0^. (5 . 4) 

r[mTJ 2 /T (^tt) 

5.1.1 Example 

Recall that in (I2.75P we defined fc( m ,0,p) as 



5.2 Simulation 
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Define p > by 



then p satisfies p = fc( m ,e,p)- Now suppose that, ct pviori^ Vtt 

is a VG random variable 

with density f T m,e,p \x) as given by (I2.74p . We then have 

/t I 21 ) 

= e^(^ M )- 2mt /S M ^-0^ (5.7) 
for £ G [0,T). The transition law of {V^r} is 

= e ^-^ 2 (k (m ,e, P r 2mit ~ s) fl-Xy-x)<iy 

= ftf' P \y-x)dy, (5.8) 

and similarly 

Q [V TT E dy I V sT = x} = ft'^iy - x) dy. (5.9) 

Over the time period [0,T], {V^r} is then a two-parameter VG process (since p is a 
function of 9 and m). The scaled LRB {uVtr/p} is a three parameter VG process 
with parameters m, 6, and a. We have shown that an asymmetric VG process can be 
constructed as a random bridge of a symmetric VG process. 



5.2 Simulation 

In this section we will assume that we can generate sample paths for Brownian and 
gamma bridges (see Avramidis et al. Q|, and Ribeiro & Webber 79]), gamma and GIG 
random variates (see Devroye 26j), and random variates from the law of Vtt- 

From the constructions (15. 2p and (15. 3p . we have two natural ways of simulating 
paths of VGRBs. Indeed, all we need to do is be able to sample from the distributions 
of Yjp and Ey, given the value of Vtt- Since the conditional distribution of is GIG, 
this causes no problems. To sample from the conditional distribution of Y T , we note 
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the following: If Vtt < then, from ( 15. 4| . the density of Yt is proportional to 

t {y> -v TT} (y(VTT + y)) mT - l e- 2my < t {y> -v TT} y 2mT - 2 e- 2m y 

txl {y>-v T T}9j n _j_(y), (5-10) 

2m 

where gl (x) is the gamma density given in (I2.46p . If Vtt > then the density of Yt 
is proportional to 

MyyoMVTT + y))^- 1 ^^ < My>-V TT] {VTT + y? mT - 2 Z- 2m{VTT+y) 

<xt {y>0} g^(V T T + y). (5.11) 

2m 

So, given Vtt, if T > (2m) _1 we can simulate Yt with the following acceptance-rejection 

algorithm: 

If V TT < 

1. Generate G from the gamma distribution with density g^ 1 x (Vtt + J/)- 

_ 2m 

2. If G < —Vtt go to step 1. 

3. Generate U from the standard uniform distribution. 

4. If (1 + V TT /G) mT - 1 < U go to step 1. 

5. Return G. 
If Vtt > 

1. Generate G from the gamma distribution with density g^j_{VTT + y)- 

2m 

2. If G < Vtt go to step 1. 

3. Generate U from the standard uniform distribution. 

4. If (1 + V TT /G) 1 - mT < U go to step 1. 

5. Return G — Vtt- 



5.3 VG equity model 
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Madan et al. [67| used an asymmetric VG process to model the log-returns of a single 
stock, and priced European options on the stock. It follows from the asymptotic 
behaviour of the Bessel function K u [z] given in (I2.70p that the exponential moment of 
the VG distribution exists only if 9 + a 2 /2 < m, which we now assume is true. We 
model the stock price as 



S t = S exp [rt + L t + wt] (t > 0) , (5.12) 



where, under the risk-neutral measure, {L t } is a three-parameter VG process with 
parameters (m, 9, a), r > is the constant rate of interest, and 



w = m log 



9 a 



2 

( > 

1 

m 2m 



(5.13) 



The characteristic function of {L t } was given in (12.621) . from which it becomes apparent 
that the drift term wt ensures that the discounted stock price process {e~ rt S t } is a 
martingale under the risk-neutral measure. 

We will show that this asymmetric VG model can be recovered using the information- 
based approach when the information process is a symmetric VGRB. First we fix some 
suitably distant future date T (for example, this could be the maturity date of the 
longest-dated, liquidly-traded, European option on the stock). We then let X T be a 
VG random variable with density f^ n ' 9 ' p \x) } where p is given by (15. 6p . and set 

h(x) = S exp(rT + ax/p + wT). (5-14) 

Let the information process {£tr} be a standard VGRB with parameter m such that 
£tt = Xt- From Section I5TT1 we have that {cr^tr/p} is a VG process with parameters 
{m, 9, a}. We then have 

Err = exp(-r(T - t))E[h(X T ) | £tr] 

= S exp(rt + oi tT jp + wT) E[e^ TT ^ tT)/p ] 

= S exp(rt + cr£ tT /p + wt), (5.15) 

for t G [0, T]. Hence we have {H tT } <t<T l = {St}o<t<r- European options on the price 



process {H tT } can then be calculated using the techniques in [67]. 
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5.4 VG binary bond 



5.4 VG binary bond 



In the example of Section 14.41 we derived the price of a binary bond in the Brownian 
information model. Brody et al. Vj\ derived a similar price process, only they included 
a rate parameter a > 0. This rate parameter behaved like the volatility parameter 
in the Black-Scholes model. We can construct a similar model using a VGRB as the 
information process. 

Let Xt be an X-factor such that Q[Xt = 0] = p and Q[Xt — a] — 1 — p. Set 
h(x) — x/a, and let Ht = K{Xt) be the redemption amount of a credit-risky, zero- 
coupon bond. We assume that the market filtration is generated by a standard VGRB 
{fa} such that £tt = X T . For {fa} to be well defined, its parameter m must satisfy 
T > (2m) . The time-t price of the bond is then given by 



B tT = P tT E[h(X T ) | fa] = Ptr Q[X T = a \ fa]. 
From Section l4~4l we have 



(5.16) 



B 



tT — r tT 



PtT I 1 + C 



fa 


m(T-t)-± X m (T-t)-l/2 






o" - fa 


Km[T-t)-l/2 


y/2m{a-faf 



(5.17) 



where 



c = 2p 



mT 1 7>- 

mo -2\ 2 4 JVmT-1/2 



V2 



2 y (1 -p)r[mT- 1/2] 
See Figures 15.11 15.21 and 15.31 for example simulations. 



(5.1* 



5.4 VG binary bond 
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(c) Information process when m = 25. (d) Bond price with m = 25. 




(c) Information process when m = 100. (f) Bond price with m — 100. 

Figure 5.1: Simulations of VG information processes and bond price processes in cases 
where the bond defaults. Various values of the parameter m used. The other parameter 
values are fixed as r t = 0, p = 0.3, T — 1, and a — 1. 
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5.4 VG binary bond 




(c) Information process when m = 100. (f) Bond price with m — 100. 

Figure 5.2: Simulations of VG information processes and bond price processes in cases 
where the bond does not default. Various values of the parameter m used. The other 
parameter values are fixed as r t = 0, p = 0.3, T — 1, and a — 1. 
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Figure 5.3: Simulations of VG information processes and bond price processes in cases 
where the bond defaults. Various values of the information rate parameter a are used. 
The other parameter values are fixed asr t = 0,p = 0.3,T = l, and m = 25. 



Chapter 6 

Stable- 1/2 information 



The term 'stable process' refers here to a strictly stable process with index a G (0,2); 
thus, we are excluding the case of Brownian motion (a = 2). The use of stable processes 
for the modelling of prices in financial markets was proposed by Mandelbrot [69( in 
connection with his analysis of cotton futures. In the tails, the Levy densities of stable 
processes exhibit power-law decay As a result, the behaviour of stable processes is 
wild, and their trajectories exhibit frequent large jumps. The variance of a stable 
random variable is infinite. If a < 1, the expectation either does not exist or is infinite. 
This heavy-tailed behaviour makes stable processes ill-suited to some applications in 
finance, such as forecasting and option pricing. 

To overcome some of the drawbacks of the stable processes, so-called tempered 
stable processes have been introduced (see Cont & Tankov [24J], for example, for details). 
A tempered stable process is a pure-jump Levy process, and its Levy density is the 
exponentially dampened Levy density of a stable process. The exponential dampening 
of the Levy density improves the integrability of the process to the extent that all the 
moments of a tempered stable process exist. Tempered stable processes do not possess 
the time-scaling property of stable processes. 

In this chapter we apply stable-1/2 random bridges to the modelling of cumulative 
losses. The techniques presented can equally be applied to cumulative gains. The 
integrability of a stable-1/2 random bridge depends on the integrability of its terminal 
distribution. At some fixed future time, the nth moment of the process is finite if 
and only if the nth moment of its terminal value is finite. Thus a stable-1/2 random 
bridge with an integrable terminal distribution can be considered to be a dampened 
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stable-1/2 subordinator. In fact, the stable-1/2 random bridge is a generalisation of 
the tempered stable-1/2 subordinator. If the Levy density of a stable-1/2 subordinator 
is exponentially dampened, the resulting process is an IG process. We shall see that 
the IG process is a special case of a stable-1/2 random bridge. 

We look in detail at the non-life reserving problem. An insurance company will 
incur losses when certain events occur. An event may be, for example, a period of high 
wind, a river flooding, or a motor accident. The losses are the costs associated with 
recompensing policy-holders who have been disadvantaged by an event. These costs 
might, for example, cover repairs to property, replacement of damaged items, loss of 
business, medical care, and so on. Although a loss is incurred by the insurance company 
on the date of an event (the 'loss date' or 'accident date'), payment is generally not 
made immediately. Delays will occur because loss is not always immediately reported to 
the company, the full extent of the costs takes time to emerge, the insurance company's 
obligation to pay takes time to establish, and so forth. 

In return for covering policy-holder risk, the insurance company receives premiums. 
The premiums received over a fixed period of time should, typically, be sufficient to 
cover the losses the company incurs over that period. Since losses can take years to pay 
in full, the company sets aside some of the premiums to cover future payments; these are 
called 'reserves'. If the reserves are set too low, the company may struggle to cover its 
liabilities, leading to insolvency. Large upward moves in the reserves due to a worsening 
in the expected future development of liabilities can cause similar problems. If the 
reserves are set too high, the company may be accused by shareholders or regulators of 
inappropriately withholding profits. Hence, it is the interest of the company to forecast 
its ultimate liability as accurately as possible when deciding the level of reserves to set. 

We use a stable-1/2 random bridge to model the paid-claims process (i.e. cumu- 
lative amount paid to-date) of an insurance company. The losses contributing to the 
paid-claims process are assumed to have occurred in a fixed interval of time. Sometimes 
claims-handling information about individual losses is known, such as that contained 
in police or loss-adjuster reports. In the model presented here, such information is 
disregarded, and the paid-claims process is regarded as providing all relevant infor- 
mation. We derive the conditional distribution of the company's total liability given 
the paid-claims process. We then estimate recoveries from reinsurance treaties on the 
total liability. The expressions arising in such estimates are similar to the expectations 
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6.1 Stable-1/2 random bridge 



encountered in the pricing of call spreads on stock prices. 

We shall examine the upper-tail of the conditional distribution of the ultimate 
liability, and find that it is as heavy as the a priori tail. This has an interesting 
interpretation in the case when the insurer is exposed to a catastrophic loss. At time 
t < T, the probability of a catastrophic loss occurring in the interval [t, T] decreases as t 
approaches T. However, in some sense, the size of a catastrophic loss does not decrease 
as t approaches T, since the tail of the conditional distribution of the cumulative loss 
does not thin. 

When the a priori total loss distribution is a generalized inverse-Gaussian distri- 
bution, we find that the model is particularly tractable. We present a family of special 
cases where the expected total loss can be expressed as a rational function of the current 
value of the paid-claims process. That is, each member of the family is a martingale 
that can be written as a rational function of an increasing process. 

The model can be extended to include more than one paid-claims process. We 
consider the case where there are two processes that are not independent, and which 
have different activity parameters. We then have two ultimate losses to estimate. We 
provide expressions for the expected values of the ultimate losses given both paid- 
claims processes. The numerical computations required to evaluate these expressions 
are no more difficult than those of the one-dimensional case. We demonstrate how to 
calculate the a priori correlation between the ultimate liabilities. The correlation can 
be used as a calibration tool when modelling cumulative losses arising from related 
lines of business (e.g. personal motor and commercial motor business). 

We also describe how to simulate sample paths of the stable-1/2 random bridge, 
and how to use a deterministic time-change to adjust the model when the paid-claims 
process is expected to develop non-linearly. 

6.1 Stable-1/2 random bridge 

In this chapter we take ft{x) to be the increment density of the stable-1/2 subordinator 
as given in (12.941) : that is 




(6.1) 



6.2 Insurance model 
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Let {£tr} be a process with law LRBc([0, T], {ft}, v)- We say that {£tr} is a stable-1/2 
random bridge. Recall that v must concentrate mass where fr{x) is positive and finite. 
Hence v must be a probability law with support on the positive half-line. 



6.2 Insurance model 



We approach the non-life insurance claims reserving problem by modelling a paid- 
claims process by a stable-1/2 random bridge. The stable-1/2 random bridge is a 
suitable candidate model for a paid-claims process because it is (i) increasing, and (ii) 
tractable (in particular, its density has a simple form). We shall look at the problem 
of calculating the reserves required to cover the losses arising fro m a singl e line of 



Norbers 



gl HQ 



business when we observe the paid-claims process. Arjas |2J and 
prov.de general descriptions of the problem. England & Verrall Q survey some of the 
existing actuarial models. Buhlmann |21[ and Mikosch |74| contain related topics. The 
present work ties in with that of Brody et al. [20] who use a gamma random bridge 
process to model a cumulative loss or gain. 

The method we use has a flavour of the Bornhuetter-Ferguson model from actuarial 



science 



161 ] (see also 34|). In implementing the Bornhuetter-Ferguson model, one 



begins with an a priori estimate for the ultimate loss (the total cumulative loss arising 
from the underwritten risks). Periodically, this estimate is revised using a chain-ladder 
technique to take into account both the a priori estimate and the development of the 
total paid (or reported) claims to date. 

In the proposed model, we assume an a priori distribution for the ultimate loss. 
By conditioning on the development of the paid-claims process, we revise the ultimate 
loss distribution using the Bayesian methods described in Chapter HI In this way, we 
continuously update the conditional distribution for the total loss. This is as opposed 
to the deterministic Bornhuetter-Ferguson model in which only a point estimate is 
updated. Knowledge of the conditional distribution allows one to calculate confidence 
intervals around the expected loss, and to calculate expected reinsurance recoveries. 

The main assumptions of the model are the following: 



1. The claims arising from the line of business have run-off at time T. That is, at 
time T all claims have been settled, and the ultimate loss Ut is known. 
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6.3 Estimating the ultimate loss 



2. Ut has a priori law v such that Ut > and E[C/|.] < oo. 

3. The paid-claims process {£*t} is a stable-1/2 random bridge, and £tt = Ut- 



4. The fresi estimate of the ultimate loss is U t r = E 
natural filtration of {£<t}- 



C/7 



7T? 



where {J-f} is the 



A few remarks should be made about assumption HI First, using the natural filtration 
of {£tr} as the reserving filtration means that the paid-claims process is the only source 
of information about the ultimate loss once the measure v is set. We do not consider the 
situation where we have access to information about claims that have been reported 
but not yet paid in full (such estimates). Second, the expectation is taken 

with respect to Q, which may or may not be the 'real-world' measure. Let us call Q 
the actuarial measure. When reserving, practitioners routinely discount data before 
modelling. Discounting may adjust the data for the time-value of money or for the 
effects of claims inflation. Claims inflation, and interest rates, though understood to be 
stochastic, usually only provide a small amount of uncertainty to the forecasting of the 
ultimate loss, relative to the uncertainty surrounding the frequency and (discounted) 
sizes of insurance claims. Furthermore, it is often for practical purposes reasonable 
to assume that claims inflation and interest rates are independent of claim frequency 
and size. Hence, a stochastic reserving model may lose little from the assumption 
that interest rates and inflation rates are deterministic. We make this assumption, 
and further assume that the paid-claims process has been discounted for the effects of 
interest and inflation. 



6.3 Estimating the ultimate loss 



The time-t conditional law of Ut is 



Vl(iZ) *(Ri&r) 

w)(^r^(-4(^-g)), ( d.-) 



(6.2) 



6.4 The paid-claims process 
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Given this law, confidence intervals and quantiles for the ultimate loss are readily 
calculated. The best-estimate ultimate loss is 



U 



IT 



zv t (dz). 



(6.3) 



At time t < T, the total amount of claims yet to be paid is Ut — £tT- The amount that 
the insurance company sets aside to cover this unknown amount is called the reserve. 
The expected value of the total future payments is called the best-estimate reserve, and 
can be expressed by 

Rtr = U tT - 6t- (6.4) 

For prudence, the reserve may be greater than the best-estimate reserve. However, for 
regulatory reasons it is sometimes required that the best-estimate reserve is reported. 
The variance of the total future payments is the variance of the ultimate loss, which is 



Var 



Ut — itT 



Var 



Un 



Ft 



{z-U tT fu t {dz). 



(6.5) 



6.4 The paid-claims process 

We give the first two conditional moments of the paid-claims process. From Corollary 
13. 8. 2} we have 

E [&r | 7f ] = + ^U sT . (6.6) 

L 1 J 1 — s 1 — s 

Equation ( 16 .6p implies that the paid-claims development is expected to be linear. We 
return to this point later. Using Proposition 12.8.51 and a straightforward conditioning 
argument, we have 

E [&} =fl°°^^-c(T-t)e^^<l> [-dTz-W] | u(dz) 



t 



= — E [U T ] - c(T - t)V2n / z 3/2 $ [-cTz" 1/2 ] u(dz). (6.7) 
T Jo 

Fix s < T and define the relocated process {r]rr}s<t<T by 



VtT — 6t — CsT- 



(6.8) 
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6.5 Reinsurance 



The dynamic consistency property implies that, given £ s t, {vtr} is a stable- 1/2 random 
bridge with the marginal law of t]tt being v*{A) = v s {A + £ s t)- Then we have 



2 

sT 



T- t 
T - s 



t - s 
T-s 



E[£#|f aT ] 



c(T-t)V2ir I (z 



3 <?{T-s) 2 



c(T 



i/ a (dz). (6.9) 



6.5 Reinsurance 



An insurance company may buy reinsurance to protect against adverse claim develop- 
ment. The stop-loss and aggregate excess-of-loss treaties are two types of reinsurance 
that cover some or all of the total amount of claims paid over a fixed threshold. Under 
a stop-loss treaty, the reinsurance covers all the losses above a prespecified level. If 
this level is K, then the reinsurance provider pays a total amount (Ut — K) + to the 
insurance company. The 'aggregate L excess of IC treaty is a capped stop-loss, and 
covers the layer [K, K + L] . In this case the reinsurance provider pays an amount 
(U T - K) + -(U T -K- L)+. 

The insurance company typically receives money from the reinsurance provider 
periodically. The amount received depends on the amount the company has paid on 
claims to-date. If the insurer has the paid-claims process {£,tr}, and receives payments 
from a stop-loss treaty (at level K) on the fixed dates t\ < t 2 < ■ ■ ■ < t n = T, then the 
amount received on date £j is 



(fc,T - K)+ - (fc_ lfT - K)+. (6.10) 

Recall the following expressions from Section [2781 the stable-1/2 bridge marginal den- 
sity function 

f t \ ft(v)fT-t(z-y) fR , 

ftriy, z) = — (6.11) 

JT{Z) 

( lc 2 (Ty-tz) 2 ' 

1 ct{T-t)^{—2^W 

— ^{0<y<z} 



v 7 ^ T (y- y 2 /zf 2 



6.5 Reinsurance 
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the stable-1/2 bridge marginal distribution function 



Ftr(y;z) = / f tT (x;z)dx 



c(Ty - tz) 
^yz(z - y) 



+ | 1 - ^ ] r lrl,[i - l] : •!• 



c((2t-T)y-tz) 
\fyz{z - y) 



(6.12) 



and the stable-1/2 bridge partial first moment 



M tT (y, z)= I z f tT (x; z) dx 
'o 



t 

T 



z<$ 



c(Ty - tz) 
^yz(z - y) 



_ e 2c 2 t(T-t)/z $ 



c((2t-T)y-tz) 
y/yz(z-y) 



(6.13) 



The expected value of reinsurance payments such as (I6.10p can be calculated using the 
following: 



Proposition 6.5.1. Fix t e (0, T). At time s < t, the expected exceedence of ^ t T over 
some fixed K > is 

D gt = E[(fa-K) + \j*] 



T-s 



T-s 



roc 

+ Mk>$ bT }(K - €st) / F t - a , T - a [K - fa] z - fa) v a (d2) 

J K 

poo 

-t {K> Z sT} / M t -s,T-s{K - &rl z - £st) v„(dz). 
Jk 



(6.14) 



Proof. If K < i S T then 



E [(fcr - K) + | J*] = E [fa I J*] - K 



T-s 

Thus we need only consider the case when K > C, s t- 
The T\- conditional law of fa is 



T-t t-s 

-fa + ^ U aT ~ K. 



T-s 



(6.15) 



(6.16) 
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6.5 Reinsurance 



Hence we have 



1 f°° 

Dst = nWT-\ / (y - y)ft-s(y - 6r) dy 

ip s {K; l; sT ) J K 





1 
1 








i 








l 






POO 






\\v 


JK 


JK 



( , /v - ) / fc-j z y) u{dz) ft _ s {y _ u) dy 

Jt{z) 



K 



(y - K) tt~\ d y u ( dz ) 

Jt{z) 



Making the change of variable x = y — £ s t yields 

poo pz-£ sT 

D st = / {x + £ sT - K)ft- s ,T-s{x; z - £ sT ) dx u a (dz) 

JK JK-£, sT 



(6.17) 



<x> 



+ -K) I {1 - F t _ s>T _ s (K - z - 6r)} v s (dz 
T-t . t-s 



K 



T-s 

"OO 



t - S f°° 

ZsT + ■=, U sT -K+ (K- £ sT )F t _ StT _ s (K - £ sT ; z - £ sT ) u s {dz) 

T-s J K 

POO 

- / M t _ S:T _ s (K-Z sT ;z-U)"s(dz). (6.18) 

JK 

□ 

Suppose that the insurance company has limited its liability by entering into a 
stop-loss reinsurance contract. At time s G [0, T), the expected reinsurance recovery 
between times t and u is 

E - K)+ - (&r - K) + \J*]= D su - D st , (6.19) 

for s < t < u < T. 

Using a similar method to the calculation of D st , we can calculate the expectation 
of conditional on it exceeding a threshold. For a threshold 6 > £ s t, we find 

l&l&T | t,sT, QtT > = 



1 ~ -G F t-s,T-s(8 - &r; ^ - £ sT ) z/ s (dz) 

(6.20) 

Sometimes called the conditional value-at-risk (CVaR), this expected value is a coherent 
risk measure, and is a useful tool for risk management (see McNeil et d. Q)- Note 



6.6 Tail behaviour 
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that CVaR is normally defined as an expected value conditional on a shortfall in profit. 
Since we are modelling loss, and not profit, the risk we most wish to manage is on the 
upside. Hence, conditioning on an exceedence is of greater interest. 



6.6 Tail behaviour 

In this section we consider how the probability of extreme events is affected by the paid- 
claims development. Suppose that the line of business we are modelling is exposed to 
rare but 'catastrophic' large loss events. In this case we assume that the a priori 
distribution of the ultimate loss has a heavy right-tail. If a catastrophic loss could 
hit the insurance company at any time before run-off, then it is important that any 
conditional distributions for the ultimate loss retain the heavy-tail property. We shall 
see that in the stable-1/2 random bridge model, the conditional distributions are as 
heavy tailed as the a priori distribution. 

Assume that Ut has a continuous density p(z) which is positive for all z above some 
threshold. Then the value of Ut is unbounded in the sense that 



h > x] > 0, for all x el. (6.21) 
Define 

Tail, = ton nH Q ^7 >L } lc , (6.22) 

If Taib = oo then the tail of the future-payments distribution at time t > is not as 
heavy as the a priori tail. That is, a catastrophic loss at time t is 'smaller' than a 
catastrophic loss at time 0. If Taih = then the tail of the future-payments distribution 
is greater at time t than a priori. If < Tail t < oo then the tail is as heavy at time t 
as a priori. Using L'Hopital's rule, we have 

Tail,= l lm 



« R ^i™Rirb (6 ' 23) 



for t G (0, T). Some examples follow: 
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6.7 Generalized inverse-Gaussian prior 



• If p(z) oc l{ 2 >o}e z (exponential) then Taih = iptiM', 6t) ^ tT . 

• If p{z) oc l{ z> o}e~ z2 (half- normal) then Tail f = ^(R; £tr) e^ T . 

• If oc l {2>oH - 3 / 2 e- 1 / 2 (Levy) then Tail, = ^(R;& T ). 

This property has an interesting parallel with the subexponential distributions. By 
definition, X has a subexponential distribution if 

i^o Q[X>L] =n ' (6 - 24) 



where {Aj}" =1 are independent copies of X (see Embrechts et al. 32]). We note that 



lim ^ r «azT>L) = 

l->oo Q [Z T - Z t > L I Z t ] 

for {^t} a Brownian motion, a geometric Brownian motion or a gamma process. If 
{Z t } is a stable-1/2 subordinator, so the increments of {Z t } are subexponential, then 

Q[Z T -Z,>L|Z t ] = 7^? (6 - 26) 

6.7 Generalized inverse- Gaussian prior 

The generalized inverse-Gaussian (GIG) distribution is a three-parameter family of 



distributions on the positive half-line (see J0rgensen |57| or Eberlein & v. Hammerstein 
29l | for further details). In the present work it has so far appeared only tangentially. 
We now provide some of its properties. The density of the GIG distribution is 

/ G/G (x;V,7) = l{x>o } (|) A 2^p] xA ~ leXP R(<^ _1 +7V)) • (6.27) 

Here K v [z\ is the modified Bessel function seen in Section l2~77l The permitted parameter 
values are 

5>0, 7>0, if A > 0, (6.28) 

<5>0, 7>0, if A = 0, (6.29) 

5 > 0, 7 > 0, if A < 0. (6.30) 

For A > 0, taking the limit 5 — > + yields the gamma distribution. For A < 0, 
taking the limit 7 — > + yields the reciprocal-gamma distribution — this includes the 
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Levy distribution for A = —1/2 (recall that the Levy distribution is the increment 
distribution of stable- 1/2 subordinators) . The case A = —1/2 and 7 > corresponds 
to the inverse- Gaussian (IG) distribution. If X has density (I6.27P then the moment 
fjLk = E[X fc ] is given by 

^ = ^fffi [ -) for AGR, 5>0, 7 >0, (6.31) 



Kx[jS) V7 

r[A + k] ( 2 x '■ 



/;< \ r[A] V7V A and A > 0, 5 = 0, 7 > 0, (6.32) 
00 k < —A 

T[-X-k] (6 2 \ k k< _ x 

fi k = <J r[-A] V 2 / " and A < 0, 5 > 0, 7 = 0. (6.33) 
00 k > —A 



For convenience, we recall some facts about the IG process that appeared in Section 
12.10.11 The IG process is a Levy process with increment density 

ct 1 / 7 2 / \ 2 \ 

q t (x) = l {g > } ^ xV2 ex P { x -p) )■ ( 6 - 34 ) 

We see that qt(x) = fGio(x', —\ 1 ct,^). The kth. moment of qt(x) is 

m[ k) = \fl^ ct (- K k . 1/2 [yct], (6.35) 



7T \ 7 

for fc > 0. 



6.7.1 GIG terminal distribution 

We shall see that the GIG distributions constitute a natural class of a priori distribu- 
tions for the ultimate loss. With 7 > and c > fixed, we examine some properties of 
a paid-claims process {£tr} with time-T density faic^z] \ cT, 7). The transition law is 

tT e dy I i sT = x] = \ ft-s(y ~ x) dy, (6.36) 

ip s {K; x) 

TT€iyUlT = x] = t^, (6.37) 
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6.7 Generalized inverse-Gaussian prior 



where 



^o(d*; = faia(z; A, cT, 7 ) dz, (6.38) 



^(dz;0 = (1-£)1 { *> C} (l-e/*) 3 / 2 W^;A,cr, 7 )dz. (6.39) 



Writing 



we have 



*=uo iK^m (6 - 40) 



^ t (R;y) = «(l-| ! )e-^''y , A +^ dz 

= «(1- «)e-3^y (z + y) A +2- dz 

= ^e-^«^ T -*) jH(* + y) A+ *«r_ t (*) dz. (6.41) 
Given £ t r = 1/, the best-estimate ultimate loss is given by 

U^Wl ^ V) = ^-^AJ- . ,6.42) 

6.7.2 Case A = -1/2 
When A = —1/2 we have 

MR;y) f f n t 1 g(* - g ) f 7 2 ((y-^)-c(t-s)/7)' 

-ft-s(y-X) = 1 {y-x>0}-7= 7 w^ ex P 



= ft- 8 (y-x). (6.43) 
Hence {£tr} is an IG process. Note that in this case {£tr} has independent increments. 

6.7.3 Case \ = n-\ 

We now consider the case where A = n — |, for n e N+. For convenience we write 

gfto = /gjg(x; 1/2, ct, 7 ). (6.44) 
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Hence we have ql°\x) = qt(x). The transition density of {£<t} is then 

ELo(>fcV 



g«- a (y-g) ^ fe=UVfc/ ~ " • (6-45) 



When n = 1 this is 



M&£ f (v - x) -a (v x ) V+ " {T ~ t] 
i^x) h - s[V X} - qt - s{y X) x + ±(T-s) 

-(i-^^4jy)A(,-) 

+ (6.46) 

Thus the increment density is a weighted sum of GIG densities. We shall now derive 
a weighted sum representation for general n. We can write 

/ (z + y) n q T -t(z) dz = / ((* + x) + (y - x)) n q T _ t (z) dz 
Jo Jo 

{z + x) k q T -t(z) dz 



k=0 x / j=0 

Then we have 



E (") (y - E (%f-7 } x\ (6.47) 

i._n V / '_n W / 



ELo a)(y-^)"- fe E; =0 



m T _ t x 



■3 



Efc=0 U™T- S ^ 



However, when G N , 

g^gt-£(g) = Jgig{z; -1/2, c(t-s), 7) 
/ / (z;fc-l/2, C (f- a ), 7 ) 

'c(t - s)\ fe X fc -i/2[7c(* - s)\ 



7 / %Wi-s)l 

(k) 



(6.49) 
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6.7 Generalized inverse-Gaussian prior 



Hence we have 



(y- x) n k q t - s (y-x) =m( l _ s k) qll s k \y -x). 
Using the identity (16.50p . (I6.48P can be expanded to obtain 



iptiM', y) 

ij s (R;x) 



ft- s {y - x) = ^2wi k \x) qf\{y - x ) 



where 



wit\x) 



(3 



m 



k=0 



(n-k) sr^k fk\^(k-j) j 



t-s 



Ek /k\ 
7=0 \ j) 



m T-t X 



(6.50) 



(6.51) 



(6.52) 



E;=o C)'"r 

Notice that w^\x) is a rational function whose denominator is a polynomial of order n, 
and whose numerator is a polynomial of order k < n. Thus the transition probabilities 
°f {£tr} depend on the first n integer powers of the current value. The conditional law 
of the ultimate loss is 



11 Ck ™( n - fc ) 



dy. 



(6.53) 



sT m T-s 



We can verify that Y^k=o w st\ x ) = 1 using the fact that IG densities are closed 
under convolution. We have 



q T - s {z) = / q T -t{y)qt-s{z - y) dy, for < s < t < T. 
Jo 

For fixed n £ N+, we then have 

n 

E 



(6.54) 



fc=0 



n \ ( n -k) k 

k W T _ s 'x 



(z + x) n q T - s (z) dz 

(z + x) n / q T - t {y)qt-s{.z -y)dydz 
Jo 

poo 

qr-tiy) j (z + x) n q t - s (z-y)dzdy 
■Jy 

POO 

qr-tiy) I (z + y + x) n q t - s (z)dzdy 



qt- s {y) 



k=0 



dy 



E 

k=0 
n 

E 

fc=0 



m 



(n-k) 



(y + x) k qt- s (y) dy 



m 



(n-k) 

t-s 2 < 

j=0 



(6.55) 
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which gives 



k=0 



(n—k) sr^k fk\ (k—j) 



(>r-7 fcj £= () 



TYlrrp . X 



j=0 \j)"*T-t 



- 1. (6.56) 



6.7.4 Moments of the paid-claims process 

The best-estimate ultimate loss simplifies to 



n+l (n+l\ (n+l-k) c k 



U tT = 



k=0 V k ) 



m 



T-t StT 



Era /ra\ 
k=0 \k) 



(n—k) £]~ 
m T-t MT 



For example, when n = 1 we obtain 

c(T - 0(1 + 7c(T - t)) + 2 7 2 c(T - t)6r + 7 3 ^t 



17, 



7 2 c(T _ t)+7 3 6T 



By similar calculations, we have 



n+m (n+m\ {n+m—k) j-k 



E[[/-|&r] 



En+m /ra+m\ 
fc=0 I fc J 



T-i 



tr 



Era /n\ 
fc=0 VaJ 



for m G N + , 



rn 



T-t MT 



and 



E 



tr 



Z^fc=0 \k)" L T-t StT 



tr 



(6.57) 



(6.58) 



(6.59) 



exp(|a 2 6T-(T-t)(7-7)), (6-60) 



for < a < 7, where 7 = a/7 2 — a 2 , and is the fcth moment of the IG distribution 
with parameters 5 — ct and 7 = 7. 



6.8 Exposure adjustment 

We have seen that 

Efor] = |e[C/ t ]; (6.61) 

thus in the model so far the development of the paid-claims process is expected to be 
linear. This is not always the case in practice. In some cases the marginal exposure is 
(strictly) decreasing as the development approaches run-off. This manifests itself as 

j^E[6t] < 0, (6.62) 
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6.8 Exposure adjustment 



for t close to T. A straightforward method to adjust the development pattern is through 
a time change. We describe the marginal exposure of the insurer through time by a 
deterministic function e : [0, T] — > 1R + . The total exposure of the insurer is 

cT 

e(s)ds. (6.63) 



We define the increasing function r(t) by 

r(t) = T Jo T \ [ . (6.64) 
Jo ^{s)ds 

By construction r(0) = and r(T) = T. 

Now let r(t) determine the operational time in the model. We define the time- 
changed paid-claims process {<Q~ T } by 

& = e(r(*),r), (6.65) 

and set the reserving filtration to be natural filtration of {^}. Then we have 

Jo £ ( s ) ds 
J T e(s) ds 

and 



E &t] = ;°t , ; ; E[t7 r ] (6.66) 



Se[&] = 7#7TT^)- (6 - 67) 

Jo e ( s ) ds 

6.8.1 Craighead curve 

Craighead {25] proposed fitting a Weibull distribution function to the development 
oattern of paid claims for forecasting the ultimate loss (see also Benjamin & Eagles 
1^|). In actuarial work, the Weibull distribution function is sometimes referred to 
as the 'Craighead curve'. To achieve a similar development pattern we can use the 
Weibull density as the marginal exposure: 

e(t) = -(x/a) b - 1 e ix/a) \ (6.68) 
a 

for a, b > 0. Then the time change r(£) is the renormalised, truncated Weibull distri- 
bution function 

1 _ e (t/*) b 



6.9 Simulation 
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See Figure I67T1 for plots of this function. When b < 1, r'(t) is decreasing. Under such 
a time change, the marginal exposure is decreasing for all t G [0, T\. When b > 1, r'(£) 
achieves its maximum at 

fb-l\ l/b 

r = a VJ ' ( 6 - 7 °) 

and r'(t) is decreasing for t > t*. Thus if T > t* then the marginal exposure is 
decreasing for i 6 [f,T]. If T < £* then the marginal exposure is increasing for 

te[o,T]. 




0.2 0.4 0.6 0.8 1.0 



Figure 6.1: Plots of the truncated Weibull time change for various parameters, and 
with T — 1. The expected paid-claims development of the model will have the same 
profile as r(t) (scaled by E[Z7t]). Hence, under any one of the above time changes, 
when t is close to T the marginal exposure falls (i.e. J^E^^] < 0). 

6.9 Simulation 

We consider the simulation of sample paths of a stable-1/2 random bridge. First, we 
can generalise (12.1031) to 

y + -( z -y)ll + — Z =), (6.71) 

where < s < t < T, and Z is a standard Gaussian. We can then generate a discretised 
path {£(ti, T)}?" , where = iT2~ n , by the following recursive algorithm: 
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6.9 Simulation 



1. Generate the variate £(T, T) with law u, and set £(0, T) = 0. 

2. Generate £(f ,T) from £(0,T) and £(T,T) using identity (IBTTTj) . 

3. Generate £(j,T) from £(0,T) and £(~,T), and then generate £(^r, T) from 
|(|, T) and £(T, T). 

4. Generate £(f ,T), £(f ,T), £(f ,T), £(f ,T). 

5. Then iterate. 




(c) c = 7 (d) c = 10 



Figure 6.2: Simulations of the paid-claims process {£tr} (bottom line) and the best- 
estimate process {^t} (top line). Various values of the activity parameter c are used. 
A priori, the ultimate loss Ut has a generalized Pareto distribution (GPD) with density 
fcpoix) = l{x>i} (l + ^ip) 5 • (This is the GPD with scale parameter a = 1, location 
parameter /x = 1, and shape parameter £ = 1/4.) 



6.10 Multiple lines of business 
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6.10 Multiple lines of business 



We shall generalise the paid-claims model to achieve two goals. The first is to allow 
more than one paid-claims process, and allow dependence between the processes. The 
second is to keep the dimensionality of the calculations low with a view to practicality. 
The following results can be applied to the modelling of multiple lines of business or 
multiple origin years when there is dependence between loss processes. 

6.10.1 Two paid-claims processes 

We consider a case with two paid-claims processes, but the results can be extended to 
higher dimensions. In what follows, we set f£(x) = ft{x) as given by ( 16. ip . and ff T (x) = 
ftr{x) as given by (16.12p . Here we have introduced the superscript to emphasise the 
dependence on c. Let {S(t, T*)} be a stable- 1/2 random bridge with terminal density 
p(z) = v(dz)/dz, and with activity parameter c. Fix a time T < T*, then define two 
paid-claims processes by 




S(t, T*) 

k 2 S(Xt + T,T*) -k 2 S(T,T*) 



(0 < t < T), 
(0 < t < T), 



(6.72) 
(6.73) 



where A = T*/T — 1, and k = c 2 /(cA) for some c 2 > 0. The density of £, 



is given by 




(6.74) 



and the density of £, 



is 




(6.75) 



(6.76) 



(6.77) 
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Here (I6.75P follows after a change of variable, (16.76P follows from the definition of 
ftr{y\z) given in (12.2ip . and (16.771) follows from the functional form of ftT(y,z) for 
stable-1/2 bridges given in (l2"T9"9"]) . It follows from Corollary KE1\ that is a 

stable-1/2 random bridge with law LRBc([0, T], {f^},p^ 1 '(z)dz), and (combined with 
the scaling property of stable-1/2 bridges) has law LRBc{[0,T], {f£ 2 },p( 2 \z)dz). 

The conditional, joint density of (Q^ , k~ 2 ^) is 

$ e d Vl , k~ 2 $ e dy 2 1 = fc -2 £$ = ^ = 

/r.-(i+A)t( z - (2/1 + 2/2)) 1 

-p(^) dz S / t _ s (|/i - xi) d.2/1 f X (t-s){y2 ~ x 2 ) dy 2 , 



: =X!+x 2 /t*-(1+A)s( 2: l 3 ' 1 + X 2 



(6.78) 



for < s < i < T. Then we have 



SsT — x l) A SsT ~~ x 2 



/r*-(i+A)t( z 2/) 



2=xi+a;2 JT*-(1+A)s 



d.i/ 



/(i+A)(t-»),T*-(i+A) 8 ( z - (xi + x 2 );y- (x 1 + x 2 ))p(z)dz\dy; (6.79) 



and, given = xi and 2 C s t = x 2, the marginal density of is 



-2t(2) 



(1) 



2/1 ^ 



2 = X1+X2 



/i C _ SiT ._ ( i+A)s(2/i - 2 - (&i + x 2 ))p(z) dz, 



(6.80) 



and the marginal density of k 2 ^ is 



2/2 H> 



POO 

J z=x\ 



f. 



+X2 



\(t-s),T*-(l+\)s 



(2/2 - ar 2 ; 2 - (^1 + ^2)) p{z) dz. 



(6.81) 



6.10.2 Correlation 

The a priori correlation between the terminal values is well defined when the second 
moment of v is finite. The correlation can be used as a tool in the calibration of the 
model. Assuming that K[S(T*,T*) 2 ] < 00, the correlation is defined as 



E 



t(2) 



— E 



E 



t(2) 



E 



E 



E 



t(2) 



(6.82) 



-E 



t(2) 
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We shall calculate each of the components of (16.82p separately. First, we have 



Noting that 



E 



E[S{T,T*)) = —E[S{T*,T*)). 



= k 2 (S(T*,T*) - S(T,T*)) 
l = k 2 S{T* -T,T*), 



(6.83) 



(6.84) 



we have 



E 



t(2) 



k 2 E[S(T* — T, T*)\ 

k 2 (l-^) E[S(T*,T*)]. 



The second moments of £^ and £jf^ follow from (16.71) . and are given by 



.(2) 



E 



(6 



(1) 
TT 



— E [S(T*,T*) 2 ] - (T* — T)Ct* 



and 



E 



fc ( 1 — 



T 



E [S(T*,T*) 2 ] - k 4 TC T * 



where 



Cn 



cV2tt 



,3/2 



e c -£-$[-cT*z- 1/2 } p{z)dz. 



(6.85) 



(6.86) 



(6.87) 



(6.{ 



The final term required for working out the correlation is the cross moment. This is 



E 4t4t = k 2 E[S(T,T*)(S(T*,T*) - S(T,T*))] 

= k 2 E [S(T, T*)S(T*, T*)} - k 2 E [S(T, T*) 2 ] . 

The first term on the right of (16.891) is 

fT( x )fr*-T(y ~ x ) 



(6.89) 



£; 2 E [S(T,T*)S(T*,T*)] = k 2 




xy 



'0 JO 

-oo POO 




xyfe T *{x;y) dxp(y) dy 



dxp(y) dy 



o Jo 



T 



k J y p(y) d v 

eL E [S(T*,T*n 



(6.90) 
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The second term on the right of (16.891) is given by (16.861) . Hence we have 



E 



k 2 (T*-T)C T *. 



(6.91) 



The expression for the correlation follows by putting (16 .83 p . (16. 85 p . f !6.86j> . (I6.87p . and 
(ESI]) together. 



6.10.3 Ultimate loss estimation 

We now estimate the terminal values of the paid-claims processes. At time t < T, the 



best-estimate ultimate loss of , } (or, indeed, }) depends on the two values £. 



(2), 



-(1) 
tT 



and The best-estimate ultimate loss of {^ T '} is 



t#>=E 



StT ~~ x l' StT ~~ x 2 



E [5(T, T*) | S(t, T*) = x u S(T + At, T*) - 5(T, T*) = £r 2 x 2 ] 

E [5(T + At, T*) | S(t, T*) = x x , S{T + At, T*) - 3{T, T*) = k- 2 x 2 ] - k' 2 x 2 

E [S(T + At, T*) | S(t, T*) = x x , S{{1 + A)t, T*) - S(t, T*) = k~ 2 x 2 ] - k~ 2 x 2 

(6.92) 

E [S(T + At, T*) | S((l + A)t, T*) =x x + k~ 2 x 2 ] - k~ 2 x 2 (6.93) 
(E [S(T*,T*) | 5((1 + A)t,T*) = X! + A;" 2 x 2 ] - AT 2 x 2 ) 



T-t 



T*-(1 + X)t 



T*-(T- t) 



;i + A)t 



(6.94) 



Equation (I6.92p holds since reordering the increments of an LRB yields an LRB with 
same law, (I6.93P follows from the Markov property of LRBs, and (I6.94p follows from 
(16. 6p . We also have 



r 

E [S{T*, T*) | S{{1 + A)t, T*) = Xl + k~ 2 x 2 ] = 

Jo 



zp t (z) dz, 



(6.95) 



where 



Pt{z) 



{z>Xl+k~ 2 X2} 



3/2 



z — (x x + k~ 2 x 2 ) 

c 2 ( (T* - (1 + A)t) 2 T* 2 



x exp — 



2 \z - (xx + k~ 2 x 2 ) 



p(z), (6.96) 
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for K a constant chosen to normalise the density Pt(z). Similarly, the best-estimate 
ultimate loss of {Qt} * s 



U§ = k 2 ^_^ + x) t ( E [S(T*,T*) | S((l + X)t, T*) = x 1 + k~*x 2 ] - Xl ) 

T — t 

+ -. r— x 2 . (6.97) 

T* - (1 + \)t K ' 

To compute both and U^) we need to perform at most two one-dimensional in- 
tegrals (the integral we need is f )6.95p . but we note that pt(x) includes a normalising 
constant K — which is found be evaluating a second integral). We are saved the com- 
plication of performing double integrals. 

To extend these results to higher dimensions we can split the 'master' process 
{Str} into more than two subprocesses. Regardless of the number of subprocesses 
(i.e. paid-claims processes), all of the best-estimate ultimate losses can be computed 
by performing at most two one-dimensional integrals. This makes such a multivariate 
model highly computationally efficient. 

It should be noted that the case where all the subprocesses are identical in law (in 
the example above, this is the case when k — 1), many of the results can be extended 
to an arbitrary master LRB using Corollary 13.8.41 Importantly, the computational 
efficiency can be achieved for an arbitrary LRB. 



Chapter 7 



Cauchy information 



In this chapter, as in the last, we examine the random bridges of a stable process. We 
look here at LRBs based on a driftless Cauchy process, which is strictly stable with 
index a — 1. We call these LRBs Cauchy random bridges (CRBs). The third moment 
of the (one-dimensional) marginal distribution of a CRB does not exist for t e (0, T), 
regardless of the choice of terminal distribution, since the corresponding moments of a 
Cauchy bridge do not exist. Integrability is obtainable for moments of order less than 
three. Hence, conditioning a Cauchy process on its time-T distribution can moderate 
its wild behaviour. This moderation cannot be as severe as exponentially-dampening 
the process's Levy density, which would produce a tempered-stable process. We pro- 
vide expressions for the first two moments of a CRB (when they exist). These follow 
from the analysis of Cauchy bridges in Section 12.91 We then outline an algorithm 
for the simulation of CRBs. The algorithm uses the representation of a Cauchy pro- 
cess as a Brownian motion time-changed by a stable-1/2 subordinator. Returning to 
information-based pricing, we consider the valuation of a binary bond when the infor- 
mation process is a CRB. We price a call option on the binary bond price from first 
principles. This proves tractable, and a closed-form expression is derived for the call 
price. Plots of simulations of bond price processes are provided at the end, which show 
the influence of the activity parameter c. 
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7.1 Cauchy random bridges 

As in (IZT2B . we define 

m = wtcvy (7 - 1} 

which is the density of a driftless Cauchy process. Let {Ztx} be a random bridge of 
a driftless Cauchy process such that the marginal law of Ztt is v - That is, {Z t rp} is 
a CRB with activity parameter c. Since Jt{x) is bounded and positive on the whole 
real line, Ztt can be restricted to take values in any Borel set of K. For fixed k > 0, 
let {Z t kT} be a CRB with activity parameter c such that Z^t^t l = k Ztt, i-e. 

Q[Z fcTjfcT < *] = Q[& Z fcT < z] = K(-°0> */*])■ ( 7 - 2 ) 

Then it follows from Proposition 12.4.11 that 

{/c _1 Z fet ,feT} < t < T = {Ztr} <t<T ■ ( 7 - 3 ) 
If E[Z TT ] < oo then we have 

E [Z tT | Z sT = x] = + ^-E[Z TT I = x], (7.4) 

and if E[Z TT ] < oo then, using f l27L30]) . 

E [Z 2 T | Z sT = x] = j^* 2 + ^Z^( E [ Z tt I = x] + c 2 (T - s)(T - *)), (7.5) 
for < s < t < T. From (17.41) and (17. 5p it is straightforward to show that 

Vai[Z tT | Z sT ] = 

{T Z! ){t ' S) (Z sT - E[Z TT I Z sT ]f + ^Var[Z TT | Z sT ] + c 2 (t - a)(T - t). (7.6) 
{1 — s) z 1 — s 

The ^^-conditional marginal law of Ztt is 

j) s (dz; Z sT ) 
M^Zst) 

z 2 +c 2 T 2 /j \ 



i/ s (dz) 



f°° z 2 +c 2 t 2 (d y 

Then the ^^-conditional marginal density and distribution functions of Z t T can be 



written 



/oo 
ft-s,T-a(y - Z s t; 2 - Z sT ) u s (dz), (7A 
-oo 
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7.2 Simulation 



and 

/CO 
F t -s,T-s{y ~ Z sT ; z - Z sT ) u s (dz), (7.9) 
-oo 

respectively. Here ftr(y; z) is given by (I2.124p . and Ftr{y \ z) is given in the statement 
of Proposition 12.9.11 

7.1.1 Example 



OO 

i=— oo" 



Consider the case where Ztt can only take values in the countable set {a^} 
Writing q{i) = z/({a,}), the transition law of {Z tT } is given by 

tT «= dy | z, T = „] = ^- ? )( f - s )E.^(^W ( _ 2 + c2(( _ <)V 10) 

TT = <Zj Z sT = V = r= 7T\ ( . v (7.11) 



and 



where 

«*» - { ^?/ + %_ tr (7 ' i2) 

7.2 Simulation 

The Cauchy process is a Brownian motion time-changed by an independent stable- 1/2 
subordinator. For {M^(t)} a Brownian motion and {St} a stable-1/2 subordinator, the 
joint density of the pair (St, W(St)) is 

„-z 2 /(2y) rp 

(y ,, W1{ „ >01 £__^e-^> (7 , 3) 
Then, conditional on W(St) = z, the density of St is 



i -(z 2 +c?T*)/(2y) 

/ !->■ l{y>0} 

Integrating (17. 14p yields 



V ^ Mv>o}^ 2 + c 2 T 2 )- . (7.14) 

z y 



T > y | W(S T ) =z] = exp ( - ^^^ ) for y > 0. (7.15) 



7.3 Cauchy binary bond 
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Given W(St), St is the reciprocal of an exponential random variable, and we have 



s T = 1 , (7.i6) 

log U 

where U is a standard uniform random variable. Conditioned on W(St), the time- 
change process {S t }o<t<T is stable-1/2 random bridge. 

An algorithm for the simulation of stable-1/2 random bridges appeared in Section 



z 2 + c 2 T 2 



16.91 If we can generate a value of Ztt from its a priori distribution then we can 
construct a CRB using the following algorithm: 

1. Generate a standard uniform random variate U, and then set St = ~(Z TT + 
c 2 T 2 )/{\ogU). 

2. Generate a path {S t } < t <T of a stable-1/2 random bridge to St- 

3. Generate a path {/3(i)}o<t<i of a Brownian bridge with the terminal value 0. 

4. Return {{S t /S T )Z T T + y/Sr~ (3(S t / S T )}o<t<T- 



Schehr & Majumdar |84j] describe a general algorithm for the simulation of a sym- 
metric stable bridge using Monte Carlo methods. This method could be adapted for 
the simulation of Cauchy random bridge paths. 



7.3 Cauchy binary bond 

We revisit the information-based binary bond model, this time taking the information 
process {£tr} to be a CRB. The X-factor Xt is the redemption amount of a credit-risky, 
zero-coupon, binary bond, and we set £tt — Xt- A priori, we have Q[Xt = k ] = p 
and Q[Xr = ki] — 1 — p, where k < k\. The ^-conditional distribution of £tt = Xt 
is given by 

orp tit n (^ , I 1 -p)(c 2 T 2 + kj)(c 2 (T - 1) 2 + (fc - £) 2 ) V 1 

For convenience, we define 

ft(0)=Q[&T = fc o |&r], (7-18) 
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and 

9t(l) = Q[6rr = h I & T ] = 1 - q t (0). (7.19) 

We define the function A(t, x) : R + x R -» R by 

A(t,x) = P tT E[X T \£ tT = x]. (7.20) 

Then the process {Pot A(£, £tr)}o<t<T is a martingale, and the bond price process {B t r} 
is given by 

B tT = A(t^ tT ). (7.21) 
After a straightforward calculation, we find 

where 

ao = pfc (c 2 T 2 + fc 2 )(c 2 (T - t) 2 + A; 2 ) + (1 - p)^(c 2 T 2 + A; 2 )(c 2 (T - t) 2 + A; 2 ), 

ttl = -2fc fci {p(c 2 T 2 + kl) + (1 -p)(c 2 T 2 + A; 2 )) , (7.23) 

a 2 = pk (c 2 T 2 + kl) + (1 - p)h(c 2 T 2 + A; 2 ), 

and 

O = p{c 2 T 2 + A; 2 )(c 2 (T - t) 2 + A; 2 ) + (1 - p){c 2 T 2 + A. 2 )(c 2 (T - t) 2 + A: 2 ), 

Pi = -2 (A^(c 2 T 2 + A; 2 ) + h(l - p)(c 2 T 2 + kl)) , (7.24) 

/3 2 =p(c 2 T 2 + A; 2 ) + (l-p)(c 2 T 2 + fc 2 ). 

7.3.1 Call option price 

We shall calculate the price of a call option on the bond price. We approach the 
problem from first principles, as opposed to referring to the general formula derived in 
Section B~3J 

At time s < t, the price of a call option on the bond price Btj< is 

C st = P st E[(B tT -K) + \U] 

= P st E[(A(t,£ tT )-K)+\£ sT ], (7.25) 
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where K is the strike price. Defining the set A t by 



A t = {x:A(t,x) >K}, 



(7.26) 



the call option price is 



P st E[t { s tTeAt} (A(t,&)-K)\£ sT ] 

P st E [t {itTeAt} A(t, 6t) I U] ~ PstK Q[i tT e A t | £ sT ]. 



(7.27) 



Note that when discussing the price of an option on an asset in Section I4.3[ we used 
the notation B t instead of A t . We have changed to A t here to avoid confusion with the 
bond price B t T- We proceed by finding an explicit expression for the set A t . 

For fixed t, classical calculus reveals that A(t, x) is maximised and minimised at 
the points 



respectively. As a function of x, A(t,x) is decreasing on the intervals {—oo,x\ and 
[x, oo), and is increasing on the interval [x,x]. If K > A(t, x) then C s t = 0, and if 
K < A(t,x) then C s t = B sT . See Figure I7TT1 for an example plot of x \- > A(t,x). 




(7.28) 



and 




(7.29) 




X 



X 



Figure 7.1: An example plot of A(t, function of x. 
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We define 



K*= lim h{t,x) = P tT ^. (7.30) 

x—^±oo p2 



When K G (x, x)\{K*}, the equation 

A(t,x)=K (7.31) 
is quadratic in x, and has the two real solutions 



2(P tT a 2 -K(3 2 ) ■ {L6Z) 

When K > K* we have x + < x~ , and when K < K* we have x + > x~ . When K = K* 
we set 

x+ = aiAs-QfrA) (? 33) 

x" = +oo; (7.34) 

in this way x + is the only solution of equation (17.31 j) . Then we can write 

(x + ,x~) iovK>K\ 
A t = r " (7-35) 

forf^fT. 

It follows from ( 173]) that 

l 

Q&r e A t 1 6r] = V<**} + Yl - > ( 7 - 36 ) 

i=0 

where 

= Fi_ s>T _ s (x ± - £ sT ; fci - U)- (7-37) 

We also have 

E [l teT6M A(*, £ fT ) I U] = [ A(t, y) Q[£«r G dy | £ sT ] 

= PtT / ^ ft-s,T- 8 (y - 6st; 2 - 6t) dy i/ s (dz). 

J 2:= — oo J 

(7.38) 
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Noting that 

/ fts.T-s (y-£sT',z- £ sT ) dy = 
J A t 

^{k<k*} + F t - 8 ,T-s{x~ - £ sT ; z - ^ sT ) - F t „ s>T _ s (x + - £ sT ; z - £ sT ), (7.39) 
equation (17.381) simplifies to 

E[l {itTeAt} A(t^ tT )\U] = 

i 

1 { k<k*}P*Hs,U) + P tT J2 k i1t(i) (F- t (i) - F+(z)) . (7.40) 

i=0 

The call price is then 

i 

C st = 1 { jc<jc. } (S«t - iVO + - P s t#) (Fr t (i) - . (7.41) 

i=0 

Simulations of binary bond prices in this model are given in Figures 17.21 and 17.31 
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1.0 r 

in 




(c) Information process when c = 10. (f) Bond price with c = 10. 

Figure 7.2: Simulations of Cauchy information processes and bond price processes in 
cases where the bond defaults. Various values of the activity parameter c are used. 
The other parameter values are fixed as r t = 0, p = 0.3, and T — 1. 
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Figure 7.3: Simulations of Cauchy information processes and bond price processes in 
cases where the bond does not default. Various values of the activity parameter c are 
used. The other parameter values are fixed as r t = 0, p = 0.3, and T — 1. 



Chapter 8 



Normal inverse- Gaussian 
information 



This chapter is a mirror of Chapter [5J only we now consider the NIG process instead 
of the VG process. Since the VG process and the NIG process are similar, we find that 
the NIG random bridge (NIGRB) is very similar to the VG random bridge. To avoid 
too much repetition, we briefly list the contents of this chapter: 

• It is shown that an asymmetric NIGRB is a scaled standard NIGRB. 

• An algorithm is provided for the simulation of NIGRB sample paths. 

• An exponential NIG equity model is shown to be a special case of a single-factor 
information-based model when the information process is an NIGRB. 

• The price of binary bond is derived in a single-factor information-based model. 

• Simulations of binary bond price processes are plotted. 



8.1 Normal inverse- Gaussian random bridge 



Let {l^r} be a standard NIGRB, so in the transition law (13.14p . we take ft{x) 
which we defined in (12.145)) to be 



ft (a \x) 



a 2 te ah 



a 



VaH 2 + 



fi a \x) 
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We assume that Ytt has marginal law v. Since < fj?'(x) < oo for all x G R, Ypr 
can be restricted to take values in any Borel set of 1R. 

8.1.1 Example 

Fix c > and 9 G M, and define 

fc 2 = c ~We 2 + c 2 , a 2 = cV# 2 + c 2 . (8.2) 
Suppose that {^*t} is a standard NIGRB with parameter a and terminal density 



Equation ( 18. 3 p is the time-T density of an NIG process with parameters c, 9, and k. 
From (12.1461) and (I2.147p . we can write 



We then have 



/ (a t\ fr-t{ z Ac,8,k) f \ j 
W^O = — „ M , N fr ( z ) dz 



= fMf)( z _ ^ dZj (8 .5) 

for t G [0, T). For < s < t < T, the transition law of {5^t} is 

Q[Y tT e d y\Y s ^x^^fn(y-^y 

= f^ M {y-x)d V . (8.6) 

Similarly 

Q [Y TT G dy | Y sT = x} = ft?\y ~ x) dy. (8.7) 

Over the time period [0,T], {^r} is a two-parameter NIG process (since A; is a func- 
tion of 9 and c). The scaled LRB {aY tT /k} is a three parameter NIG process with 
parameters c, 9, and a. 
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8.2 Simulation 

Recall that a standard NIG process can be written as {W(X t )}, where {W t } is a 
standard Brownian motion, and {X t } is an independent IG process. The joint density 
of (X T ,W(X T )) is 



cTe c2T 1 / „ 2 _, 1 ,\ ex P - 



i-» t {y> o} C 4= 3/2 exp (- c ^(T 2 y 1 + y) 



2y 



Given = z, up to a normalisation the density of Xt is 

/ i f c 2 T 2 + z 2 \ \ 
y i y t {y>0} y- 2 exp -- + c 2 y )). (8.9) 



2 V 2/ 

Hence the conditional distribution of X? is generalized inverse-Gaussian with density 
fcwiy] -1, Vc 2 T 2 + z 2 , c), where 

/ G/G (x; A, 5, 7 ) = l { ,>o } 2K x ^ xX ~ lex P ("K^" 1 + 7^)) • (8-10) 

Furthermore, given VF(Xr) the time-change {X t }o<t<T is a stable-1/2 random bridge 
(since an IG random bridge is a stable-1/2 random bridge). We can simulate a sample 
path of {5^t} by the following algorithm: 

1. Generate a variate Ytt from the probability law v. 

2. Generate a variate X?t from the GIG distribution with parameters A = — 1, 



5 = ^/c 2 T 2 + Y 2 T , and 7 = c. 

3. Simulate a path {X^} of a stable-1/2 random bridge with parameter c, and 
terminating at value X?t at time T. 

4. Simulate a standard Brownian bridge {P{t)}o<t<i, terminating at value at time 
1. 



5. Return {X^Ytt I 'Xtt + yXTT(3(X t T/XTT)}o<t<T- 



We note that Ribeiro & Webber [80( apply related bridge-based constructions to the 
simulation of the NIG process. 



8.3 NIG equity model 
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8.3 NIG equity model 

In a similar way to the VG equity model of Section 15.31 we can use a three-parameter 



NIG process to model the log-returns of a single stock (see Ribeiro & Webber [80]). It 
follows from the asymptotic behaviour of the Bessel function K v [z\ (see (12.701) ) that 
the exponential moment of the NIG distribution exists only if c 2 > cr 2 + 29, which we 
assume holds. We model the stock price as 

S t = S exp [rt + L t + wt] {t > 0), (8.11) 

where, under the risk-neutral measure, {L t } is an NIG process with parameter set 
{c, 9, cr}, r > is the constant rate of interest, and 



w = c 



Vc 2 -a 2 -29- c 2 . (8.12) 



Since E[e Lt ] = exp(— wt), the drift term wt ensures that the discounted stock price 
process {e~ rt S t } is a martingale under the risk-neutral measure. 

We now show how to recover this model using the information-based approach when 
the information process is a symmetric NIGRB. We let Xt be a NIG random variable 
with density ftf' ' k \x), where k = k(c,9) is given by (I8.2p . and set 

h{x) = exp(rT + ax/k + wT). (8.13) 

Let the information process {£tr} be a standard NIGRB with parameter a given by 
(18. 2p . such that £tt = Xt- Then {uC,tT/k} is an NIG process with parameter set 
{c, 9, a}. Furthermore 

H tT = exp(-r(T - t))E[h(X T ) | &] 

= exp(rt + aitr/k + wT) E[e a(?TT - ?tT)/fc ] 

= exp(rt + a£ tT /k + wt), (8.14) 

for t G [0,T]. Hence we have {H tT } < t < T = {S t } <t<T- 



8.4 NIG binary bond 

We price a binary bond in an NIG information model including a rate parameter a > 0. 
We let X T be an X-factor such that Q[X T = 0] = p and Q[X T = a] = 1 - p. Setting 
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h(x) = x/a, we suppose that Ht = H{Xt) is the redemption amount of a credit-risky, 
zero-coupon bond. We assume that the market filtration is generated by a standard 
NIGRB {itr} such that £tt — X T . The time-t price of the bond is then given by 



B tT = P tT E[h(X T )\£ tT } 
= P tT Q[X T = a\£ tT } 



From Section I4.4[ we have 



-1 



B 



tr 



PtT\l+ 



P 



tr 



fi a) (P)f!&{<r-ZiT) 1 -P, 
1 + 71 



'aPiT-ty + ia-fa) 2 R i 






a*(T-ty + e tT Kl 


a^a 2 {T-t) 2 + {<j-£ tT ) 2 



where 



7 



P 



a 2 T 2 K x [ay/a 2 T 2 + <r 5 



1 - p V a 2 T 2 + a 2 K x [a 2 T] 
See Figures 18. l\ 18. 2^ and 18.31 for example simulations. 



(8.15) 
(8.16) 



.17) 
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(e) Information process when c = 100. (f) Bond price with c = 100. 



Figure 8.1: Simulations of NIG information processes and bond price processes in cases 
where the bond defaults. Various values of the parameter c used. The other parameter 
values are fixed as r t = 0, p = 0.3, T — 1, and a = 1. 
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Figure 8.2: Simulations of NIG information processes and bond price processes in cases 
where the bond does not default. Various values of the parameter c used. The other 
parameter values are fixed as r t = 0, p = 0.3, T — 1, and a — 1. 
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0.0 0.2 0.4 0.6 0.8 1.0 



(c) Information process when a = 10. (f) Bond price with a = 10. 

Figure 8.3: Simulations of NIG information processes and bond price processes in cases 
where the bond defaults. Various values of the information rate parameter a are used. 
The other parameter values are fixed asr i = 0,p = 0.3,T = l, and c = 10. 



Chapter 9 



Poisson information 

In this chapter we present some properties of a discrete random bridge, namely, the 
Poisson random bridge (PRB). PRBs are counting processes, and we shall see that 
they are Poisson processes with a state-dependent intensity. 

For a PRB, the distribution function of the waiting time to the next jump takes 
an interesting form. It can be written as an elementary transformation of the proba- 
bility generating function of the PRB's final value. By differentiating this distribution 
function, we find that that the intensity of the PRB is a linear function of the PRB's 
expected final value. 

We consider the case when the terminal value of a PRB has a negative binomial 
distribution. We find that all of the increments of this PRB have a negative binomial 
distribution. The intensity of this process is a linear function of the current value. 
When a jump occurs, this causes the intensity to jump, and in the periods between 
jumps the intensity decreases smoothly. We see therefore that this process exhibits a 
kind of contagion in the sense that every jump instantaneously increases the probability 
of another jump happening in some fixed future time interval. A PRB whose terminal 
value has a log-series distribution also has negative binomial transition probabilities. 
It can be considered to be a limiting case of a PRB with a negative binomial terminal 
distribution. 

A mixed Poisson distribution is a Poisson distribution with a mixed mean. Thus, 
the mixed Poisson distribution is a distribution on the non-negative integers. When 
the terminal distribution of a PRB is mixed Poisson, then every increment of the PRB 
has a mixed Poisson distribution. Such a process is called a mixed Poisson process. If 
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the gamma distribution is used to mix the Poisson mean, the negative binomial PRB 
is recovered. Since the log-series distribution is a distribution on the positive integers, 
it is not a mixed Poisson distribution. In fact, if the terminal distribution of a PRB 
takes values in a strict subset of Mo, then it is not a mixed Poisson process. Hence, the 
class of PRBs does not coincide with the class of mixed Poisson processes. 

The simulation of sample paths of a PRB is straightforward using the order statistics 
property of the jump times of the process. We shall provide example simulations of a 
PRB with its conditional expected terminal value and intensity. 

Using a construction similar to that of the compound Poisson process, the com- 
pound Poisson random bridge is a PRB with random jump sizes. We shall derive an 
explicit expression for the marginal characteristic function of this process. Although 
not explored in the present work, these compound processes could be applied to prob- 
lems in finance and insurance (as an alternative to the compound Poisson process). 

At the end of this chapter, we provide an application of PRBs to finance. We 
consider the problem of pricing an nth-to-default credit swap. The buyer of such a 
contract pays a continuous premium in exchange for a lump-sum payment on the date 
of the nth default from a basket of credit risky assets. Modelling the default times of 
the credit risky assets as the jump times of a PRB we are able to provide an explicit 
expression for the value of the swap, with the additional simplifying assumptions that 
(a) the interest rate is constant, and (b) premiums are paid continuously. 



9.1 Poisson random bridge 

Let {N tT } be an LRB with law LRB v ([0,T], {Q t },P), where 

e~ xt (Xt) k 

Qt(k) = k E No, (9.1) 

for some constant A > 0. We call {N^x} a Poisson random bridge. One can con- 
sider {N t r} to be a Poisson process conditioned to have the probability mass function 
P : No — > K+ at time T. We assume that 

E[iV TT ] = kP ( k ) < 00 • ( 9 - 2 ) 

k 

We can use the results from Section [2. Ill and Chapter [3] to derive various properties 
of {N tT }. It follows from ( I3.27P that the N sT - conditional probability mass function of 
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Ntt is 



Corollary 13.8.21 gives 



ik-N aT y. i 1 t) p ( k ) 



P s (k) = ±1 _ . (9.3) 

Z^k=N sT (k-N sT )\ \ L TV r \^> 



E [N tT | N sT ] = ?—^Nst + [Ntt \N sT }. (9.4) 

If E[iV|. T ] < oo, then conditioning on the terminal value Ntt an d recalling equation 
(12TT58D gives 



E[N? t ]=E[E[N? t \N tt ]] 

^(l-^E^ + ^yE^]. (9.5) 



(l-L) FjNrrrp] + ( - ) 

T 

This result can be generalised to the following: 
Proposition 9.1.1. If^2 k k 2 P(k) < oo, then 



+ yjly , Z [^t I iV sT ] + f ^ J E [iV| T I iV sT ] . (9.6) 

Proof. Fix s < T and define the process {??iT}s<t<T by rjtT = — -/VsT- It follows 
from the dynamic consistency property that, given N s t, {ijrr} is an LRB with terminal 
probability mass function P*(i) = P s (i + N s t)- With the understanding that E S [A] = 
E[A I N sT ], we have 

E s [N? T ] =E S [(N sT + Vo>) 2 ] 

= N 2 sT + 2N sT E s [r] tT ] + E s [rf T } 



N* T + 2N sT ^E s [ VTT ] + £±^E s [riTT) + (^) 2 E s [r7 2 



TTi 



N* T + 2N sT ^ s (E S [N TT ] - N sT ) + (E s [iV TT ] - N sT ) 

+ {^f [E S [N TT ] - 2N sT E s [N TT ] + N* T ) 
(^) 2 ^T + |^^(2E,[iV TT ]-l)iV sT 

+ T^^ANtt] + (|^) 2 E s [iV| r ]. (9.7) 

□ 
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Counting processes are characterised by the distribution of their jump times. We 
shall see that the distribution function of the ith jump of the PRB is an infinite sum 
of weighted beta probabilities. This is a consequence of the order statistics property 
of the jump times of Poisson processes. From this result it becomes apparent that 
the distribution of the next-jump time of a PRB can be expressed in terms of the 
conditional probability generating function of its terminal value. We shall then use the 
next-jump time distribution to write the intensity of the PRB in a form that highlights 
the impact of new information. 

Denote the probability generating function of P by Gp(z), i.e. 



G P (z) = E[z Ntt ] = J2 z k P(k). 



(9i 



k=0 



Note that Gp(z) < oo for < z < 1, so the definition of Gp is non-degenerate. Indeed, 
we have 

gP(o) 



P{k) 



k\ 



(9.9) 



Furthermore, since we have assumed that E[iVV r ] < oo, it follows that G' P (z) exists for 
< z < 1, with 



G' P (1) = lim G' P (z) = J"kP(k) = E[N TT ] 



(9.10) 



k=0 



The probability generating function of N t r is 



E [z NtT ] = E [E [z NtT | N TT \] 



E H-L + ± z 
GpH-^+'-z 



Ntt 



(9.11) 



for z < 1. Define the zth jump time of {Ntr} by 



T i = mf{te[0,T]:N tT = i}. 



(9.12) 



We adopt the convention that inf = oo; if Ntt < i then Tj = oo. From equation 
f )2.163p . the distribution function of the ith jump time Tj can be written in terms of 



126 



9.1 Poisson random bridge 



the regularized incomplete beta function I z [a,b] in the form 

oo 

Q [Ti < t] = [Ti < x | N TT = k] P(k) 

k=0 

oo 

= ^2 1 {i<i<k}h/T[i, k-i + l)P{k), 

for < t < T. We also have 

Q[Ti = oo] 



k=0 



(9.13) 



TT 



< I 



(9.14) 



k=l 



The distribution of Tj is mixed in the sense that it has a point mass at oo given by 
fl9.14p . and a continuous part with density 



, k—i 



00 (±y- L (i-±y 

k=0 1 



Note that 



I z [l,k} = l-(l-z) k , 



from which it follows that 

oo 

Q[Ti <t] = J2 



k=l 



P(k) 



fc=0 



P(fc) 



1-Gp( 1-- 



for i G [0,T]. The distribution function of 7\ is given by 



FT 1 {t) = < 



t < 0, 



1-Gp(1-±) 0<t<T, 



1 - P(0) 
1 



T < t < oo, 
t = oo. 



The density of the continuous part of the distribution is 



t \-> 1{0<*<T}^ G' P f 1 - — V 



(9.15) 



(9.16) 



(9.17) 



(9.18) 



(9.19) 
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and there is a point mass -P(O) at t = oo corresponding to the case where there are no 
jumps. For t G [0,T], we have 

tr = 0] = 1 - Q[Ti < t] = G P (l - I). (9.20) 



T , 

Since {A^t} has stationary increments, we also have 

Q[N u+t>T - N UtT = 0]=G P (l-^j, (9.21) 

where t, u satisfy 0<u<u + t<T. We can use the dynamic consistency property 
to find the conditional distribution of the waiting time until the next jump. We define 
the process {r] tT }s<t<T by rj tT = N tT - N sT . Then, given N sT , {q tT } is a PRB with 
terminal probability mass function 

P*(k) = P s (k + N sT ). (9.22) 

The probability generating function of P* is 

oo 

G P *(z) = J2 zkp *( k ) 

oo 



k=0 

oo 



= z n ~ NsT Ps(n) 

n=N sT 

= z- N ^G Ps (z). (9.23) 
Denoting the ith jump time of {r] tT } by T^ v \ we have 
Q[W i+1 < 1 1 N sT = i] = Q[Tj v) < t] 

= i - G P .m) 

= ±-mr NsT G P m), ^ 

for t G [s,Tj. This is equivalent to 

Q[iW - N uT > | N sT ] = 1 - (2^±a) ~ NsT G Ps (^±) , (9.25) 
for < s < u < u + 1 < T. 
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Proposition 9.1.2. The intensity of the PRB is 

E[iV TT | N sT ] - N sT 



Proof. 



A, 



. Q[N s+tiT - N sT — 1 1 N sT ] 



T-s 



(9.26) 



1 - 



lim 

t-*o+ 



— lim 
d 

du 



T-(s+t) 
T-s 



i (T-u\~ N sT n (T-u\ 



G' Ps (l)-N sT G Ps (l) 

T-s 
E[N TT | N sT ] - N s t 

T-s 



(9.27) 
□ 



Remark 9.1.3. The intensity of a counting process is related to survival probabilities 
of jump times by 



Q[T 1+NsT > 1 1 N sT \ = E 
where < s < t < T. 



exp ( — / \ u _ du 



N. 



sT 



(9.28) 



9.1.1 Example: Negative binomial prior 



Consider the case where P is a negative binomial probability mass function, i.e. 

'k + m — 1 



P(k) 



k 



p m (l-p)\ for k e N , 



(9.29) 



where m > and < p < 1 are constants. For m E N + , the negative binomial distri- 
bution gives probabilities relating to the number of failures of independent Bernoulli 
trials. If {Bi}^ is a sequence of outcomes from independent Bernoulli trials, each 
with probability of success p, then 

P(k) = Q[B k+m is the mth failure] (9.30) 
= Q[#{Sj : Bi = 0,1 < i < k + m} = m] . (9.31) 

The probability generating function for the negative binomial distribution is 



G P (z) = 



P 



1 - (l-p)z_ 



(9.32) 
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For the remainder of this chapter, if x\ is written and x ^ No, then we adopt the 
convention that x\ = T[x + 1]. 

Let {N t jn} be a PRB where Ntt has the probability mass function P given in (19.2 
Then we have 

Qr-tik-j) 



Qr(k) 

l{0<j<fc} 



-P(k) 



and 



XPTJ(k 



f^l (1 _ ±) k -> (k + m-l\ _ k 

j(k-j)\ [ T > V fc / ' 



(9.33) 



fc=— oo 



(m-1)! 
P 



'(i-p^EKi-lOa-*)]' 

fc=j 



_j fk + m — 1 
k-j 



e tx j\ 



p+Ui-p) 



i — p 



p + |(i -p) 



j + m — 1 



n=0 



n + j + m — 1 



e Uj| 



P 



p+|(l-p) 



1 — p 



p+|(l-p) 



j + m — 1 



(9.34) 



Hence the transition probabilities of {Nt?} are given by 

Q[iv tT = j | iV sT = i] = ^kMtiA Qt _ s (j - i) 



1 



{0<i<j} 



+ m - 1^ 


rp+i(i 


-p)" 




r ^(i-p) " 






Lp + |(i 


-p)_ 




_p + |(i-p)_ 





(9.35) 



and 

Q[iV Tr = j | N sT = i] 



1 



{0<i<j} 



E fc <MM) 
j + m — 1 



[p + f(l-p)] J+m [(l-f)(l-p 



(9.36) 



Remarkably, all the increments of the process {N t T}o<t<T have negative binomial dis- 
tributions. Note that 



Q[iV tT = o] 

as expected. By definition, we have 



P 



LP + U 1 -P) 

Ps(j) = Q[N t T=j\N sT 



Gp(l — ^) , 



(9.37) 



(9.38) 
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and so 



GpAz) 



oo 

E 

k=N sT 



Z k P s {k) 



JV„' 



Y^z k P s (k + N sT ) 



k=0 



-i m+N sT 



l-(l-p)(l-±) z 

It is then straightforward to verify that 



®[N t 



tr 



N 



sT 



I N. 



N sT 



-i m+AT sT 



T-t 

T- s 



T- S/ 

as expected. At time s, the conditional expectation of N?t is 



E[AT rr | N sT ] = N s t + (N sT + m) 



p+UI-p) 



7V sT + m(l-f)(l-p) 



and the intensity of the process {N t r} is 



E[iV TT | N s t] - N sT (N sT + m)(l-p) 



(9.39) 



(9.40) 



(9.41) 



(9.42) 



T — s pT + s(l — p) 

We see that a jump in the process {N t r} causes a jump in its intensity process. 
Biihlmann 2l|, 2.2.4] proposes a positive contagion model where the intensity of a 
Poisson process is a linear function of its current level, whic h bears som e similarity 
with the negative binomial PRB. One stark difference is that iBuhlmannl assumes the 
linear form for the intensity process in an ad hoc manner, whereas the intensity of the 
PRB is a consequence of Bayesian updating. 



9.1.2 Example: Log-series prior 

We shall see that a PRB with log-series terminal distribution has similar dynamics to 
the a PRB with a negative binomial terminal distribution. Notably, after the first jump 
(which occurs with probability 1) the increment distributions of the PRB are negative 
binomial. 
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For p G (0, 1), log-series distribution has probability mass function 

-1 (1 -p) k 



P(k) 



k G N+. 



(9.43) 



logp k 

Assume that {N t r} is a PRB with the terminal probability mass function given in 
([S32D - Note that N TT > 1. When N sT = 0, the transition probabilities of {N tT } are 
given by 

p)] 



( log[p + |(l 



tT 



j | iV sT = 0] = < 



log[p + #(l -p)] 



(J = o), 



[ jlog[p+ f (1 



p)] 



f — s / 



l-p) 



and 



Q[iVrr = j I iV sr = 0] 



P+ ?(l-p) 



(9.44) 



iiog[p + f (l-p)] 

When iV S T > the transition probabilities are similar to negative binomial case: 



(9.45) 



Q[N tT =j\N ; 



sT 



Q[N : 



TT 



3 I N, 



sT 



j-i 
3-1 
3 -i 



p + £(l-_p) 
p) 



P 



_T 

U 1 



¥(i-p) 



.p+H 1 



p) 



[ P + i(i- p )]'[(i _±) (i-p 



(9.46) 



(9.47) 



for i G N+ and j > i. Comparing (IQ5]) and (KT37I1 to (I9~3"5j) and (19T36D we see that 
{N t j<} can be considered a PRB with negative binomial terminal distribution in the 
limiting case m — > 0. 

The probability generating function of the log-series distribution is 

log[l-(l-p)z] 



Gp(z) 



\ogp 



The distribution function of the first jump time of {iV^} is 



Q[Ti <t] = l 



{t>0} 



1 



{0<t<T}" 



log[p+ |(1 -p)] 



logp 



Then 7\ has a continuous distribution with density 

— (1 — p) 

{0<t<T} _|_ — p)) logp 

Given iV S T = 0, the distribution function of T\ updates to 



Q[Ti < t 1 N sT = 0] = l {t > s} - l {s <t< T} - 



Iog[l-^(l-p)(l-i)] 



(9.48) 



(9.49) 



(9.50) 



(9.51) 



log[p+ | (1 -p)] 

If iVs^ > is given, then the jump time distributions are recovered from the negative 
binomial case by setting m = 0. 
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9.2 Mixed Poisson processes 

A wide and tractable class of PRBs consists of the mixed Poisson processes. Grendell 



46| provides a thorough exposition of these processes (see also Lundberg [63]). If every 
increment of a process has a mixed Poisson distribution (i.e., a Poisson distribution 
with a mixed mean) then it is a mixed Poisson process. We shall see that when the 
terminal distribution of a PRB is a mixed Poisson distribution, then the PRB is a 
mixed Poisson process. 

Let {Ntt} be a PRB with terminal probability mass function P(k) given by 

rpk POO 

P{k)= e k e- eT n{6)de, (9.52) 

ki Jo 

for some probability density 7r : IR + — > IR + : We say that Ntt has a mixed Poisson 
distribution. If is a random variable with density ir(9), then Ntt can be interpreted as 
a Poisson random variable with the unknown parameter BT. Without loss of generality, 
we my assume that A = 1 in ( 19. ip . so 

t k e- 1 



Then we have 



Qt(k) = - w - (9.53) 



Q T -t{k - j)P{k) 
Mk;j) = ^ (9.54) 

' 1 9 k e- eT n(0)d9, (9.55) 



o 



(k-j)\ 

for j, k G Nq satisfying j < k. It is straightforward to show 

VVfe') =e* / 9 j e- et 7t(9) d9. (9.56) 
Then the transition probabilities of {A^t} are given by 



/r = J I ^sT = l] = T /, ., Qt-s(j 



J °° 9ie- et ir(9) d9 (t - *VM P -(*-*) 



(9.57) 



f™9 i e- 6s Tr{9)d9 (j - i)\ 
where i, j G N satisfy z < j, and s, t satisfy < s < t < T; and 

TT = k\N sT = i} = ^t s{k '% - (9.5f 



f™e k e- eT ir(6)de (T-s] 



k-j 



/ °° ffie-0"ic(0) d9 (k-j)\ 



(9.59) 
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where i, k G No satisfy % < k, and < s < T. When n(6) = 5^(9) for some [l > 0, 
the transition probabilities reduce to those of a Poisson process with parameter fi. So, 
given the value of G > 0, we have 

Q[N tT = j | N sT = i, 6] = 11 !>_.* , (9.60) 

for i < j and s < t < T. 

We shall show that the transition distributions of {Ntr} are mixed Poisson distri- 
butions. First, we find the conditional distribution of G given the information to date. 
For < ti < t 2 < ■ ■ ■ < t n < T, the joint law of ({N tuT }i, N TT , Q) is given by 



Q [N tliT = h, N t2tT = h, . . . , N tniT = k n , N T = K,Qe d6] = 

QT-t n (K - fen) nr=i Qu-U-Ah - fcj-i) (0T) K e- 9T 
Qt(K) K\ 

where to = and ko = 0. Using the Bayes theorem, we have 



n(9)d9, (9.61) 



Q [N TT = K,Q e d9\ N tl , T = h, N t2>T = k 2 ,..., N tn , T = k r , 

'-tt(6) d6 



} T -t n {K-k n ) (eT) K e- » T 
Q T (K) K\ 



Q T -t n {K -k n )9 K e- eT 7r(9) d6 
Jo°° T,K= kn Qt-u(K - k n )9^n(9) W 



(9.62) 



This implies that 



O r Q r A9 I t n } YZLsrr ~ N tT )9 k e-^{9) d9 

^ L ' 1 J r EZ NtT Qr-t n (k - N tT )9*e-^(9) d9 

9 N < T e~ et n(9) d9 



/°° 9 N ^e-^7r(9) d9'' 



(9.63) 



10 

where {J 7 ^} is the filtration generated by {N tT }. We write 

9 NtT e- et n(9) 

= jr^e-»4)d<r (9 - 64) 

thus 7Tt(0) is the Ff- conditional density of 9. From equation ( J9,60p . we then have 

Q[N tT = j | N sT , = 9} Q[6 G d9 | N sT ] 

=o 

= [ . I r r. / ^-e-^Tr.^d^. (9.65) 
U - N sT) ] - Jo 

Hence, all the increments of {N tT } have a mixed Poisson distribution. 
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Remark 9.2.1. The negative binomial PRB encountered in the previous section can 
be recovered by taking 

1/ m \ Tfl / rji 

I P T \ nm-l _____ I P T 



In this case 9 is a gamma random variable. 

Remark 9.2.2. The log-series distribution is not a mixed Poisson distribution (note 
that -P(O) = for the log-series distribution). Hence, if the terminal distribution of a 
PRB is log-series then the PRB is not a mixed Poisson process. Indeed, if P(k) = 
for any k G No, then a PRB with terminal mass function P is not a mixed Poisson 
process. 

The iV S T-conditional distribution of Ntt is 

°° (rp _ \k~N sT poo 

nN TT | n sT ] = J2 k (k _ N )\ / e k - NsT e- e{T - s) -Ks(e)de 

k=N.~ ^ sT '' ^° 



sT 



6 7T s (e) d6 

= E[Q\N sT ]. (9.67) 
Then the intensity of {N t -r} is given by 

E\e\N sT ]-N sT _ i fj o °°e^e-^(0)de _ N \ 



T-s T - s V J °° 9 N ^ e - es 7i(9) d9 

9.3 Simulation 

An efficient method for the simulation of a PRB is to generate the jump times. The 
case Ntt — is trivial since we have N t T = for all t e [0, T]. When Ntt > we set 
the jump times to be given by 



E1+Ntt TP 
1=1 E 3 



Ti = JTIZL - for 1 < i < Ntt, (9.69) 



where the E^s are independent, identically distributed exponential random variables 
(with arbitrary parameter). See Figure I9TT1 for example simulations. 
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Figure 9.1: Top row: Two simulations of a Poisson random bridge with a negative 
binomial terminal distribution are plotted (bottom line) with the conditional expecta- 
tion of the terminal value (top line). Bottom row: The intensity process of the Poisson 
random bridge is plotted. The negative-binomial parameters are m = 8, p = 0.4. 

9.4 Compound Poisson random bridge 

Compound Poisson processes are Levy processes. If the jump size distribution of a 
compound Poisson process is discrete then associated discrete LRBs can be defined, and 
if the jump size distribution of a compound Poisson process is sufficiently smooth then 
associated continuous LRBs can be defined. We call such LRBs compound-Poisson 
random bridges (C-PRBs). Hence, a C-PRB can be considered to be a compound 
Poisson process conditioned to have a particular marginal distribution at time T. The 
a priori distribution of the number of jumps and the distribution of the jump sizes 
depend on the choice of terminal distribution. The jumps of the C-PRBs occur at the 
jump times of a PRB that is adapted to the filtration of the C-PRB, and in general 
the jump sizes of a C-PRB are not independent. Progress can be made in the analysis 
of C-PRBs, but dependence between the jump-size and jump-time distributions makes 
the path tricky. 
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9.4 Compound Poisson random bridge 



A simpler but more tractable class of processes can be constructed by allowing the 
jump sizes of a PRB to be random and independent. We call such processes compound 
Poisson random bridges (CPRBs), and trust that the absence of a hyphen is enough 
to avoid confusion with C-PRBs. It is CPRBs that we consider for the rest of this 
section. 

Let the process {Ntr}o<t<T be a PRB, and let Ntt have arbitrary marginal prob- 
ability mass function P <C Qt- Let {A^} be a sequence of independent, identically 
distributed random variables with characteristic function xx- We assume that the 
XiS are independent of {N tT }. We then define the CPRB {Y tT } by 



Y 



tT 



5> 



(9.70) 



i=i 



If Xi 7^ then {Y tT } jumps at the same times as {N tT }, but the jump sizes of {Y tT } 
are random. The characteristic function of Y t r is 



E [e iaYtT ] = E 



E 



N tT 

oxp } in Xj 

i=i 

N tT 

l[E[exp(iaX t )\N tT } 

.i=i 



Using the dynamic consistency property, this can be generalised to 



(9.71) 



E [( 



Y sT ,N sT ] =e k 



= e 



E 



E 



exp la 



i=N sT +l 



N, 



sT 



N tT 



II E[exp(iaX,)|iV sT ,iV tr ] 

-N sT +l 



e iaY ^E[ X x(a) 



e iaY ° T G P * ( 1 - 



NtT-N* 



t-s 



+ 



N sT ] 

t - i 



jctY s T 



T-s T-s 

J t — s 

T-s r T-s 



Xx(a) 



N„ 



(9.72) 
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Counting processes compounded by random jump sizes are common in the insurance 
literature, and we refer the reader to Buhlmann 
for further details and references. 



211 ] . Lundberg [63|, and Mikosch 



9.5 nth-to-default baskets 

We now consider the pricing of an nth-to-default credit swap. The swap provides 
protection against defaults occuring in a basket of K homogeneous credit risks. In 
return for a continuously paid premium, the buyer receives an amount 1 — R on the date 
of the nth default, if this default occurs before date T. The buyer pays the premium 
until the earlier time of the nth default and T. Here R represents the recovery rate on 
the credit risks, and is for simplicity assumed constant. We also assume that interest 
rates are constant, i.e. r t = r, for all t. 

We assume that the defaults occur on the jump times of the PRB {N t r} with 
terminal probability mass function P : {1, 2, . . . , K} — > R + . Denoting the nth jump of 
{N tT } by T n , the value of the swap to the buyer is 



V = E 



»T n AT 

(9.73) 



/ e- rt e 9i dt+(l-i?)l {Tn < T} e- rT " 
Jo 



where q is the premium rate. This can be rewritten as 

V = -^(E [l {Tn <T } e {q ~ r)Tn } + e(9_r)T ® l T n > T] - 1) + (1 - R) E [t { T n <T}^ rTn ] • 

(9.74) 

Note that 

n-l 

Q [T n > T] = Q [N TT <n] = J2 P(Q. (9.75) 

k=l 

Recalling the expression for density of T n given in f l9.15p . we have 

/ TB[n,k-n+l) 



rT ( t_\ f-\ t_\ 



k=n 
K 



= ^M[n,k + l,-rT]P(k), (9.76) 

k=n 

where Kummer's function M[a,(3,z] was given in f l2.53p . The initial value of the 
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contract is then 



ra-1 



V = — !— I J2 M K k + !> (? - r)T]P(A;) + e (9 " r)T ^ " 1 

^ r \fc=n fc=0 / 

K 

+ (1-R) ^M[n,k + l,-rT]P{k). (9.77) 

k=n 

Typically, the premium rate q would be set to ensure that Vq = 0, so the swap has zero 
initial value. At time s < T, if the nth default has not yet occurred, the value of the 
swap to the buyer is 

*T»AT 



V S = E 



e -^ d t + (l- J R)l {Tn < T} e- rT " 



AT, 



sT 



(9.78) 



By the dynamic consistency property, if we know N s t we can define an LRB by 



r ltT = N tT -N sT (s<t<T). 



(9.79) 



Then we can write 



V a 



q — r 



-E 



1 {T^ <T}' 



iV. 



sT 



g — r 
+ (1 -P)E 



[T n >T|iV sT ] 



q — r 



-rT 



(-7) 



AT, 



sT 



(9.80) 



In a similar way to the calculation of Vo, we find 

K 

(j - r 



V s = -^— r J2 M[n - N sT , k-N sT + 1, (g - r)(T - S )]P.(A;) 

n— 1 -. 

£ - — 



k=n 
1 



,(?-r)T 



— r * — ' ' ' q — r 

* k=N aT H 

K 

+ (1-R)J2 M[n - N sT , k-N sT + 1, -r(T - s)]P 8 (k), 



(9.81) 



k=n 



where P, is the ^^-conditional probability mass function of Ntt, and is given by ( 19. 3p . 



Appendix A 
Some integrals 



A.l Proof of Proposition 12.8.2 



Proposition. For y G [0,z], the distribution function of the random variable Sj>j) is 
given by 



Ftr(y; z) = $ 



c(Ty - tz) 



^yz(z-y) 
Proof. We need to show that 



... f l_p) ( -.2-'l7-/, : ([) 



c{(2t-T)y-tz) 
^yz(z-y) 



(A.l) 



1 Ct{T-t) [ yeX P{-2 uz(z-u) 



1 c 2 (Tu-tzf 



V2n T J (u-u 2 /zf 2 
c{Ty - tz) 



du 



y/yz(z - y) 



c{{2t-T)y-tz) 
^yz(z-y) 



First we define 



x\u) 



c(Tu — tz) 



■ (A.2) 



(A.3) 



^Juz{z — u) 

Inverting this function is a matter of finding the positive root of a quadratic equation, 
and it gives 

:{2cHT + x 2 z) + xz^Ac 2 t{T -t)z + x 2 z 2 



u(x) 

Differentiating x{u) yields 

x'(u) 



2(c 2 T 2 + x 2 z) 
(T -2t)u + tz 



2: J {u-u 2 /zf/ 2 ' 
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(A.4) 



(A.5) 
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A. 2 Proof of Proposition I21T51 



Writing 



we have 



a(x) = yJk?t{T -t)z + x 2 z 2 , 



(A.6) 



t(T - t) 



2z 



T (T — 2t)u(x) + tz 

At{T-t){c 2 T 2 + x 2 z) 



Ac 2 tT 2 (T -t)+ T 2 x 2 z + T(T - 2t)xa(x) 

T 2 a{x) 2 — (T — 2t) 2 x 2 z 2 
T 2 a(x) 2 + T(T -2t)xza(x) 
(Ta(i) + (T — 2t)xz){Ta{x) — (T — 2t)xz) 

Ta{x){Ta{x) + (T — 2t)xz) 
! (T - 2t)xz 
Ta(x) 



Then (1A.5j) and (1A.7|) give 



(T - 2t)x{u)z 



x'(u) 



ct{T - t) 



T^Ac 2 t{T - t)z + x(u) 2 z 2 J " ^ T(m-m 2 /2;) 3 / 2 ' 
So making the change of variable x = x(u) on the left hand side of (IA.2[) gives 

(T - 2t)xz 



(A.7) 



(A.8) 



27T 

$[x(y)] 



1 - 



Ty/Ac 2 t(T - t)z + x~ : 



2t 
T 

2t 



2^2 
X 



dx 



-x 2 /2 



^Ac 2 t(T -t)/z + x 2 V2n 



dx 



$[x(y)] - 1 1 - - ] e 2c2 *( T -*)^ ($ [^4cH(T-t)/z + x(y) 2 ] - l) 



2f 



$[x(y)] + 1 - - e 2c t{T - t)/z $ -^^(T-^/^ + x^) 2 



c(Ty - tz) 
^yz(z-y) 



2t 



+ (!___] e 2c 2 t(T-t)/z $ 



c({2t-T)y-tz) 
^yz(z-y) 



(A.9) 
□ 



A. 2 Proof of Proposition 12.8.5 



Proposition. Define the incomplete first moment of S" tT by 



M tT {y;z) 



uf tT (u;z)du (0<y<z). 



(A.10) 



A. 2 Proof of Proposition [233] 
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Then we have 



M tT (y; z) = -z{§ 



c(Ty - tz) 



and the second moment of is given by 



_ e 2c?t(T-t)/z $ 



c((2t-T)y-tz) 
y/yz(z-y) 



, (A.ll) 



E 



S, 



^z 2 ll- c(T - t)e^J^ $ 1-cTz- 1 ' 2 ] 



Proof. Making the change of variable f]A.3[) . we have 

t(T-t) [ x{y) 2zu(x) 



uftfr(u;z) du 



T 



-x 2 /2 



(T - 2t)u(x) + tz 



dx. 



Recalling equation (IA.7j) . we have 
t(T-t) 2zu{x) 



t{T-t)z 

~f (T - 2t)u(x) + tz ~ T(T - 2t) 

tz 



2 - 



2tz 



T{T - 2t) 



(T - 2t)u(x) + tz 
(T - 2t)xz 



T -2t + 



a(x) 



t 



1 + 



xz 



a(x) 



Then we have 



u ftr(u; z) du 



t(T - t) 
T 



2zu(x) 



-x 2 /2 



dx 



, (T - 2t)u(x) + tz y/fr 

z<&[x(y)]+ —== 
' 1 ^/4c 2 t(T -t)/z + x 2 V^F 



T 



x 2 /2 



dx 



t 



--z \ $[x(y)] - e 2c t ^ T ~^ z $ -v/4c 2 t(T -t)/z + x{y) 



T 



ciTy - tz) 
sjyzjz - y) 



e 2c 2 t(T-t)/z $ 



c((2t-T)y-tz) 
\Jyz{z - y) 



as required. 

Changing the variable in the second integral of the proposition gives 



Z 2r ( \ j t(T-t) 

u ftr(u; z) du = — - — 

1 x 



2zu(x) 2 e x2 l 2 



(T - 2t)u(x) + tz 



dx. 



(A.12) 



(A.13) 



(A.14) 



(A.15) 



(A.16) 
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A. 3 Proof of Proposition I21TT1 



Now, using (1A.6j) and ( 1A.14j) . we have 

t(T - t) 2zu(x) 2 
T (T - 2t)u{x) + tz 



T 

t 
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Proof. Writing a = ct and b = c(T — t), we need to calculate the integral 
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The distribution function is then 
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Substituting for a and b yields the required result. 
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Then we have 
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